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INTRODUCTION 


The  design  of  a helicopter  drive  system  presents  a difficult 
challenge.  The  continually  escalating  requirement  for  improved 
power -to -weight  ratio  capability  is  further  complicated  by  the 
additional  demands  for  extended  service  life,  improved  relia- 
bility/maintainability, better  survivability/vulnerability  and 
reduced  vibration/noise  levels.  Faced  with  such  an  array  of 
design  constraints,  the  helicopter  drive  system  design  team 
must  consistently  push  technology  beyond  the  state  of  the  art. 

In  addition  to  these  technological  advancements,  a major  goal 
for  helicopter  improvement  is  to  reduce  operating  cost.  Since 
30  percent  of  the  total  helicopter  direct  maintenance  cost 
is  associated  with  the  drive  system,  drive  system  cost  is  a 
significant  part  of  the  total  operating  cost  (Figure  1) . 

Service  life  duration  of  the  major  transmission  components  is  a 
prime  contributor  to  the  operating  cost.  The  most  effective 
way  to  reduce  the  drive  system  cost  is  to  extend  the  average 
service  life  of  major  transmission  components  (gears,  bearings, 
splines,  retention  hardware,  interface  connections  and  joints, 
and  lubrication  system  components) . Increased  life  will  reduce 
spares  procurement,  maintenance  requirements,  depot  facilities, 
and  labor. 

Considerable  research  has  been  devoted  to  investigating  and 
improving  individual  transmission  components  such  as  gears, 
bearings,  and  lubrication  systems.  Conversely,  helicopter 
transmission  housings  had  not  previously  received  the  attention 
necessary  to  define  fully  and  optimize  their  functional  require- 
ments, and  a gap  in  transmission  technology  existed  in  this 
area.  Complex  structures  such  as  rotor  transmission  housings 
used  in  current  helicopters  are  not  generally  designed  for 
stiffness  characteristics,  however  they  have  high  strength  mar- 
gins for  safe  life  and  seldom  exhibit  gross  structural  failures. 
However,  since  the  housing  provides  structural  support  to  the 
internal  components,  its  physical  characteristics  significantly 
affect  overall  transmission  performance  and  life  in  terms  of 
internal  bearing  capacity,  gear  capacity,  fretting,  misalign- 
ments, and  load  maldistributions.  Therefore,  the  full  benefits 
of  advancements  achieved  in  component  technology  cannot  be 
realized  in  practice  until  the  housing  has  been  optimized. 

A major  constraint  in  meeting  cost  improvement  goals  is  the 
deflections  of  the  present  housings  under  operating  load 
conditions.  Housing  deflections  under  load  have  been 
identified  as  a cause  of  accelerated  wear  and  surface  deteri- 
oration of  gears,  bearings,  splines,  retention  hardware,  and 
interface  connections  and  joints.  In  addition,  these 
deflections  result  in  excessive  gear  misalignment  which 
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greatly  increases  the  dynamic  force  (excitation)  at  the  gear 
meshes,  resulting  in  vibration  and  noise  generation.  Reduction 
in  the  magnitude  of  these  housing  deflections  by  structural 
optimization  and/or  advanced  materials  offers  great  promise 
for  both  prolonging  the  life  of  the  transmission  components 
and  substantially  reducing  the  vibration  and  noise  level 
without  a system  weight  penalty. 

To  continue  to  improve  transmission  analysis  capability,  a 
detailed  understanding  of  the  structural  and  thermal  aspects 
of  the  transmission  housing  must  be  developed  and  integrated 
with  the  existing  component  analyses.  This  effort  for  further 
research  in  the  stress  and  dynamic  characteristics  of  a heli- 
copter transmission  housing  uses  the  finite  element  technique. 
The  effects  of  thermal  growth,  structural  deformation  and 
operational  housing  life  can  be  accurately  studied.  Along 
with  structural  design,  fabrication  and/or  modification,  and 
testing,  this  program  also  attempts  to  analyze  the  housing  for 
fail-safe/safe-life  design. 


BACKGROUND 


To  provide  an  understanding  of  the  configuration,  functional 
requirements,  and  design  criteria  for  a helicopter  transmission, 
a brief  description  is,  included  here.  A typical  contemporary 
helicopter  main  transmission  housing  is  composed  of  three  major 
parts  with  essentially  separate  functions:  the  upper  cover, 
the  ring  gear,  and  the  case.  An  example  of  this  configuration, 
the  CH-47C  forward  rotor  transmission  is  shown  in  Figure  2. 

The  upper  cover  supports  the  rotor  shaft  and  provides  lugs  for 
mounting  the  transmission  to  the  airframe.  Hence,  the  rotor 
system  loads  are  transmitted  through  the  upper  cover  into  the 
airframe.  The  case  contains  and  supports  the  main  bevel  gears 
and  may  also  include  a tail  rotor  or  sync  shaft  drive,  lube 
pump,  or  accessory  drives.  The  transmission  may  also  have  a 
separate  sump  for  containment  of  the  lubricant,  as  does  the 
CH-47C,  or  it  may  use  an  integrally  closed  lower  portion  of  the 
case  for  this  purpose.  The  stationary  ring  gear,  which  con- 
nects the  upper  cover  and  case,  contains  the  planetary  gear 
system.  The  entire  housing  also  performs  the  functions  of 
sealing  in  the  lubricant,  providing  passages  for  lubricant 
delivery,  protecting  critical  transmission  components,  and 
dissipating  heat  by  conduction  and  radiation.  Figure  3 shows 
the  transmission  case  in  detail,  since  much  of  the  work  herein 
is  concentrated  on  analysis  of  the  case. 

The  upper-cover  design  criteria,  in  order  of  importance, 
include  ultimate,  fatigue,  and  crash  load  conditions.  The 
gear  case  design  criteria  include  stiffness  for  gear  mounting 
and  fatigue  loads  in  certain  areas.  The  ring  gear  must 
provide  adequate  strength  to  react  the  planetary  gear  loads 
and  to  support  the  case  and  must  provide  sufficient  stiffness 
to  maintain  planet/ring  gear  tooth  alignment.  The  requirement 
for  oil  containment  exists  for  all  of  the  housing  parts. 

Because  of  the  complex  geometry  and  the  many  functions  that  the 
typical  housing  must  perform,  transmission  housings  have  tradi- 
tionally defied  analysis.  Due  to  the  lack  of  analytical 
methods  for  predicting  and  optimizing  its  many  functional  load 
paths,  transmission  housings  have  evolved  as  generally  thick- 
walled  cast  or  forged  structures  made  of  aluminum  or  magnesium. 
With  the  increasing  power,  size  and  weight  of  helicopter 
transmissions,  the  resulting  over-design  produces  undetermined, 
but  probably  substantial,  weight  penalties.  Figure  4, 
which  presents  a weight  breakdown  for  a typical  helicopter 
transmission,  shows  that  the  housing  weight  is  24%  of  the 
total  transmission  weight.  It  is  recognized  that  minimum 
wall  thickness  in  nonstructural  areas  is  limited  by  the 
casting  process.  Nevertheless  improvements  in  weight  can  be 
realized  in  the  thicker  structural  support  regions. 
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Figure  3.  CH-47C  Forward  Rotor  Transmission  Case  — a.  Left  Side  View 


Figure  3.  Continued  — b.  Right  Side  View 
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Figure  3.  Continued  — d.  Aft  View 
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Figure  3.  Continued  — f.  Bottom  View 
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Figure  4.  Example  of  Main  Transmission  Weight  Breakdown 


Considerably  more  significant  application  of  such  a structural 
load  path  analysis  can  be  made  to  advanced  transmission  concepts 
employing  fabricated  housings,  composite  materials,  and  other 
advanced  concepts,  which  will  permit  greater  design  flexibility. 
For  example,  major  load  paths  could  be  selectively  reinforced 
while  the  thickness  of  non-load  carrying  regions  of  the  housing 
wall  could  be  reduced  to  the  minimum  necessary  for  containment 
of  the  lubricant. 

A basic  requirement  of  a helicopter  transmission  is  dimensional 
stability  of  each  gear-mesh,  which  implies  dimensional 
stability  of  the  housing  and  bearing  mounting  locations. 

Because  of  the  influence  of  the  housing  stiffness  on  the 
transmission  internal  components,  it  must  be  optimized  during 
the  development  of  an  advanced  technology  transmission  system. 
Also,  initial  knowledge  of  the  dynamic  stability  of  the 
system  must  be  obtained  to  avoid  resonances  and  the  accompany- 
ing vibration/noise. 

Analytical  evaluation  of  the  load-carrying  capacity  of  bevel 
gears  involves  assumptions  regarding  the  nature  of  the  tooth 
bearing  for  the  specific  gear  mountings  under  load.  A uniform 
stress  distribution  across  the  tooth  and  rigid  mounting  is 
typically  assumed.  Unless  these  assumptions  are  accurate, 
actual  stresses  may  vary  considerably  from  the  calculated 
values,  resulting  in  a possible  life  reduction.  During  manu- 
facture of  bevel  gears,  the  desired  tooth  contact  pattern  is 
established  from  observation  of  the  pattern  obtained  under 
light  load  in  a gear  tooth  pattern  checker.  It  is  possible 
to  vary  the  length,  width,  and  position  of  a tooth  bearing 
by  selection  of  grinding  wheel  diameters  and  grinding  machine 
settings.  This  procedure  is  known  as  developing  the  tooth 
bearing  contact  pattern. 
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with  gear  mountings  that  are  rigid,  the  behavior  of  the  tooth 
bearing  under  load  is  usually  more  predictable,  and  thereby 
can  be  developed  using  previous  experience.  However,  in  the 
case  of  aircraft  applications,  the  mounting  designs  are 
markedly  different  in  that  rigidity  is  sacrificed  in  favor  of 
weight  reduction.  Therefore,  it  is  seldom  possible  to 
accurately  predict,  during  the  design  stage,  the  size  and 
position  of  the  tooth  contact  pattern  required  at  no  load  in 
order  to  obtain  the  desired  bearing  pattern  in  the  final  gear 
mountings.  A study  of  the  mounting  design  and  operating 
conditions  together  with  a judgement  based  on  experience  must 
be  utilized  to  establish  the  initial  tooth  bearing.  From 
this  point,  the  development  of  the  final  tooth  bearing  is 
accomplished  by  actual  trial  of  the  gears  in  their  final 
mountings . 

Predicted  improvements  in  the  load  capacity  for  gears  and 
bearings  may  be  offset  in  practice  by  poor  load  distribution 
resulting  from  misalignment  caused  by  the  deflection  of  mounting 
surfaces  within  the  housing  (Figure  5) . The  detrimental 
effects  of  misalignment  on  gear  teeth  is  documented  by  the  AGMA 
(Reference  1) . Gear  tooth  bending  and  surface  contact  stresses 
are  proportional  to  load  distribution  factors  (I^,  Cm) . These 
factors  evaluate  the  effects  of  nonuniform  load  distribution. 


Figure  5.  Measured  and  Theoretical  Stresses  in  Thin  Rimmed  Bevel  Gears 


1.  American  Gear  Manufacturers  Association  Standard  210.02. 
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They  are  dependent  upon  several  factors  including  gear  mesh 
misalignment  due  to  housing  distortion  caused  by  loads  and 
thermal  variations  (Figures  6 and  7) . The  effect  of 
different  rates  of  misalignment  is  shown  in  Figure  8,  F 
represents  the  face  width  having  100-percent  contact  for'^^a 
given  tangential  load  and  misalignment  error.  Uneven  load 
distribution  caused  by  significant  shaft  misalignment  will 
result  in  tooth  pitting  and  scuffing  failures.  Gear  mesh 
misalignment  is  also  important  from  the  aspect  of  vibration/ 
noise  generation  (Figure  9,  Reference  2). 

Shifting  of  the  gear  tooth  contact  pattern  may  also  be  caused 
by  differential  thermal  growth.  When  two  bevel  gears  are 
properly  mounted,  and  given  that  the  gears  have  been  properly 
manufactured  and  the  bearing/gear  positioning  shim  stack 
"heights"  are  properly  assembled,  the  cone  centers  are  coin- 
cident (at  room  temperature) . Since  the  gears  and  bearings 
are  made  of  steel  and  the  housing  and  bearing  cartridges  are 
made  of  magnesium,  differential  thermal  expansion  (thermal 
coefficient  of  mag  approximately  twice  that  of  steel)  would 
cause  the  relative  positions  of  the  bevel  pinion  and  gear  to 
change.  The  cone  centers  would  not  be  coincident  at  operating 
temperature  and  the  tooth  pattern  and  stress  distribution 
across  the  tooth  would  change. 

Housing  deflection  also  has  detrimental  effects  on  the  per- 
formance of  the  shaft  seals.  The  operational  effectiveness 
and  life  of  a helicopter  transmission  seal  is  dependent  upon 
a number  of  interdependent  factors  such  as  rubbing  velocity, 
pressure,  temperature,  dimensional  and  finish  conditions  of 
the  shaft  surface,  nature  of  the  sealed  medium,  housing 
deflections,  and  shaft  deflections.  All  of  these  factors 
become  proportionally  critical  in  high  speed  applications. 
"Leakers"  pose  hazards  including  fire,  personnel  injury  due 
to  slippery  surfaces,  and  objectionable  appearance.  Severe  seal 
leaks  can  result  in  depletion  of  the  oil  supply.  Studies  of 
advanced  helicopter  transmission  concepts  indicate  that  shaft 
speeds  will  be  increased  and  time  between  overhauls  for  trans- 
missions will  be  extended.  These  requirements  will  have  a di- 
rect effect  on  shaft  seal  design  by  imposing  more  stringent 
performance  requirements  and  extended  useful  life.  Reduction 
of  deflections  and  vibrations  will  help  achieve  this  goal. 

Experience  to  date  has  been  associated  with  reduced  bearing 
lives  due  to  varying  system  stiffness.  Present  bearing  life 
equations  assume  that  the  bearing  is  rigidly  supported, 
operates  under  no  misalignment,  and  operates  under  a constant 
and  uniform  load.  In  helicopter  applications,  these  assumptions 
are  not  true  and  therefore  calculated  lives  are  not  precise. 


2.  George  C. , GEAR  NOISE  SOURCES  AND  CONTROLS,  Detroit  Diesel 
Allison,  Division  of  General  Motors  Corporation. 
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PINION 


Figure  6.  Example  of  a Pinion  and  Gear  Misalignment  Under  No  Load 


(REFERENCED 

Figure  7.  Load  Distribution  Across  Face  Width  for  Various  Contact  Conditions 
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FORCE  - ERROR  COEFFICIENT  (C.) 


^"1  (REFERENCED 

Figure  8.  Spur  Gear  Load  Distribution  Factor  — 


Figure  9.  Influence  of  Tooth  Alignment  Error  on  Gear  Noise 


Early  bearing  failures  have  resulted  from  shaft  misalignment 
(edge  loading)  and  nonuniform  housing  support  (local  hard 
spots).  Present  methods  of  bearing  analysis,  which  include 
complex  computer  programs,  have  very  limited  capabilities  of 
evaluating  the  effect  of  structural  shapes  and  flexibility  on 
bearing  life.  An  example  of  this  type  of  problem  has  been  ex- 
perienced during  the  testing  of  a large  swashplate  bearing 
(Reference  3) . This  bearing  was  sized  by  the  best  available 
computer  programs  and  modified  in  an  attempt  to  account  for 
structure  deflection.  After  239  hours  of  testing,  two  bearings 
failed  due  to  distress  caused  by  the  ball  riding  over  the  edge 
of  the  inner  and  outer  races.  Investigation  of  the  failure 
revealed  that  relatively  large  structure  deflection  (0.016  inch 
ring  separation)  resulted  in  high  local  loading  and  excessive 
race  depth  requirements.  The  severe  effect  of  misalignment  on 
bearing  life  is  further  shown  in  Figure  10. 

In  order  to  reduce  this  type  of  failure,  special  consider- 
ation must  be  given  to  the  uniform,  rigid  support  of  critical 
components.  Reduced  shaft  and  housing  deflection  will  result 
in  better  performance  and  life  of  gears,  bearings,  and  seals. 
Therefore,  analytical  methods  must  be  developed  which  will 
permit  evaluation  of  important  design  paramete"o,  allow  trade- 
studies,  and  provide  guidance  to  designers.  The  specific 
topics  being  investigated  under  this  program  and  the  areas  of 
anticipated  benefit  are  summarized  in  Table  1. 


0 0.001  0.002  0.003  0.004  0.005  0.006  0 0.001  0.002  0.003  0.004  0.005  0.006 

TANGENT  OF  MISALIGNMENT  ANGLE  TANGENT  OF  MISALIGNMENT  ANGLE 

Figure  10.  Effect  of  Misalignment  on  Service  Life  (REFERENCE  4) 

3.  Lenski,  J.  W. , HLH  SWASHPLATE  FAILURE  ANALYSES,  Boeing 
Vertol  Inter-Office  Memorandum  8-7462-1-40,  dated  10  July 
1974. 

4.  Leibensperger,  R.  L.,  WHEN  SELECTING  A BEARING  LOOK  BEYOND 
CATALOG  RATINGS,  Machine  Design,  April  1975,  Pages  142-147 
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TABLE  1.  AREAS  OF  INVESTIGATION 
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ANALYTICAL/COMPUTER  METHODOLOGY 


For  the  work  conducted  under  this  contract,  the  finite 
element  technique  using  the  NASTRAN  computer  program  was 
selected  for  structurally  analyzing  the  housing  stresses, 
dynamic  response,  deflections,  and  thermal  effects,  as  well 
as  strength,  weight,  and  fail-safe/safe  life. 

The  development  of  this  approach  results  in  an  advanced 
tool  for  analyzing  a transmission  housing  through  the  use 
qf  finite  element  techniques.  Since  NASTRAN  is  used,  the 
results  indicate  the  feasibility  of  using  NASTRAN  as  a 
versatile  transmissions  design  tool,  specifically  in  the  areas 
of : 

1.  Thermal  distortions,  cooling  requirements. 

2.  Stress  Analysis 

A.  Maneuver  Loads 

B.  Crashworthiness  simulation  using  inertia  relief 
capability  of  NASTRAN  (Rigid  Format  2) . 

3.  Dynamic  analysis  due  to  n per  rev  hub  excitations. 

4.  Weight  reduction  through  optimization  schemes  such  as 
strain  energy  distribution  or  fully  stressed  design 
analysis. 

5.  Design  of  conceptual  transmissions. 

6.  Reduction  of  bearing  misalignment  by  developing  a 
stiff er  case  and  evaluating  geometry  and  slope  changes. 
Reduction  of  load  maldistribution  across  gear  teeth 
and  improved  gear  life  by  similar  means. 

7.  Consideration  of  f ail-safe,'3afe  life,  survivability  and 
vulnerability  by,  for  example,  simulating  a crack  in  the 
mode 1 . 

8.  Assess  radar  cross-section  by  evaluating  scale  plots  of 
the  transmission  drawn  at  any  specified  orientation. 

9.  Effects  of  nonlinear  behavior  of  bearings. 


The  contractor  has  utilized  pre-  and  post-processor 
computer  programs  compatible  with  NASTRAN  to  improve  its 
utility.  Examples  of  such  programs  are  SAIL  II  (an  input 
data  generator)  and  S-83  (strain  energy  analysis), 
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The  automatic  generation  of  grid  point  coordinates, 
element  connectivity,  etc.,  is  essential  for  complex 
models.  Boeing  has  developed  a sophisticated  finite 
element  input  capability  for  use  with  NASTRAN  entitled 
SAIL  II  (Structural  Analysis  Input  Language) . This 
preprocessor  allows  the  user  to  take  advantage  of  any 
pattern  that  occurs  in  the  data  by  making  available 
straight-forward  techniques  for  describing  algorithms  to 
generate  blocks  of  data.  Grid  points  and  element  connec- 
tions may  be  generated.  This  program,  although  proprietary 
to  the  contractor,  has  been  utilized,  and  is  available  for 
purchase  by  industry.  An  alternative  for  nodal  generation 
and/or  connectivity  would  be  user  generated  WATFOR  computer 
programs  which  would  punch-out  the  NASTRAN  input  bulk  data 
cards.  Although  the  contractor  has  used  SAIL  II  and  feels 
that  this  program  is  more  versatile,  other  users  of  NASTRAN 
can  conduct  the  same  work  by  alternate  methods.  The  Boeing 
Vertol  SAIL  II  computer  program  will  be  compatible  for  use 
with  NASTRAN  Level  16. 


For  ease  of  identification  a complex  model  is  typically 
subdivided  into  several  regions  and  the  grid  points  in 
each  region  are  labeled  with  a specific,  but  arbitrary,  series 
of  numbers.  Although  these  grid  point  numbers  act  only  as 
labels,  they  effect  the  bandwidth  of  the  stiffness  and 
mass  matrices.  In  order  to  minimize  the  matrix  bandwidth 
for  most  efficient  running  of  NASTRAN,  the  BANDIT  computer 
program  (Reference  5)  can  be  used  to  automatically  renumber 
and  assign  internal  sequence  numbers  to  the  grid  points. 

The  output  from  BANDIT  is  a set  of  SEQGP  cards  that  are 
then  included  in  the  NASTRAN  bulk  data  deck,  and  which  relate 
the  original  external  grid  numbers  to  the  internal  numbers, 
numbers . 


After  reviewing  many  of  the  structural  optimization  methods 
in  existence  (see  Appendix  A) , strain  energy  and  stress- 
ratio  resizing  techniques  were  employed.  For  applications 
such  as  helicopters  where  weight  is  critical,  it  is  more 
appropriate  to  evaluate  the  strain  density  (strain 
energy/volume)  distribution  within  a structure  which 
provides  guidance  for  structural  optimization.  A strain 
density  analysis  for  dynamics  applications  has  been 


5.  Everstine,  G.,  BANDIT  - A COMPUTER  PROGRAM  TO  RENUMBER 
NASTRAN  GRID  POINTS  FOR  REDUCED  BANDWIDTH,  Naval  Ship 
Research  and  Development  Center  Technical  Note  AML-6-70, 
February  1970. 


developed  by  Boeing  Vertol  under  ARO  Contract  DAHC04-71- 
C-0048  (Reference  6) . Expanding  upon  this  work,  a post- 
processor program  (S-83)  compatible  with  NASTRAN  has  been 
developed  for  analysis  of  the  strain  density  distribution 
throughout  a structure  for  dynamics  application  and  is  based 
on  the  concept  that  for  a given  load  condition,  a uniform 
strain  density  under  distortion  yields  a maximum  stiffness/ 
minimum  weight  structure.  This  program  uses  stress  data 
output  by  NASTRAN,  calculates  the  strain  density  of  each 
element,  and  tabulates  these  from  highest  to  lowest.  By 
employing  the  ALTER  feature  of  NASTRAN,  a checkpoint  tape 
containing  the  stresses  for  each  element  is  generated.  The 
post-processor  program  uses  the  data  stored  on  the  check- 
point tape  to  calculate  the  strain  density  of  NASTRAN  plate 
elements  and  tabulates  the  elements  in  descending  order  of 
strain  densities.  The  structural  elements  with  the  highest 
strain  density  are  the  best  candidates  for  effective 
optimization  since  a minimal  weight  change  will  yield  a 
maximum  benefit.  By  locally  altering  the  housing  wall  to 
change  the  mass  and  stiffness  in  these  areas  of  high  strain 
density,  the  structure  can  be  optimized.  This  strain 
density  distribution  concept  can  also  be  utilized  static- 
ally to  identify  structural  load  paths  and  to  evaluate  the 
efficiency  of  the  housing  structural  design  (stiffness/ 
weight) . The  S-83  strain  density  computer  program  can  be 
used  for  both  static  and  dynamic  analyses  (including 
thermal  effects) . Under  modal  distortion  for  a given 
natural  frequency,  a uniform  strain  density  yields  a 
maximum  lowest  eigenvalue/minimum  weight  structure.  The 
computer  model,  pre-  and  post-processor  programs,  and  all 
computer  data  and  format  are  compatible  with  NASTRAN 
(commercial  or  Schwendler  version) . The  capabilities  of 
NASTRAN  for  automated  plotting  and  analysis  of  composite 
materials  are  also  utilized. 

In  order  to  be  able  to  realistically  determine  thermal 
distortions/stresses  for  a new  or  conceptual  transmission 
case,  it  is  necessary  to  calculate  the  output  of  the 
heat  sources  which  are  the  gear  meshes  and  bearings. 

Since  this  capability  does  not  exist  in  NASTRAN,  the 
results  of  the  existing  Boeing  Vertol  computer  program,  which 
have  been  extensively  used  and  thoroughly  validated  for 
the  analysis  of  bearings  (S04)  and  gears  (R20  bevel,  R23 
spur  and  helical),  were  applied  to  calculate  power 
dissipation  and  the  resultant  heat  generation.  Using 
these  results  as  the  driving  potential  for  the  thermal 
analysis,  NASTRAN  was  then  run  for  a steady-state  heat 


6.  Sciarra,  J. , USE  OF  THE  FINITE  ELEMENT  DAMPED  FORCED 

RESPONSE  STRAIN  ENERGY  DISTRIBUTION  FOR  VIBRATION  REDUC- 
TION, U.S.  Army  Research  Office  - Durham,  Final  Report 
Contract  DAH-C04-71-C-0048 , July  1964. 
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solution  (heat  balance)  using  the  housing  model.  Cooling 
was  provided  by  modeling  the  oil  flow  and  incorporating 
it  into  the  solution.  Nodal  temperatures  were  punched  out 
automatically  on  cards  in  a format  compatible  with 
NASTRAN  input.  A thermal  distortion/stress  analysis 
followed  using  the  same  basic  model.  This  was  accomplished 
using  Rigid  Format  1 (Static  Analysis)  of  NASTRAN. 

The  use  of  computer  programs  that  are  proprietary  to 
Boeing  Vertol  presents  no  gap  in  the  final  product  to  be 
delivered  under  this  contract.  The  bevel  gear  program 
(R20)  is  based  on  Gleason  commercially  available  dimension 
sheet  methods,  and  the  spurAielical  gear  program  (R23)  is 
based  on  standard  AGMA  information.  The  rolling  element 
bearing  analysis  program  (S04)  was  purchased  by  Boeing 
Vertol  from  Mr.  A.  B.  Jones,  bearing  consultant,  and  is 
also  commercially  available  to  other  users.  Furthermore, 
most  industrial  and  commercial  potential  users  of  the 
finite  element  structural  analysis  method  have  their  own  gear/ 
bearing  analyses. 

The  present  NASTRAN  transmission  model  uses  elements  with 
a solid  homogeneous  cross-section  which  allow  for  both 
membrane  and  bending  stiffness.  The  material  properties 
are  isotropic.  These  conditions  represent  the  solid  cast 
structure  of  current  helicopter  transmissions.  However,  the 
following  options  exist  within  NASTRAN  to  allow  analysis  of 
structures  made  of  composite  materials : 

• A solid  homogeneous  cross-section  element  with  anisotropic 
material  properties. 

• Sandwich  plate  elements  that  can  reference  different 
materials  for  membrane,  bending,  and  transverse  shear 
properties,  each  of  which  may  have  either  isotropic  or 
anisotropic  material  properties. 

• A general  element  whose  properties  are  defined  directly 
by  the  user  in  terms  of  influence  coefficient  or 
stiffness  matrices. 

The  particular  element  chosen  is  dependent  upon  the  type 
cf  composite  material  to  be  represented.  The  NASTRAN 
model  may  also  be  modified  to  analyze  fabricated  truss 
type  structures. 

The  flowchart  presented  in  Figure  11  sximmarizes  the  overall 
transmission  housing  analysis.  The  previously  existing  work 
and  the  items  developed  herein  have  been  identified,  and  it 
is  shown  how  the  work  herein  contributes  to  the  overall  design 
goals . 
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Figure  1 1.  Flow  Diagram  of  Transmission  Housing  Analysis 


DES IGN  ANALYSIS  OF  TRANSMISSION  HOUSING 


I . DEFINITION  OF  MODEL 


The  contractor  selected  the  forward  main  rotor  transmission 
of  the  CH-47C  helicopter  for  consideration  under  this 
program.  A NASTRAN  finite  element  model  of  the  CH-47 
forward  rotor  transmission  constructed  by  Boeing  Vertol 
(under  USAAMRDL  Contract  DAAJ02-7 4-C-0040 ) for  vibration/ 
noise  reduction  studies  and  correlated  with  data  from  the 
HLH/ATC  noise  reduction  program  (References  7 and  8)  was 
used  and  improved  for  the  thermal  and  stress  analyses 
herein.  This  model  is  described  briefly  below  and  in  more 
detail  in  References  9 and  10. 

The  transmission  model  was  used  both  as  a baseline  concep- 
tual housing  for  generalized  studies,  along  with  simpler 
models  to  illustrate  cooling  fins  and  bearing  races/gear 
interfaces  for  heat  transfer  analyses,  as  well  as  being 
used  for  a specific  analysis  of  the  operating  conditions 
experienced  by  an  actual  CH-47C  transmission.  This 
generalized  conceptual  housing  model  was  used  to  address 
various  transmission  design  areas  and  to  develop  a general 
tool  for  specific  application  to  the  redesign  of  the 
CH-47C  forward  transmission.  The  conceptual  baseline  was 
used  to  develop  the  design  constraints  necessary  to  meet 
the  goals  of  this  contract  and  to  demonstrate  the  analytical 


7.  Hartman,  R.  M. , and  Badgley,  R. , MODEL  301  HLH/ATC  TRANS- 
MISSION NOISE  REDUCTION  PROGRAM,  Boeing  Vertol  Company, 
USAAMRDL  TR  74-58,  Eustis  Directorate,  U.S.  Army  Air 
Mobility  Research  and  Development  Laboratory,  Fort  Eustis, 
Virginia,  May  1974,  AD784132. 

8.  Hartman,  R.  M. , A DYNAMICS  APPROACH  TO  HELICOPTER  TRANS- 
MISSION NOISE  REDUCTION  AND  IMPROVED  RELIABILITY,  Paper 
Presented  at  the  29th  Annual  National  Forum  of  the  American 
Helicopter  Society,  Washington,  D.C.,  May  1973,  Preprint 
No.  772. 

9.  Sciarra,  J.  J.,  Howells,  R.  W. , and  Lemanski , A.  J.,  HELI- 
COPTER TRANSMISSION  VIBRATION  AND  NOISE  REDUCTION  PROGRAM, 
USAAMRDL  Contract  DAAJ02-74-C-0040 , Interim  Report, 

October  1975. 

10.  Howells,  R.  W. , and  Sciarra,  J.  J. , FINITE  ELEMENT  ANALYSIS 
USING  NASTRAN  APPLIED  TO  TRANSMISSION  VIBRATION/NOISE 
REDUCTION,  NASA  TMX-3278,  September  1975. 
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technique  developed  through  design  of  a housing  to  meet 
the  design  constraints  selected^  Several  areas  of  potential 
improvement  needing  further  consideration  have  been 
ident  d.f  ied ; 

• Housing  Material  - Desirable  Properties 

- Inherently  high  stiffness  and  strength. 

- Good  high  temperature  strength  to  provide  gear 
and  shaft  support  even  when  running  under 
emergency  lube  conditions  and  under  operating 
temperatures  of  350Of. 

- Use  material  forms  like  composites  or  rolled 
sheet  which  lend  themselves  to  very  thin  sections 
in  areas  where  oil  containment  is  the  main 
requirement . 

- The  housing  material  strength/elongation  proper- 
ties should  accept  the  hydraulic  ram  effect 
that  cracks  castings  after  ballistic  impact. 

Should  not  support  combustion  as  magnesium  does. 

- Inherent  cooling  capabilities. 

• Housing  Construction  and  Geometry 

- Minimize  interfaces,  doubled  thicknesses,  and 
bolted  joints . 

Simplify  the  rotor  load  path  and  minimize  bending 
as  far  as  possible. 

- Locate  gears  as  rigidly  as  possible  by  reacting 
their  loads  through  a direct,  short  and  rigid 
load  path.  Locate  bearings  to  minimize  bending 
loads  on  the  support  structure. 

Insure  rigid  support  locations  at  housing/ 
rotating  component  interfaces. 

The  above  design  features  and  the  increasing  perform- 
ance demands  placed  upon  helicopter  transmissions  have 
led  to  interest  in  the  possible  use  of  advanced  com- 
posite materials  for  transmission  housings.  Although 
composites  have  been  applied  with  much  success  for  many 
aircraft  parts,  their  use  for  transmission  housings  is 
new.  Many  design  areas  must  be  addressed:  load  path 
definition  for  orientation  of  non-isotropic  material 
properties;  fabrication,  tooling  and  machining  methods; 
thermal  properties,  etc. 
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A USAAMRDL/BHC  program  (Reference  11)  was  previously  con- 
ducted to  evaluate  carbon/epoxy  as  a possible  transmission 
housing  material.  Machining  of  the  carbon/epoxy  material 
was  found  to  be  difficult,  and  grinding  had  to  be  used 
almost  exclusively  to  finish  the  composite  housing.  Also, 
the  thermal  conductivity  and  structural  integrity  of  the 
carbon/epoxy  housing  was  found  to  be  poor.  Stiffness  was 
slightly  improved  as  evidenced  by  gear  development  test- 
ing. The  lack  of  success  of  this  particular  program  should 
not  deter  the  pursuit  of  composites  as  a housing  material. 
Conversely,  the  problems  encountered  during  the  program 
indicate  the  need  for  accurate  analysis  and  careful  design 
as  well  as  the  development  of  acceptable  fabrication 
procedures . 

For  example,  structural  failures  occurred  at  bonded  joints 
in  tension.  Such  a design  does  not  utilize  the  properties 
of  composites  properly  and  must  be  avoided.  Another  ex- 
ample is  the  concern  over  the  low  thermal  conductivity  of 
nonmetal  composites.  Graphite/polyimide  has  a lower 
thermal  conductivity  than  magnesium  by  a factor  of  about 
25  to  1,  thus  heat  dissipation  from  the  housing  will  be 
less  than  that  of  a magnesium  equivalent  and  areas  of 
higher  temperature  will  be  more  localized.  Since  convec- 
tion cooling  through  a magnesium  case  accounts  only  for  an 
estimated  15%  of  the  total  heat  rejection,  the  estimated 
weight  penalty  for  the  cooling  system  resultant  is  not 
prohibitive . 

The  second  probable  effect  of  the  thermal  conductivity 
difference  is  in  the  areas  immediately  adjacent  to  the 
bearings.  Bearing  heat  will  not  be  carried  away  as  fast  as 
it  is  with  a magnesium  housing.  The  coefficient  of  thermal 
expansion  of  molded  polyimide  with  chopped  graphite  rein- 
forcement is  6 X 10  ® inch/inch°F,  which  is  very  close  to 
steel  (6.5  x 10  ®).  The  effect  on  the  bearing  should  be 
beneficial  under  conditions  of  both  normal  operation  and 
lube  interruption.  That  is,  the  temperature  of  the  outer 
race  should  remain  closer  to  the  inner  race  under  all  con- 
ditions and  thus  eliminate  the  temperature  gradient  across 
the  bearing,  which  is  responsible  for  bearing  seizure  under 
abnormal  conditions.  Also,  the  housing  will  not  expand 
away  from  the  outer  race  when  the  thermal  growth  coeffi- 
cients are  matched.  This  will  also  tend  to  maintain  a 
constant  internal  clearance  and  stabilized  operating 
conditions  in  both  normal  and  abnormal  lubrication  situa- 
tions. The  conclusion  is  that  molded  chopped  graphite 


11.  Battles,  Roy  A.,  DYNAMIC  TESTING  OF  A COMPOSITE  MATERIAL 
HELICOPTER  TRANSMISSION  HOUSING,  Bell  Helicopter  Company, 
USAAMRDL  TR  75-47,  Eustis  Directorate,  U.S.  Army  Air 
Mobility  Research  and  Development  Laboratory,  Fort  Eustis, 
Virginia,  September  1975,  AD  A015521. 
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composite  bearing  supports,  with  thermal  coefficients 
matched  to  bearing  materials,  will  improve  the  ability  to 
operate  at  design  conditions  as  well  as  at  abnormal  lube 
conditions . 

The  finite  element  analysis  considered  herein  is  an 
essential  tool  necessary  to  handle  the  composite  material 
properties  most  efficiently.  Such  analyses  in  combination 
with  advanced  materials  promise  to  yield  helicopter  trans- 
missions with  substantially  improved  performance  and  lower 
life-cycle  cost. 

By  utilizing  the  existing  model  of  the  CH-47C  forward  rotor 
transmission,  the  contractor  has  been  able  to: 

• Cost-effectively  conduct  the  thermal/stress  analysis 
presented  herein. 

• Concentrate  effort  on  the  application  of  the  model  to 
derive  useful  information  rather  than  on  model 
building . 

• Directly  apply  results  of  other  related  contracts. 

Further  reasons  for  selecting  this  model  include; 

• Model  validation  has  provided  confidence  in  its 
accuracy. 

• Use  of  a widely  accepted  and  throughly  validated 
computer  program  (NAS TRAN) . 

• Extensive  computer-generated  plotting  capability 
used  to  debug  the  model. 

• Cross-checking  of  the  model,  design  drawings,  and 
hardware . 

• Good  correlation  of  the  model  and  hardware  weights 
(Table  2) . 

• Test  data  is  available  from  previous  Boeing  Vertol 
programs  for  model  validation  and  correlation. 

- Housing  displacements  from  HLH/ATC  noise  reduc- 
tion program  (Reference  7) . 

- Housing  temperatures  from  thermal  mapping  program 
(Reference  12) . 

• Hardware  is  available  for  further  testing  and 
modification. 

12.  Tocci,  R.  C. , Lemanski,  A.  J.,  and  Ayoub,  N.  J.,  TRANS- 
MISSION THERMAL  MAPPING,  Boeing  Vertol  Company,  USAAMRDL 
TR  73-24,  U.S.  Army  Air  Mobility  Reisearch  and  Development 
Laboratory,  Fort  Eustis,  Virginia,  May  1973,  AD  767875. 


33 


TABLE  2.  SUMMARY  OF  CH-47  FORWARD  TRANSMISSION  HOUSING 
NASTRAN  MODEL. 
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The  Boeing  Vertol  CH-47  forward  rotor  transmission  housing 
is  composed  of  three  major  sections;  upper  cover,  ring 
gear,  and  case  (including  sump) . The  upper  cover  provides 
lugs  for  mounting  the  transmission  to  the  airframe  and 
transmits  the  rotor  system  loads.  The  case  contains  and 
supports  the  main  bevel  gears.  The  ring  gear,  which 
connects  the  upper  cover  and  case,  contains  the  planetary 
gear  system.  This  natural  division  of  the  housing  was 
adhered  to  for  ease  of  modeling  (Figure  12) . 

The  geometric  grid  points  for  the  model  were  defined  from 
design  drawings  and  by  cross-checking  on  an  actual  housing. 
CQUAD2  (Quadrilateral)  and  CTRIA2  (Triangular)  homogeneous 
plate  (membrane  and  bending)  elements  were  used  to  connect 
the  grid  points  and  to  build  the  NASTRAN  structural  model.  A 
Boeing  Vertol  preprocessor  program  (SAIL  II  - Structural 
Analyses  Input  Language)  for  the  automatic  generation  of 
grid  point  coordinates  and  structural  element  connections 
was  used.  This  preprocessor  allows  the  user  to  take 
advantage  of  any  pattern  which  occurs  in  the  data  by 
providing  straightforward  techniques  for  describing  algo- 
rithms to  generate  blocks  of  data.  The  extensive  computer 
generated  plotting  capability  of  NASTRAN  was  used  to  de- 
bug the  structural  model. 

For  ease  of  identification  the  housing  was  subdivided  into 
several  regions,  and  the  grid  points  in  each  region  were 
labeled  with  a specific,  but  arbitrary,  series  of  numbers. 
Although  these  grid  point  numbers  act  only  as  labels,  they 
effect  the  bandwidth  of  the  stiffness  and  mass  matrices. 

In  order  to  minimize  the  matrix  bandwidth  for  most  efficient 
running  of  NASTRAN,  the  BANDIT  computer  program  (Reference 
5)  was  used  to  automatically  renumber  and  assign  internal 
sequence  numbers  to  the  grid  points.  The  output  from 
BANDIT  is  a set  of  SEQGP  cards  that  are  included  in  the 
NASTRAN  bulk  data  deck  and  which  relate  the  original  external 
grid  numbers  to  the  internal  numbers . 

The  model  includes  grid  points  representative  of  the 
structure  where  the  shafts  are  supported  by  their  bearings. 
These  grid  points  are  used  to  apply  the  dynamic  and 
static  loads  to  the  housing.  Each  geometric  grid  point 
has  six  possible  degrees  of  freedom  (three  translational 
and  three  rotational) . To  conveniently  evaluate  the  motion 
normal  to  the  housing  surface,  numerous  local  coordinate 
systems  were  defined  and  oriented  such  that  the  displace- 
ments and  accelerations  calculated  at  each  grid  point 
could  be  referred  to  as  a coordinate  system  having  an  axis 
normal  to  the  housing  surface.  One  degree  of  freedom, 
rotation  about  the  normal  to  the  surface,  was  constrained 
since  the  stiffness  for  this  component  is  undefined  for 
NASTRAN  plate  elements.  The  other  two  rotational  degrees 
of  freedom  were  omitted.  All  translational  degrees  of 
freedom  were  retained  to  accurately  represent  the  motion 
of  the  actual  housing.  The  model  parameters  are  summarized 
in  Table  2. 
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CUT-AWAY  VIEW 


Figure  13.  CH-47C  Forward  Transmission  Internal  Components 


II.  EXTENSION  OF  MODEL  FOR  INTERNAL  COMPONENTS 


The  original  transmission  model  included  the  complete 
housing  structure  (Figure  12)  as  well  as  the  internal 
dynamic  components  (Figure  13) . The  model  of  these  internal 
components  was  used  to  calculate  the  dynamic  forces  that 
were  subsequently  applied  analytically  to  the  housing  model 
at  the  bearing  locations.  The  structural  aspects  of  the 
internal  components  were  not  considered.  Since  for  static 
load  conditions  the  internal  components  provide  additional 
structural  constraints  on  the  housing,  the  effect  of  these 
internal  components  must  be  accounted  for  in  the  thermal 
and  structural  models.  The  effects  of  bevel,  sun  lower  and 
upper  planetary  gears,  and  supporting  structure  were  con- 
sidered. It  is  generally  not  necessary  to  model  these  com- 
ponents in  full  detail,  since  only  their  gross  effect  on  the 
housing  is  desired.  A simple  beam  representation  is  adequate. 


Including  the  internal  components  in  the  thermal  model 
was  neither  necessary  or  desirable.  The  bearing  outer 
races  (steel  o<.s  - 7.5  X 10“^  in/in/°F) , which  are  the 
case/internal  component  interface,  are  press  fit  into  the 
housing.  Elevated  temperatures  cause  the  magnesium  case 
(o<r^  = 15.0  X 10“^  in/in/°F)  to  expand  away  from  the  outer 
race  and  may  even  result  in  a "floating"  fit  at  operating 
temperature.  This  situation  was  experienced  in  Reference 
12,  where  it  was  necessary  to  key  the  outer  races  of  the 
spiral  bevel  pinion  bearings  in  their  case  liners  to  pre- 
vent rotation  that  would  be  permitted  by  increased  clearances 
caused  by  thermal  expansion/growth.  This  condition,  plus 
the  internal  tolerance  within  a bearing,  precludes  the  trans- 
mittal of  radially  outward  thermal  loads  into  the  housing. 
Furthermore,  it  is  not  possible  for  the  bearing  races  to 
impose  radial  restraint  upon  the  housing  expansion.  Thus, 
no  representation  of  the  internal  components  in  the  radial 
direction  is  necessary.  In  the  axial  direction,  thermal 
growth  of  the  shafting  is  absorbed  in  the  form  of  reduced 
gear  backlash.  Thus,  no  axial  loads  are  generated  unless 
the  temperature  exceeds  that  necessary  to  reduce  the  back- 
lash to  zero.  In  such  a situation,  the  housing  loads 
would  be  of  little  interest  since  the  gears  would  distress 
and  fail.  Figures  14,  15  and  16  indicate  the  changes  in 
bevel  gear  backlash  and  root  clearance  due  to  elevated 
temperatures.  Allowing  for  these  effects,  a high  tempera- 
ture transmission  must  be  designed  to  maintain  adequate 
backlash  as  well  as  bearing  clearances. 

The  only  significant  property  of  the  internal  components 
in  the  static  stress  model  is  the  radial  resistance  'of 
the  bearing  outer  races  to  compressive  (radially  inward) 
forces  acting  on  the  housing. 
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CHANGE  IN  BACKLASH (aB) 


AXIAL  MOVEMENT  ALONG  OF  MEMBER  (Z) 


Figure  14.  Change  in  Backlash  Due  to  Axial  Movement 
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Figure  15.  Change  in  Backlash  Due  to  Temperature  Change 


Figure  16.  Change  in  Root  Clearance  Due  to  Temperature  Change 
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bearing  races  offer  no  resistance  to  radially  outward 
forces  on  the  housing.  A problem  arises  here  due  to  the 
nature  of  the  beam  element.  No  capability  exists  in 
NASTRAN  for  a beam  which  will  act  in  compression  only. 

Thus,  a beam  model  of  a bearing  race  will  also  act  to 
impose  unwanted  restraint  on  the  housing  directed  radially 
inward.  This  could  be  circumvented  by  first  analyzing  the 
housing  model  without  including  the  beams  representing  the 
races  and  thereby  defining  those  housingA>earing  interfaces 
with  radially  inward  deflections.  A beam  model  could  then 
be  inserted  at  these  points  to  resist  the  radially  inward 
forces  and  the  analysis  could  be  rerun.  Figure  17  indi- 
cates the  be a ring /ho using  interfaces. 


Figure  17.  Shaded  Lines  indicate  Bearing  Outer  Race/Housing  Interfaces 
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III.  THERMAL  DISTORTION  AND  THERMAL  STRESS  ANALYSIS 


1.  General 

A complete  thermal  analysis  must  consider  thermally  in- 
duced deflections  and  stress  as  well  as  heat  transfer 
characteristics.  The  capability  for  predicting  the 
thermal  distortion  of  a transmission  housing  under  op- 
erating conditions  is  important  from  the  aspects  of 
load  capacity  improvement  and  avoidance  of  premature 
component  distress.  The  sensitivities  of  bevel  gear 
load  capacity  to  tooth  pattern  and  of  bearing  perfor- 
mance to  mounting  tolerances,  both  of  which  are  affect- 
ed by  the  dimensional  stability  of  the  housing,  are 
well  known.  The  thermally  induced  stress  distribution 
in  the  housing  must  be  known  and  then  superimposed  upon 
the  stress  due  to  static/dynamic  loads.  Only  with  such 
a detailed  knowledge  of  the  housing  stress  can  a 
structural  design  be  optimized  for  maximum  efficiency. 
Furthermore,  in  line  with  the  long  range  objective  of 
an  integrally  lubricated  transmission,  a thermal  analy- 
sis is  required  for  predicting  the  heat  flow  character- 
istics of  the  transmission  housing.  Since  such  a 
transmission  will  operate  at  higher  temperatures  than 
current  conventionally  cooled  systems,  the  effects  of 
thermal  distortion  and  stress  will  become  more  critical. 
Such  studies  also  bear  directly  on  the  development  of 
transmission  capability  to  endure  relatively  long  peri- 
ods of  emergency  hot  running  after  loss  of  lubricant 
(Figure  18) . The  sequence  of  development  required  to 
arrive  at  an  integrally  lubricated  transmission  is 
shown  in  Figure  19.  Test  results  have  also  indicated 
that  thermal  growth  may  be  a contributor  to  premature 
gear  distress  at  present  operating  temperatures. 

Thermal  growth,  being  proportional  to  both  temperature 
change  and  component  size,  becomes  more  critical  with 
larger  transmissions  (e.g.,  HLH)  and  higher  temperatures 
(e.g.,  sealed-lube  systems).  Thus,  advanced  trans- 
mission designs  will  need  improved  thermal  analyses  to 
alleviate  such  problems. 

The  investigations  documented  herein  indicate  the  rela- 
tive contribution  to  distress  of  transmission  compon- 
ents when  subjected  to  differential  thermal  growth. 

The  results  of  this  analysis  indicate  the  work  required 
to  arrive  at  acceptable  drive  system  mission  reli- 
abilities when  differential  thermal  growth  takes  place 
in  hot  running  helicopter  transmissions  for  continuous 
and  emergency  periods  of  time.  The  finite  element 
transmission  housing  model  discussed  in  a previous 
section  is  ideally  suited  for  application  to  a thermal 
analysis.  Figure  20  describes  the  thermal  model  of 
the  housing. 
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CALENDAR  YEARS 


Figure  18.  Emergency  OH-Off  Capability  Goals 


OIL  LINES 


EXISTING  CH-47C  LUBE  SYSTEM 
(VULNERABLE  AREA - 
11.4SQFT) 


EVOLUTION  OF  A 
SURVIVABLE 
TRANSMISSION  SYSTEM 


OIL  LINES  ^ r i 
BY-PASS  VALVE 


95% IMPROVEMENT 


NON-LUBE  CAPABILITY 
(VULNERABLE  AREA  - 0.6 SO  FT) 


OIL  PUMP 


ARMOR 


80%  IMPROVEMENT 


BY-PASS  LUBE  SYSTEM 
(VULNERABLE  AREA- 
4.1  SO  FT) 


ARMORED  SUMP 
(VULNERABLE  AREA - 
2.3  SO  FT) 


Figure  19.  Sequence  of  Development  Required  to  Arrive  at  a Factory 
Sealed  Advanced  Technology  Transmission 
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UPPER  COVER 

GRID  POI NTS  = 160 
ELEMENTS  = 202 


■ 

t 


AUTOMATIC  GENERATION  OF  GRID 
POINTS  AND  ELEMENT  CONNECTIONS. 
USE  SAME  GRIDS  FOR  STRUCTURAL 
AND  THERMAL  ANALYSIS. 

HOUSING  MADE  OF  PLATE  ELEMENTS 
METAL  OR  COMPOSITE  MATERIAL. 
THERMAL  DISTORTION  AND  THERMAL 
STRESS. 

GEAR/BEARING  ANALYSIS  TO  DETER- 
MINE HEAT  GENERATION,  USE  AS 
DRIVING  POTENTIAL  FOR  MODEL. 

HEAT  TRANSFER  BOUNDARY  CONDI- 
TIONS AT  EXTERNAL  SURFACES 
DETERMINED  FROM  ENVIRONMENT. 
THERMAL  (HEAT  TRANSFER)  ANALYSIS 
DEFINES  TEMPERATURE  FIELD. 

DEFINE  RESPONSE  (DEFORMATION)  AT 
GEAR  MESHES  AND  BEARINGS  DUE  TO 
HEAT  GENERATED. 

STRUCTURAL  OPTIMIZATION. 
VALIDATION  OF  ANALYTICAL  TOOLS 
FOR  APPLICATION  TO  ADVANCED 
TRANSMISSION  DESIGN. 


Figure  20.  Model  for  Thermal  Analysis 
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Steady-state  heat  transfer  analyses  may  be  accomplished 
with  the  existinq  finite  element  model  of  the  CH-47C 
forward  transmission  housing  using  the  NASTRAN  program. 
The  NASTRAN  heat  flow  capability  analysis  is  compatible 
with  the  structural  analysis  capability.  Using  the 
same  model  (grid  points,  coordinate  system,  elements, 
constraints  and  sequencing)  as  used  for  a structural 
analysis,  the  steady-state  temperatures  at  each  grid 
point  may  be  calculated  and  provided  directly  in 
punched  form  for  later  use  in  a thermal  stress/defor- 
mation analysis  (static  structural  analysis) . This 
NASTRAN  housing  model  may  also  be  modified  to  handle 
a thermal  analysis  of  fabricated  housings  and/or 
composite  materials  such  as  metal-matrix,  which  will 
be  particularly  valuable  for  application  to  advanced 
transmission  development. 

Input  to  the  model  for  a thermal  analysis  requires 
knowledge  of  the  thermal  boundary  conditions  at  spec- 
ified locations  on  the  housing.  The  driving  potential 
for  a thermal  analysis  is  the  power  dissipation  by 
bearings  and  gears  at  operating  conditions.  This 
originated  from  other  Boeing  Vertol  computer  programs. 
The  temperature  distribution  obtained  from  the  steady- 
state  temperature  analysis  can  then  be  input  to  NASTRAN 
to  provide  a thermal  stress  (static)  analysis.  Thermal 
stresses  and  distortions  result.  The  thermal  distor- 
tions obtained  are  of  particular  use  in  determining  the 
adequacy  of  clearances  for  a longer  life  transmission. 
The  thermal  analysis  is  outlined  in  Figure  21. 

A sample  NASTRAN  thermal  stress  analysis  presented 
below  shows  good  correlation  and  indicates  the  feasi- 
bility of  using  NASTRAN  for  thermal  analyses.  The 
model,  a rectangular  plate,  is  given  a temperature 
gradient  which  causes  internal  loads  and  elastic 
deflections.  The  temperature  load  is  constant  in  the 
y direction  and  symmetric  about  the  y-axis.  Since 
membrane  elements  are  used  to  model  the  structure, 
it  is  necessary  to  remove  all  rotational  degrees  of 
freedom  and  translational  degrees  of  freedom  normal  to 
the  plate.  The  symmetric  boundary  conditions  were 
modeled  by  constraining  the  displacements  normal  to 
the  planes  of  symmetry.  Figure  22  shows  a comparison 
of  stresses  predicted  by  NASTRAN  and  the  experimentally 
measured  stresses. 
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Figure  21.  Flow  Diagram  of  NASTRAN  Thermal  Analysis 
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Figure  22.  Comparison  of  NASTRAN  and  Experimental  Stresses  for  Free 
Rectangular  Plate  With  Thermal  Loading 


13.  NASTRAN  User's  Manual  Level  15,  June  1972 
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2 . Uniform  Temperature  Testing 


The  initial  phase  of  the  thermal  work  was  the  appli- 
cation of  a uniform  temperature  distribution  over  the 
housing  model.  This  was  done  analytically  using 
NASTRAN  (Rigid  Format  1)  for  temperatures  representa- 
tive of  current  operating  (160°F) , projected  operating 
(350°F) , and  projected  loss-of-lubricant  emergency 
operating  (700°F)  conditions,  and  required  a total  of 
40  minutes  execution  time  on  an  IBM  370  computer. 

Although  the  thermal  stresses  are  invariant,  boundary 
fixity  was  found  to  influence  the  point -by -point 
distortions  and  to  make  interpretation  difficult. 

Studies  were  made  to  resolve  this,  and  the  housing  was 
analyzed  for  thermal  stresses  and  distortions.  The 
thermally  induced  deformations  were  defined  (Figure  23) . 

For  model  validation,  the  baseline  housing  was  heated 
to  temperatures  in  the  range  of  160°-400°F.  By  measuring 
selected  dimensions  of  the  housing  at  normal  and 
elevated  temperatures,  this  testing  experimentally 
determined  the  thermal  distortion  of  the  transmission 
case.  Correlation  with  the  dimensional  changes  pre- 
dicted analytically  by  NASTRAN  provided  confidence  in 
the  thermal  model  for. further  analysis. 

Thermocouples  were  installed  on  the  CH-47  forward 
rotor  transmission  housing  at  five  locations  dispersed 
over  the  housing  and  representative  of  various  wall 
thicknesses  occurring  in  the  l^ousing.  The  points  for 
dimensional  inspection  were  defined  to  provide  repre- 
sentative distortions  and  to  be  readily  accessible  for 
measurement  under  elevated  temperature  conditions. 

The  thermocouple  locations  and  points  for  dimensional 
inspection  are  shown  in  Figure  24. 

After  measuring  and  recording  the  room  temperature 
dimensions,  the  housing  was  placed  in  a 204  kW  walk-in 
type  electric  oven  (Figure  25)  and  heated  to  each  of 
the  following  three  approximate  target  temperatures: 
160®F,  280“F,  and  400“F.  The  temperature  at  each 
thermocouple  was  recorded  continuously  on  paper  tape. 
After  the  temperature  in  the  oven  stabilized,  the 
housing  was  allowed  to  "soak"  at  the  elevated  tempera- 
ture for  approximately  one  hour  to  insure  a thorough 
and  uniform  heating.  During  soaking  the  temperatures 
indicated  by  the  five  thermocouples  remainej  constant 
and  within  a +5°F  band,  indicating  a uniform  heating 
of  the  housin. 

After  soaking,  t>  specified  housing  dimensions  were 
measured.  The  how  ig  remained  in  the  oven  during 
measuremert  in  order  to  minimize  heat  loss,  but  the 
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Figure  23.  CH-47  Forward  Transmission  Case  Uniform  Temperature  Analaysis,  Static 
Deformation  Due  to  Elevated  Temperatures 


INSPECTION  \ ^ 
POINT  NO.  3 

(GP2000+2180:  R ONLY) 


Figure  25.  CH-47C  Forward  Rotor  Transmission  Housing  in  Dispatch  Oven 

oven  doors  were  partially  open  and  the  oven  was  turned 
off.  During  the  measurement,  the  housing  temperature 
was  continuously  recorded  on  paper  tape  and  the  tape 
was  annotated  to  show  the  temperature  at  the  time  each 
dimensional  check  was  made  (Figure  26) . Typically, 
about  four  dimensions  could  be  measured  at  the  160 °F 
range  and  about  two  dimensions  at  the  400 °F  range  before 
the  housing  temperature  decreased  substantially.  The 
experimental  setup  is  shown  in  Figures  27  through  30. 

The  instruments  used  for  the  measurements  are  shown  in 
Figure  31.  Table  3 summarizes  the  equipment  used  in 
the  testing  and  indicates  the  reliability  cf  the 
dimensions  taken.  When  evaluating  the  reliability  of 
the  dimensional  data  taken  at  elevated  temperatures, 
consideration  must  be  given  to  the  facts  that  some 
growth  of  the  inspection  gages  occurred  when  taken  into 
the  oven  and  also  the  inspector  had  to  work  in  an 
uncomfortably  hot  environment. 

The  experimental  data  obtained  is  plotted  in  Figure  32 
as  the  change  in  linear  dimensions  versus  temperature. 
Also  shown  in  the  figure  are  the  theoretical  changes 
in  the  dimensions  as  predicted  both  by  the  NASTRAN 
thermal  analysis  and  by  a simple  linear  thermal  expan-  ' 
sion  calculation.  The  agreement  of  the  data  and 
analytical  methods  confirms  the  validity  of  the  model, 
which  can  thus  be  used  with  confidence  to  predict 
deformations  of  the  housing. 
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Figure  26.  Control  and  Data  Monitoring  Panel  for  Dispatch  Oven 


Figure  27.  Measurement  Procedure  — Bar-Type  Dial  Indicator  Gage 
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Figure  28.  Measurement  Procedure  — Outside  Micrometer 


Figure  29.  Measurement  Procedure  — Depth  Gage 
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Figure  30.  Measurement  Procedure  — Vernier  Caliper 


BAR-TYPE  DIAL 
INDICATOR  GAGE 

- . \ 


Figure  31.  Measuring  Instruments 
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TABLE  3.  EQUIPMENT  USED  TO  PERFORM  UNIFORM  TEMPERATURE  TESTING 
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Figure  32.  Transmission  Housing  Dimensiorial  Changes  as  a Function  of  Temperature 


3 . Gear/Bearinq  Thermal  Simulation 


The  temperature  distribution  throughout  a new  transmission 
is  obtained  during  the  design  phase  from  heat  transfer 
analyses  for  both  the  housing  and  internal  components. 

The  heat  output  of  the  gears  and  bearings  are  the  heat 
sources  within  the  transmission  and  are  the  forcing 
functions  for  the  thermal  analysis.  Since  the  develop- 
ment of  a completely  rigorous  analysis  of  gear/bearing 
heat  generation  was  beyond  the  scope  of  this  contract, 
within  the  transmission  and  are  the  forcing  functions 
for  the  thermal  analysis.  Since  the  development  of  a 
completely  rigorous  analysis  of  gear/bearing  heat 
generation  was  beyond  the  scope  of  this  contract, 
approximate  methods  were  applied.  Nevertheless,  the 
methods  used  for  the  component  heat  generation  yielded 
results  which  were  within  6%  of  independent  calculations 
based  on  oil-in/oil-out  cooling. 

The  dynamic  and  kinematic  parameters  output  by  existing 
Boeing  Vertol  computer  programs  (S04  rolling  element 
bearings,  R20  bevel  gears,  and  R23  spur/helical  gears), 
in  conjunction  with  information  derived  from  a review 
of  References  14,  15,  and  16,  were  used  to  predict 
temperatures.  Figure  33  defines  the  data  flow  for  the 
heat  generation  analysis.  The  results  of  this  work 
were  used  for  correlation  with  measured  thermal  data 
and,  in  the  next  section,  as  input  to  the  finite 
element  model  for  a thermal  analysis  of  the  housing. 

A thermal  distortion/stress  analysis  was  accomplished 
using  NASTRAN  Rigid  Format  1 (Static  Analysis) . 

Figure  21  indicates  the  overall  scheme  of  the  thermal 
analysis . 

The  computer  programs  used  to  predict  gear  temperatures 
are  based  on  an  analysis  of  the  temperature  at  a gear 
mesh,  or  the  flash  temperature,  which  is  used  to  deter- 
mine the  average  temperature  generated  over  several 
mesh  cycles.  It  had  previously  been  assumed  that  the 
power  loss  per  mesh  was  between  .5%  and  .75%.  The  .5% 
estimate  gave  an  error  of  9%  based  on  a comparison  of 
oil-in/oil-out  heat  dissipation.  Using  the  computer 
program  results,  however,  the  power  loss  calculated  for 


14.  Buckingham,  E. , ANALYTICAL  MECHANICS  OF  GEARS. 

15.  Harris,  T.,  ROLLING  BEARING  ANALYSES. 

16.  Rumbarger,  J.,  and  Filetti,  E.,  HIGH  TEMPERATURE 
HELICOPTER  MECHANICAL  TRANSMISSION  TECHNOLOGY,  Franklin 
Institute  Research  Laboratories,  Phila. , Pa.,  FIRL 
Final  Report  F-C3229,  December  1971. 
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Figure  33.  Flow  of  Heat  Generation  Analysis 


the  bevel  mesh  was  .2%  for  a better  system  correlation 
of  6%.  As  an  alternative,  a gear  temperature  formula- 
tion from  Reference  17  was  also  attempted: 


cos  0^  cos  ^2 

Efficiency  = 

cos  0^  cos  Vi  cos  2 + f 

where  = Pinion  = 7.55/2  Large  End  Pitch  Diameter 
R2  = Gear  = 13.278/2  Large  End  Pitch  Diameter 
n = RPM  Pinion  = 7460 
V = Pitch  Line  Velocity  = fpm 

Vg  = Sliding  Velocity  Between  Basic  Racks,  fpm 
0 = Pressure  angle  of  basic  rack  = 22.5° 


= Shaft  angle  = 99° 

= Assumed  coefficient  of  friction  f = 0.01 


V = .5236  n = (.5236)  (3.775)  (7460)  = 14745.36 

For  any  shaft  angle  Vg  = V sin^/cos  = 16069.  fpm 

= 14,745.36  sin  99°/cos  25° 

EFF  = (.92388)  (.906308)^ 

Q:os  (22.5)J  (.906308)^  + f sin  99° 


EFF  = .75887  

.75887  + (.01)  (.987688) 


= 98.7% 


HEAT  GENERATED  = (.013)  (3600)  (42.42)  = 1758  BTU/MIN 


17.  Buckingham  and  Ryffel,  DESIGN  OF  WORM  AND  SPIRAL  GEARS, 
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18. 


A computer  program  in  the  "BASIC" language  with  sample 
results  (BTU/MIN)  for  various  assumed  coefficients  of 
friction  is  also  included (Figures  34  and  35).  Reference 
18  was  also  utilized. 

10  INPUT  PI,  P2,  P3,  P4 
20  A=COS (PI) *COS (P2) *COS (P3) 

30  FOR  1=1  TO  50 
40  E=0.01*I 

50  E=E*SIN(P4) 

60  E=E+A 
70  E=A/E 

80  PRINT  P1;P2;P3;P4;E 
90  E=l-E 
100  B=E*3600*42.42 

110  F=0.01*I 

120  PRINT  "COEFF  OF  FRICTION="F ; "BTU/MIN="B 
130  NEXT  I 
140  END 

AVERAGE  POWER  LOSS  = 0.769HP  = 32 . 590BTU/MIN 

Figure  34.  "Basic”  Computer  Program  for  Heat  Generated  by  Spiral  Bevel 
Gears  (Hewlett  Packard  Minicomputer) 


722.5 

,25 

. ,25. ,99. 

22.5 

25  25 

99 

0.988488388 

COEFF 

OF 

FRICTION  = 

= 0.01 

BTU/MIN=  1757.961316 

22.5 

25  25 

99 

0.977238794 

COEFF 

OF 

FRICTION= 

0.02 

BTU/MIN=  3475.909312 

22.5 

25  25 

99 

0.966242373 

COEFF 

OF 

FRICTION= 

0.03 

BTU/MIN=  5155.194742 

22.5 

25  25 

99 

0.955490674 

COEFF 

OF 

FRICTION= 

0.04 

BTU/MIN=  6797.10824 

22.5 

25  25 

99 

0.944975617 

COEFF 

OF 

FRICTION= 

0.05 

BTU/MIN=  8402.88363 

22.5 

25  25 

99 

0.934689474 

COEFF 

OF 

FRICTION= 

0.06 

BTU/MIN=  9973.701007 

22.5 

25  25 

99 

0.924624852 

COEFF 

OF 

FRICTION= 

0.07 

BTU/MIN=  11510.689.64 

22.5 

25  25 

99 

0.914774669 

COEFF 

OF 

FRICTION= 

0.08 

BTU/MIN=  13014.93069 

22.5 

25  25 

99 

0.905132146 

COEFF 

OF 

FRICTION= 

0.09 

BTU/MIN=  14487.4597 

22.5 

25  25 

99 

0.895690784 

COEFF 

OF 

FRICTION= 

0.1 

BTU/MIN=  15929.26905 

Figure 

35. 

Sample  Computer  Output  for  Various  Assumed  Coefficients  i 

Faires,  V.,  DESIGN  OF  MACHINE  ELEMENTS,  The  MacMillan 
Company,  New  York,  4th  Edition,  1965. 
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The  output  from  computer  program  S04  has  been  applied 
to  predict  the  temperatures  of  rolling  element  bearings. 
Some  work  has  been  done  in  this  area  by  Rumbarger  and 
Filetti  (Reference  16) . Since  the  heat  generated  by 
the  bearings  is  comprised  of  two  parts,  the  shearing 
of  the  interface  oil  film  and  by  the  friction  of  the 
rolling  elements,  two  separate  calculations  must  be 
combined  (Figure  36) . Also,  the  heat  generation  is  a 
function  of  load,  sliding  velocity,  and  coefficient  of 
friction.  Sample  procedures  for  calculating  the  heat 
generated  by  both  angular  contact  ball  bearings  and 
roller  bearings  are  outlined  in  Appendix  B. 


Lubrication  characteristics  of  the  CH-47C  forward 
transmission  have  been  determined.  Table  4 indicates 
the  oil  flow  for  forced  convective  cooling  of  the 
internal  components,  and  Figure  37  indicates  the  orifice 
locations  numbered  according  to  Table  4.  The  stabi- 
lized thermal  power  at  the  bearings  and  gears  act  as 
heat  sources  for  a NASTRAN  analysis  including  conduc- 
tivity, natural  and  forced  convection,  and  radiation 
(Stefan-Boltzmann  Law)  (Figures  38  and  39) . Resulting 
overall  temperatures  can  then  be  used  in  a NASTRAN 
static  analysis  in  order  to  determine  thermal  distor- 
tions and  stresses. 

Thermal  power  has  been  determined  for  all  gear  meshes 
(Figure  40).  Samples  of  the  computer  output  for  the 
CH-47C  forward  transmission  are  shown  in  Figure  41 
(bevel  gear  program  R20)  and  Figure  42  (spur  and 
helical  gear  program  R23) . 

A tabulation  of  the  heat  generated  by  the  CH-47C 
forward  transmission  internal  components  is  provided 
in  Table  5. 
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r notion  Heat  Generation  Versus 
Load;  218  Angular-Contact  Ball 
Bearing  10,000  RPM,  5 Cent!- 
stokes  OH,  Jet  Lubrication 


RAYMOND  J.  DRACO 

BEVEL  GEAR  SURFACE  LOAD  CAPACrTY  ANALYSIS 
CONTACT  STRESS  - G-PACTOR  SCORING  HAZARD 


CH-47C  FWD  BEVEL  MESH  RUN  FOR  J.  S.  4/14/76 


NORMAL  PRESSURE  ANGLE 
SHAFT  ANGLE 
MEAN  SPIRAL  ANGLE 
NUMBER  OF  TEETH 

TRANSVERSE  DIAMETRAL  PITCH  OUTER 
MAXIMUM  SURFACE  FINISH  (MICROINCHES) 
EFFECTIVE  FACE  WIDTH 
PITCH  DIAMETER 
PITCH  ANGLE 

PITCH  LINE  VELOCITY  FPM 
REDUCTION  RATIO 

CONE  DISTANCE  OUTER  = 7.19384 


PINION  GEAR 

22.50000 

99.00000 

25.00000 

29.00000  51.00000 

3.84100 

20.00000 
2.18800 

7.55012  13.27779 

31.65230  67.34764 

14745.01000 
1.75862 

MEAN  = 6.09984 


STRESS: 

POWER 

SPEED 

TORQUE 

TANGENTIAL  TOOTH  LOAD 
CONTACT  RATIO  - FACE 

- PROFILE 

- MODIFIED 

INERTIA  FACTOR 


3600.000 

7460.00000  4241.96000 

30414.20000  53487.03000 

8056.61700 
1.50491 
1.25868 
1.96189 
1.01942 


KO  - 1,000  KV  - 1,000  KM  = 1,100  KS  = 0,714  CS  = 1,000  CF  = 
ELASTIC  COEFFICIENT 

PROFILE  CURVATURE  RADIUS  - PITCH  POINT 

- CRITICAL  POINT 

LENGTH  OF  CRITICAL  LINE  OF  CONTACT 
LOAD  SHARING  RATIO 

DISTANCE  BETWEEN  TOOTH  MIDPOINT  AND  CRITICAL  POINT 
CONTACT  GEOMETRY  FACTOR 
MAXIMUM  CONTACT  STRESS 


1.000 

2800.00000 

1.69785  6.59948 

1.87387  6.42346 

1.544348 
1.000000 
-0.017148 
0.094451 
211021.3000 


SCORING  HAZARD: 

THERMAL  CONSTANT,  Cl 
THERMAL  FACTOR,  CG 

PRCS'ILE  CURVATURE  RADIUS  AT  CRITICAL  POINT 
LENGTH  OF  CRITICAL  LINE  OF  CONTACT 
LOAD  SHARING  RATIO 

DISTANCE  FROM  TOOTH  MIDPOINT  TO  CRITICAL  POINT 
GEOMETRY  FACTOR,  G 
TEMPERATURE  RISE  AT  CRITICAL  POINT 
AVERAGE  POWER  LOSS  ■ 7,3178  H.P.  - 310,6389  BTU/MIN 


\ 


HEAT  GENERATED 


41.00000  41.00000 

1.00000 

1.97082  6.32651 

1.46105 
0.70570 
0.18285 
0.00178 
215.50340 


Figure  41.  Gear  Computer  Program  (R20)  for  Calculating  Heat  Generated  by 
Spiral  Bevel  Mesh.  Sample  Output  CH-47C  Forward  Transmission 
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TAN. 

SCORNG 

FRICT. 

TEMP. 

FINAL 

FILM 

INSTANTANEOUS 

LOAD 

FACTOR 

COEFF . 

RISE 

TEMP. 

THCK. 

POWER  LOSS 

LBS. 

F DEG 

F DEG 

MICRO 

HP 

BTU/MIN 

0. 

0.0026 

0.0600 

0.00 

200.00 

0.000 

0.00 

0.00 

4773. 

0.0009 

0.0600 

18.64 

218.64 

12.245 

1.23 

52.24 

4773. 

0.0000 

0.0600 

0.00 

200.00 

13.624 

0.00 

0.00 

4773. 

0.0008 

0.0600 

16.91 

216.91 

15.053 

1.29 

54.74 

682. 

0.0020 

0.0600 

9.15 

209.15 

22.210 

0.47 

20.05 

0. 

0.0026 

0.0600 

0.00 

200.00 

0.000 

0.00 

0.00 

682. 

0.0023 

0.0600 

10.61 

210.61 

13.603 

0.39 

16.52 

1364. 

0.0020 

0.0600 

15.59 

215.59 

12.821 

0.69 

29.46 

2046. 

0.0017 

0.0600 

18.17 

218.17 

12.536 

0.92 

38.80 

2728. 

0.0014 

0.0600 

19.00 

219.00 

12.438 

1.05 

44.50 

3410. 

0.0012 

0.0600 

18.41 

218.41 

12.436 

1.10 

46.50 

4773. 

0.0009 

0.0600 

18.63 

218.63 

12.246 

1.23 

52.22 

4773. 

0.0007 

0.0600 

13.74 

213.74 

12.589 

0.93 

39.25 

4773. 

0.0005 

0.0600 

8.99 

208.99 

12.935 

0.62 

26.17 

4773. 

0.0002 

0.0600 

4.38 

204.38 

13.282 

0.31 

12.99 

4773. 

0.0000 

0.0600 

0.10 

200.10 

13.632 

0.01 

0.31 

4773. 

0.0002 

0.0600 

4.47 

204.47 

13.984 

0.32 

13.73 

4773. 

0.0004 

0. 0*^00 

8.72 

208.72 

14.338 

0.64 

27.26 

4773. 

0.0006 

0.0600 

12.86 

212.86 

14.694 

0.97 

40.93 

4773. 

0.0008 

0.0600 

16.91 

216.91 

15.053 

1.29 

54.72 

4092. 

0.0010 

0.0600 

18.58 

218.58 

15.725 

1.39 

58.85 

3410. 

0.0012 

0.0600 

19.21 

219.21 

16.481 

1.39 

59.08 

2728. 

0.0014 

0.0600 

18.73 

218.73 

17.358 

1.31 

55.38 

2046. 

0.0016 

0.0600 

17.06 

217.06 

18.428 

1.12 

47.68 

1364. 

0.0018 

0.0600 

14.00 

214.00 

19.859 

0.85 

35.92 

682. 

0.0020 

0.0600 

9.15 

209.15 

22.210 

0.47 

20.05 

Figure  42.  Typical  Computer  Output  for  Calculating  Thermal  Power 
Generated  by  Gear  Teeth 


TABLE  5.  HEAT  GENERATED  BY  COMPONENTS 
MESHES  (GEARS,  FIGURE  40)  IDENTIFICATION 


6 

66.1 

396.6 

UP 

Sun-Planet 

4 

120.7 

482.8 

LP 

Sun-Planet 

6 

25.7 

154.2 

UP 

Planet 

-Ring 

4 

32.6 

130.4 

LP 

Planet 

-Ring 

1 

1145 

Spiral  Bevel 

2309.05 

BTU/MIN 

BEARINGS 

(FIGURE  40) 

FRICTION 

VISCOUS 

TOTAL 

TYPE 

Pinion  No 

. 16 

34.8 

23.8 

58.6 

Roller 

Pinion  No 

. 15 

109.5 

25.6 

134.9 

Ball 

Pinion  No 

. 13 

128.0 

40.7 

168.7 

Roller 

Lower  Sun 

No. 

17A 

15.2 

5.8 

21.0 

Upper  Ball 

Lower  Sun 

No. 

17B 

.14 

5.8 

5.9 

Lower  Ball 

Lower  Sun 

No. 

9 

55.6 

19.7 

75.3 

Roller 

LP  No.  10 

411.6 

18.0 

429.6 

Roller 

UP  No.  4 

409.8 

5.0 

414.8 

Roller 

Rotor  No. 

1 

21.2 

1.1 

22.3 

Ball 

Rotor  No . 

2 

8.7 

.6 

9.3 

Roller 

1194.5 

146.1 

1340.4 

BTU/MIN 

TOTAL 

3649  4 

= 86.03  HP 

TOTAL  HP 

= 

3600.,  % Lost  = 2.39 

EFFICIENCY 

ZS 

97.6% 

TOTAL  HEAT  GENERATED  = 4262.6  • 60  = 218964. 


The  heat  generated  by  all  the  bearings  as  well  as 
other  pertinent  data  is  shown  in  Figure  43.  In  order 
to  obtain  a comparative  relationship,  an  assessment  of 
the  heat  capacity  of  the  transmission  oil  was  made. 

The  following  oil- in/oil-out  temperatures  shown  in  Table  6 
were  obtained  by  experiment  (Reference  12) . 
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r 


9.33BTU/MIN 

PITCH  DIAMETER  =7.48  IN. 


BNG  SHAFT  FORCE  = 23,803  LB 


22.3BTU/MIN 

PITCH  DIAMETER  =9.84  IN. 

BNG  SHAFT  FORCE 
RADIAL  = 20,347  LB 
THRUST  = 26,889  LB 
BALL  DIAMETER  = 1.625  IN. 

NO.  OF  BALLS  = 17 
CONTACT  ANGLE  = 26  DEGREES 


.PITCH  DIAMETER  =4.48  IN. 

I BNG  SHAFT  FORCE  = 32,222  LV 
6 PAIRS 

414.75  BTU/MIN 


106  TEETH 
(RING) 


PITCH  DIAMETER 
= 4.71  IN. 

BNG SHAFT  FORCE 
= 14,972  LB 
4 PAIRS 
429.6  BTU/MIN 


58.6  BTU/MIN 


^=5^.  ./i 

168.7  BTU/MIN 

6 BTU/MIN  BNG  SHAFT  FORCE  = 12,881  LB 
21  BTU/MIN  «TCH  DIAMETER  = 6.6  IN. 


0 BTU/MIN 


HP 


(IN.  LBS)  (RPM) 
63,025 


Figure  43.  Bearing  Heat  Generation  CH-47  Forward  Transmission  100% 
Torque,  3600  HP 
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TABLE  6.  OIL- IN/OIL-OUT  TEMPERATURES 


T 

IN 

“F 

T 

OUT 

“F 

AT 

op 

160 

185 

24.81 

25 

Baseline 

256 

286 

24.81 

30 

Red  Line 

348 

398 

24.81 

50 

100%  Cooler  Bypass 
at  40%  Torque 

Using  AT  = 30°F  from  previous  experience  and  from 
the  above  information,  let: 

q - = coefficients  of  friction 
^cf 

Then,  q = WC  (AT)  Btu  for  oil  cooling, 

hr 

Where  C = specific  heat 

Then,  C = 0.485  Btu/lb  ®F  for  MIL-L-7808  oil 
and  W = oil  flow,  Ib/hr 

Hence,  W = gal  (24.)  * 231  in. ^ * 60  min  * .0347  lb  = 

inin  gal.  3 

^ in.  ^ 

= 11543  lb 
hr 

and  q . = (11543)  (.485)  (30°F)  = 167,945  Btu 

hr 

2 

Assuming  a surface  area  of  29  FT  , the  following  major 
types  of  cooling  occur  (Figure  44) : 

OIL  COOLING  93.5%, 

AIR  COOLING  1.7%, 

RADIATION  4.8%. 

^cf  ” 167,945  Btu/hr  = 93.5%  ‘of  total,  the  trans- 

mission total  heat  generated  = 167, 945/. 935  = 179,620 
Btu/hr . 
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From  the  calculation  of  heat  generated  shov\m  in 
Table  5,  218,964  Btu/hr  was  obtained.  Thus,  this 
heat  loss  is  apportioned  as  follows: 


AIR  COOLING  3,723  Btu/hr  (Natural  Convection) 

RADIATION  10,511  Btu/hr 

OIL  COOLING  204,730  Btu/hr 


218,964  Btu/hr 

If  a .5%  power  loss  had  been  assumed  for  the  spiral 
bevel  mesh,  one  would  have  obtained  196,044  Btu/hr 
(compared  to  179,620)  or  9%  error.  This  gives  a 97.9% 
efficient  transmission.  Reference  12  indicates  a 
98.6%  efficient  transmission.  This  again  is  good 
correlation . 

The  conceptual  design  analysis  of  a heat  generating 
bearing  race  and  gear  with  oil  cooling  was  conducted 
using  the  heat  transfer  capabilities  of  NASTRAN.  The 
model  and  resulting  temperatures  are  summarized  in 
Figure  45.  Using  a similar  procedure,  a thermal  study 
using  the  thermal  power  of  the  gear  meshes  and  bearings 
as  input  was  made  for  the  CH-47C  forward  transmission 
for  the  baseline  configuration  (100%  torque) . The 
resulting  temperatures  were  correlated  with  the  thermal 
mapping  test  results  (oil-out,  180°F  + 10%,  100% 
torque)  for  the  steady  state  temperature  distribution. 
Correlation  of  this  predicted  and  test  data  provided 
confidence  in  the  analytical  techniques.  This  predicted 
and  measured  housing  temperature  distribution  could 
then  be  utilized  as  input  to  the  housing  thermal  model, 
and  the  housing  thermal  distortion  and  stress  for  these 
temperatures  could  be  calculated. 
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BEARING  RACE 


Km„  = 88.5  BTU  1 HR  1 FT  =.123  BTU 

HR/FT/®F  60  MIN  12  IN.  MIN/IN./‘’F 


Kstee  L = 26.4/60*  1 2)  = .0367  BTU 

MIN/IN./°F 


MODEL  AND  THERMAL  CONDUCTIVITIES 


IIN. 


RESULTS 


INPUT  OIL  (.60,  SPECIFIED) 

= 150“F 

OIL  (60,  CALCULATED) 
= 193“F 

OUTPUT  OIL  (CALCULATED, 

70)  = 204“ 

OIL  (CALCULATED, 
80)  = 193“ 

INPUT  OIL  (GEAR  MESH, 
SPECIFIED)  = 150“ 


Figure  45.  Conceptual  Drawing  and  Sample  Calculations 
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TABLE  7.  MEASUREMENTS  OF  THERMAL  GROWTHS 
(CASE  ELEMENTS) 


DIMENSIONS 
OR  DIAMETERS* 


LOCATIONS* 


D1  (Dia.) 

D2  (Dia.) 

D3  (Dia.) 

d4  (Dia.) 

D5  (Dia.) 

Rotor  Shaft  Radial  Bearing  Housing 

06  (Dia.) 

Pump  Housing 

D7  (Dia.) 

Sync  Shaft 

06  (Dia.) 
D9  (Dia.) 
XI  (Dim.) 


Zl  (Dim.) 
Z2,  1 (Dim.) 
Z2,  2 (Dim.) 
Z3  (Dim.) 


Sync  Shaft  Co\4>lin9 

Point  thru  Rib  at  Cone  0 S/B  and  £ 
Shaft  at  D7  Diameter 

Diagonal : ^ S/B  Cone  and  Sump  Flange 
at  Bearing  C 

Between  B and  C Bearings  on  ^ Shaft 
and  Cone  £ S/B 

S/B  Cone  Center  and  Plane  2 

Outer  Points 

Inner  Points 

Between  Planes  4 and  5 


TABLE  8.  DIMENSIONAL  GROWTH  PARAMETER  EVALUATION 


DIMENSION 

IDENTITY 

D1 

(Mag.) 

D2 

(Steel) 

D3 

D4 

See  Note 

D5 

(Aim.  Aly.) 

D6 

(Alum.) 

D7 

(Mag.) 

D8 

(Mag.) 

D9 

(Steel) 

XI 

(Nag.) 

X2 

(Mag.) 

X3 

(Mag.) 

Zl 

(Hag.) 

22, 

, 1 (Steel) 

Z2, 

, 2 (Steel) 

Z3 

(Aim.) 

1^  (inches)  at 

Th  Hot 

1 Ic  TcCF) 

ih 

Th 

24.860 

73 

24.900 

190 

24.826 

73 

24.642 

190 

NOT  MEASURED 

24.877 

73 

24.940 

190 

10.500 

73 

10.516 

215 

5.579 

73 

5.595 

215 

12.501 

73 

12.533 

230 

12.444 

73 

12.470 

230 

2.643 

73 

2.645 

180 

15.520 

73 

15.556 

260 

10.638 

73 

10.718 

250 

9.690 

73 

9.710 

250 

7.864 

73 

7.899 

240 

4.275 

73 

4.261 

220 

3.243 

73 

3.237 

220 

15.634 

73 

15.682 

210 

4T 

T-Th-Tc 

(•F) 


COEFFICIENT  OF 
EXPANSION 

k X 106 


To  further  investigate  the  effects  of  temperature  upon 
a transmission  housing,  the  thermal  distortions  and 
stresses  were  calculated  using  the  NASTRAN  model  and 
thermal  map  data  from  Reference  12  for  the  following 
conditions : 


OIL-OUT 


185^F  + 10% 
286°F.  + 10% 


400  F + 10% 


TORQUE  (%) 
100 
100 
100 


The  temperature  within  each  oil-out  condition  was 
independent  of  torque.  The  NASTRAN  static  analysis 
(Rigid  Format  1)  was  used  to  calculate  the  thermal 
distortions  and  stresses,  and  dimensional  stability  of 
critical  housing  points  and  misalignment  effects  were 
evaluated . 

The  computer-generated  plot  in  Figure  47  shows  the 
regions  of  the  housing  where  it  interfaces  with  the 
bearings.  The  vectors  plotted  indicate  the  displace- 
ments at  each  grid  point  due  to  the  applied  temperature 
distribution  for  the  185°F  oil-out  condition.  By 
evaluating  the  distortion  of  the  bearing  interface  at 
each  end  of  the  respective  supporting  shafts  individu- 
ally and  then  evaluating  the  relative  distortion 
between  the  shaft  ends,  the  thermally  induced  misalign- 
ment of  the  pinion  shaft  and  the  bevel/sun  gear  shaft 
was  calculated.  By  comparing  the  relative  misalignment 
between  the  pinion  and  gear  shafts,  the  overall  effect 
of  temperature  upon  the  gear  mesh  alignment  may  be 
assessed. 


A NASTRAN  post-processor  program  uses  the  grid  point 
displacement  and  geometry  data  to  calculate  the  induced 
misalignments.  A sample  output  for  the  IBS^F  case  is 
provided  as  Appendix  C.  This  program  indicated  that 
the  induced  slopes  of  the  pinion  and  bevel/sun  shafts 
are  .0003  in. /in.  and  .0004  in. /in.,  respectively  (Figure 
48) . Also,  the  displacements  at  the  pinion  and  bevel 
gear  pitch  diameters  are  .005  and  .008  inch,  respectively. 
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VECTORS  INDICATE 


Figure  47.  Induced  Displacements  at  Housing/Bearing  Interface  Due  to 
Temperature  — 13^F  f85°C)  Oil  Out 


TYPICAL 

BEARING 

LOCATION 


DISPLACEMENT  OF  GEAR  ° 0.020  CM  (0.008  IN.) 


Figure  48.  Displacement  of  Internal  Components  Due  to  Imposed  Temperature 
Distribution  (Thermal  Map  Data  WSPF  (8SPCI  ) 
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At  the  400°F  condition,  the  pinion  and  sun/bevel  gear 
shaft  slopes  are  .0007  and  .0009  inch,  respectively; 
the  displacements  at  the  pitch  diameters  are  .013  and 
.015  inch  respectively.  Depending  upon  the  type  of 
shaft  support  bearings,  shaft  slopes  of  these  magni- 
tudes can  be  detrimental  to  bearing  performance.  Sim- 
ilarly, the  displacements  at  the  gear  mesh  point  are 
undesirable,  but  they  must  be  further  studied  to  estab- 
lish a quantitative  effect  on  gear  performance.  The 
thermal  deflections  and  stresses  are  summarized  in 
Table  9 . 


TABLE  9.  THERMAL  DEFLECTION  AND  STRESS  SUMMARY 


SHAFT  SLOPE  (IN. /IN.)  MESH  DISPLACEMENT  (IN.) 
LOAD  CONDITION  


(°F) 

PINION 

GEAR 

PINION 

GEAR 

185 

.0003 

.0004 

.0051 

.0080 

286 

.0005 

.0007 

.0087 

.0111 

400 

.0007 

.0009 

.0130 

.0150 

MAXIMUM  NOMINAL 

LOAD  CONDITION  (°F) STRESS  (PSI)  RANGE  (PSI) 

Uniform  Temperature  (160°F)  2000  100-  600 

185  Thermal  Map  3200  200-2500 

286  Thermal  Map 
400  Thermal  Map 
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5 • Thermal  Studies 


For  the  thermal  analysis  the  following  heat  rejection 
means  were  identified: 

• Conduction  (to  a heat  sink;  e.g.,  airframe) 

• Radiation 

• Convection 

- Natural  (free) 

- Forced  (oil) 

• Changes  of  State 

Reference  12  indicated  that  the  amount  of  heat  rejec- 
tion by  conduction,  radiation,  and  natural  convection 
was  small  in  comparison  to  the  forced  convection  due 
to  lubricant  cooling.  Little  data  is  available 
regarding  heat  absorption  due  to  change  of  state.  In 
one^test/  smoke-off  (change  of  state)  was  recorded  at 
240  F with  aggravated  smoke-off  (vaporization)  at 
300  F.  The  following  equations  were  useful  as  a 
basis  for  the  thermal  studies: 


• Radiation:  The  Stef an-Boltzmann  Law: 

q = AE.a  (T  ^ - T-^)  = AEa  ( (t,  + 460)^  - 

IT  X ^ X 

(t^  + 460)^) 

where,  q^  = heat  rejection  rate,  Btu/hr 

A = black  gray  body  (specimen  transmission) 
surface  area,  ft  = 29  ft^ 

T^  = black-gray  body  surface  temperature  (equal 
to  oil-out  temperature) , ®R 

T = temperature  of  surroundings  (test  cell  wall 
^ temperature  - ambient  air  temperature) , 560° 

e = emissivity/absorptivity  constant;  optimisti- 
cally assumed  at  0.9 

0 = Stefan-Boltzmann  constant  assumed  at  0.173 

X 10"® 
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which  becomes 


q = 4.5153 


X 10 


-8 


(t^  + 460)  - 


(t^  + 460)/*) 


• Convection; 

Natural  Convection 

^cn  = ^ (^1  “ ^2^ 


where  a = heat  rejection,  Btu/hr 
^cn 

2 

A = specimen  transmission  surface  area  = 29  ft 

1/3 

h = 0.22  (t^  - t^) 

D = specimen  transmission  average  diameter  = 2 ft 

= specimen  transmission  surface  temperature 
^ °F  or  °R 

t2  = ambient  air  temperature,  same  degree  unit 
which  becomes 


q 


cn 


6.38  (t^  - t^) 

Forced  Convention  (Oil) 


where  W = oil  flow,  8775  Ib/hr 
C =0.485  Btu/lb/°F 

t^  = specimen  transmission  oil-out  temperature,  °f 
t^  = specimen  transmission  oil-in  temperature,  °F 

Based  on  test  results  (ti  - t2)  = At  = 30°F,  and  q^f  = 
8775  X 0.485  x 30  = 127,676  Btu/hr  heat  rejected  by  oil 
for  100  percent  CH-47C  transmission  power. 
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The  same  test  indicated  an  efficiency  for  the  CH-47 
forward  rotor  transmission  of  98.6  percent.  At  100 
percent  CH-47C  forward  transmission  power  (3600  SHP) , 
the  heat  rejection  rate  is  128,000  Btu/hr  (2134  Btu/ 
min).  For  75  percent  and  50  percent  powers,  the  heat 
rejection  rates  are  96,000  Btu/hr  (1600  Btu/min)  and 
64,000  Btu/hr  (1067  Btu/min)  respectively.  The  family 
of  curves  in  Figure  49  indicates  the  three  specimen 
power  levels  and  heat  rejection  rate  versus  average 
case  temperature.  When  the  oil  cooler  is  fully  by- 
passed and  heat  rejection  is  accomplished  by  natural 
convection  and  radiation  alone,  the  case  temperatures 
theoretically  would  reach  those  points  indicated  by 
the  intersection  of  the  curves  with  the  zero  percent 
ordinate.  However,  an  analysis  of  the  aft  transmission 
test  does  not  bear  this  out.  With  the  evidence  of 
smoke-off  beginning  at  240®F  and  terminal  temperature 
at  320Of,  it  appears  that  other  means  of  heat  rejection 
are  operating.  Otherwise,  the  case  temperature  would 
reach  approximately  4300f.  The  straight  line  inter- 
section, both  the  lowest  curve  and  240 °F  and  the  zero 
percent  ordinate  at  330°F,  represents  a cutoff  of  case 
temperature  rise  to  match  test  experience. 

After  conducting  the  initial  model  development  and 
validation,  studies  were  performed  to  evaluate  a closed 
transmission  system  with  no  external  cooler.  Such 
factors  as  air  flow  and  use  of  cooling  fins  were  con- 
sidered. The  housing  model  was  modified  to  include 
cooling  fins  of  various  configuration  and  total  area. 

Fin  cooling  assuming  an  oil-out  temperature  of  200°F  v, 

was  investigated. 


Assuming  the  conservative  .75%  power  loss  figure  for 
the  spiral  bevel  mesh,  a study  of  forced  convection 
(air)  to  eliminate  the  oil-cooling  lines  (204,730  Btu/ 
hr)  was  conducted.  For  the  conceptual  study,  assume 
that  a transmission  housing  is  represented  by  a 2-foot- 
diameter  by  3-foot-high  cylinder.  The  Nusselt  number 
for  the  flow  perpendicular  to  the  longitudinal  axis, 
from  Reference  19,  is 


Nu^  = ,43  + C (Re^) 


m 


where  Re^,  the  Reynold's  number,  = 


u d 
m 


and 


u^  = mean  air  velocity  (ft/sec) 
d = cylinder  diameter  (ft) 

2 

^ = kinematic  viscosity  of  air  (ft  /sec) 


19.  Eckert  and  Drake,  HEAT  AND  MASS  TRANSFER,  McGraw-Hill, 
1959. 
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2,134  BTU/MIN 
128,000  BTU/HR 
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Theoretical  Temperature  Projection  for  Given  Heat  Rejection  Rates  and  Power 
Levels  (Based  on  CH-47  Forward  Transmission  Efficiency  of  98.6  Percent) 
(Reference  12) 


In  order  to  determine  Uu,  for  an  unfinned  transmission, 
assume  the  ambient  temperature  to  be  100°F. 

^ = 18.1  X 10“^  FT^/SEC, 

K = .01565  BTU/HR  FT  °F  (Thermal  Conductivity), 


C = .0239 


m = .805 


(From  Reference  19) 


The  Nusselt  number  may  also  be  expressed  as 


XT  _ 'hd 
*^^d  K 


2h 

.01565 


where  h = film  heat  transfer  coefficient  (Btu/hr  ft^  °F)  . 
Then  h = Nu^  = -Q^565  [.43  + C(Re,)"‘] 


.015651 


!"  .43  + 


,0239 


f ^ ^ 

ll8.1  X 


h = .0034  + .3363  (u^^^) 


.805 


Assuming  conservatively  that  the  effective  area  of  the 
cylinder  is  the  frontal  area  (Figure  50) , then, 

Q = I'  ^eff  ^ (9.42)  (200°F  - 100°F)  Btu/hr 

Here  the  cylinder  surface  temperature  is  assumed  to 
be  200°F,  and  the  ambient  temperature  = 100°F. 

This  gives 

h = 204748  = 217.4 
942 

Hence, 

= 3101.0  ft/sec  a 2114  mph 

The  Reynold's  number  in  Figure  51  is  low  (1260). 

Here  the  Reynold's  number  is  7,144,750;  at  high 
Reynold's  numbers,  the  film  heat  transfer  increases 
on  the  back  of  the  cylinder.  In  fact,  at  Re  * 50,000, 
h is  the  same  on  the  front  and  back. 


I 

CYLINDER 


Figure  50.  Effective  Area  of  a Finned  Cylindrical  Surface 


Figure  5 1.  Isotherms  Around  a Cylinder  Cooled  by  a Fluid  Flowing  Normal  to  its 
Axis,  as  Revealed  by  an  Interference  Photograph.  Re  = 1,260 
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Repeating  the  calculations  gives 


u = 1308.7  ft/sec  = 892.2  mph 
m 


If  the  shape  of  the  transmission  is  considered 
(Table  10) , then 


u = log 
m ^ 


4.142  ~ log  C - 4.04 
m 


ft  /sec 


yields ; 


CONFIGURATION 

(ft/sec) 

1 

4273 

2 

13058 

3 

6314 

4 

1328 

TABLE  10.  COEFFICIENTS  FOR  CALCULATION  OF 

HEAT  TRANSFER  FROM  CYLINDERS  WITH 
DIFFERENT  CROSS  SECTIONS  TO  AN  AIR 
FLOW  NORMAL  TO  THEIR  AXES  (Reference  19  ) 


1 

1 Cross  section 

Res 

c 

m 

Configuration 

- □ 

5,000-100,000 

0.0921 

0.67$ 

1 

o 
-*  o 

5,000-100,000 

0.222 

0.588 

2 

5,000-100,000 

0.138 

0.638 

3 

5,000-  19,500 

0.144 

0.638 

- 0 

19,500-100,000 



0.0347 

0.782 

4 
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Such  air  velocities  are  not  only  impossible,  but  they 
would  cause  a heat  increase  by  converting  internal  skin 
friction  into  heat.  As  a second  resort,  fins  may  be 
added.  In  fact,  they  can  be  combined  with  a forced 
air  flow. 

A conceptual  design  is  made  for  the  previous  cylinder 
to  find  the  number  of  fins  2 inches  long  and  1 inch  thick 
that  would  be  required  in  a forced  air  flow  of  50  ft/sec 
(34  mph) . 

(with  fins)  = 2 tt  y b K ^ ^ (^T)  Z = 

Heat  Loss  Per  Fin,  BTU/HR 
where : 

^1  ^^2^  ^ ^o  ^^1^  ^ ^ ^1  ^^2^  ^ ^o  ^^l"^  ^ 

2h  , 

Kb 


= 100  Btu/hr  ft  “F  (Conductivity  of  Magnesium) 


^1 


Y 2 


= 12  in. 

= 14  in. 

= 1 in.  (see  Figure  50) 


Let 


"^AMBIENT 

^SURFACE 


V 


= 50  ft/sec  (Forced  Convection) , 


^1  ' 


I , K-  and  K are  Bessel  functions, 
o 1 o 


Now 


Re  = (50)  (2)  = 552,486 

° .000181 

Nu,  = .43  + .0239  (552486)’®°^  = 2h 

= 1002.66  = 2h 

.01565 

h = 7.845  Btu/hr  ft^  °F 


Using  the  previous  method  of  calculation,  the  forced 
convective  cooling  of  the  cylinder  with  no  fins  is: 


= (7.846)  (18.84)  (100)  = 14781.9  BTU/HR 

With  the  addition  of  one  fin, 

3 = 12)  (7.846)  (12)  = 1.883,  IT  = 1.3722, 

(100)  (1") 

Y,  J7=  12  (1.3722)  = 1.3722,  y 1.6009 

’ 12  2-  ^ 

= (2  7r)  (!•)  (1/12)  (100)  (1.3722)  (100)  Z = 

7184.8  Z per  fin.  Now 

I 


(1.6009)  = 1.086 

^1 

(1.3722)  = .3368 

'V, 

(1.3722)  = .862 

^1 

(1.6009)  = .2406 

(1.3722)  = 1.531 

K 

O 

(1.3722)  = .2541 

Z = (1.086)  (.3368)  - ( .2406)  ( .862)  = .1584  = 2458 

(1.086) (.2541)  + (.2406) (1.531)  .6443 


Hence,  = 1766  Btu/hr  per  fin 

In  order  to  eliminate  the  oil  cooling  of  204,748  Btu/hr, 
it  would  require  116  fins.  For  a 3-foot  high,  1-inch- 
thick  fin,  this  is  impossible.  Hence,  the  air  velocity 
of  50  ft/sec  (34  mph)  must  be  increased.  On  the  other 
hand,  the  thickness  of  the  fins  could  be  decreased. 
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If 


V = 100  ft/sec 

Re,  = (100)  (2)  = 1,104,974 

.000181 

Nu,  = .43  + .0239  (1,104,972) 


.805 


= 1751.5  = 


2h 


.01565 


h 

Q. 


13.7 

25810.8  for  no  fins 


y.JT 


= 12) 

7)  (12) 

= 3.288 

' JT 

= 1. 

81 

(100) 

(1) 

= 1.81 

iJT  = 

2.117 

(2  TT  ) 

d' 

) (1/12) 

(100) 

(1.81) 

(100) 

Z 

= 9477.1 

(2. 

117) 

= 1.78 

(1 

.81) 

= 

.180 

(1. 

81) 

= 1.33 

(2 

.117) 

= 

.121 

^o 

81) 

= 2.00 

""o 

.81) 

= 

.144 

Z = (1.78)  (.180)  - (.180)  (1.33)  = .1595  = .320 
(1.78)  (.144)  + (.121)  (2.00)  .4983 

= 3033  Btu/hr  per  fin.  This  gives  68  fins  (b  = 1 ft) 
which  is  again  impossible. 

If  V = 25  ft/sec  = 17  mph, 

= 1008  Btu/hr,  203  fins,  1 in.  thick,  which  is  clearly 
impractical. 


In  summary; 
V (ft/sec) 

25 

50 

100 


FIN  COOLING 
(Btu/hr) 

1008 

1766 

3033 


NUMBER 
OF  FINS 

203 

116 

68 


HEIGHT  (ft) 

16.9 

9.7 

5.7 
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Now,  let  b = .5  in.  (fin  thickness)  and  let  V = 50  ft/sec. 

Re,  = (50)  (2)  = 552,486  as  before 
° .000181 


Hence , 

h = 7.846 

e = (2)  (7.846)  (12)  = 3.766,  ^IT  = 1.94  , 


(100) 

(.5) 

= 1.94  , 

2.263 

Hence,  = (2  tt  ) (1) 

(.5)  (100)  (1.94) 
12 

(100)  z = 

5078.9 

(2.263) 

= 2.035 

(1.940) 

= .1527 

(1.940) 

= 1.5086 

(2.263) 

= .1003 

(1.940) 

= 2.1926 

Kq  (1.94) 

= .1235 

2 = (2.035)  (.1527) 

- (.1003) 

(1.5086) 

= .1594  = 

.3383 

(2.035)  (.1235) 

+ (.1003) 

(2.1926) 

.4712 

= (.3383)  (5078.9)  = 1718  Btu/hr 


This  gives  119  fins,  which  is  again  impractical.  Some 
results  from  a computer  program  to  predict  fin  cooling 
with  forced  air  flow  across  a cylinder  (Appendix  D) 
are  summarized  in  Table  11. 


TABLE  11.  RESULTS  OF  FIN  COOLING  STUDY 


V 

(ft/sec) 

LENGTH 
(in. ) 

THICKNESS 
(in. ) 

COOLING 

(Btu/hr) 

NUMBER 
OF  FINS 

HEIGHT  OF 
FINS  (ft) 

50. 

. 2. 

.1 

1505 

136 

1.1 

V 

50. 

2. 

.2 

1629 

126 

2.1 

50. 

2. 

.3 

1676 

122 

3.05*- 

50. 

2. 

.5 

1716 

Impractical* 

50. 

2. 

1.0 

1747 

Impractical* 

20. 

2. 

.1 

782 

262 

2.2 

20. 

2. 

.2 

815 

251 

Impractical* 

20. 

2. 

.3 

20. 

2. 

.5 

20. 

2. 

1.0 

♦The  accumulative  height  of  the  fins  must  be  less  than 
the  height  of  the  cylinder  (3  feet) . 


Figure  52  shows  the  feasibility  of  using  thin  fins  with 
forced  air  convection  to  replace  external  oil  cooling 
(V  = 50  FT/SEC) . A fan  delivering  a nominal  20  ft/sec  with 
262  2-inch  fins,  .1  inch  thick,  is  in  the  realm  of  possi- 
bility also.  It  is  noteworthy  that  fin  cooling  is  insensi- 
tive to  fin  thickness  above  about  0.6  inch. 
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The  relative  importance  of  forced  convection  cooling 
(oil),  natural  convection  (air  cooling),  and  radiation 
was  determined  from  Reference  12.  For  example, 

(1)  For  oil  cooling: 

= WC  AT  = 167,800  Btu/hr 

W = Oil  Flow  = 11,532  Ib/hr  = 24  gal/min, 

C = .485  BTU/LB/°F  for  MIL-L-7808  Oil, 

A T = (Oil-Out  Temperature)  - (Oil-In 
Temperature)  = 300f 

(2)  For  natural  convection: 

= .22  (AT)"^/^  (A) 

= .22  (100)^^^  (29)  = 2961.3  Btu/hr 

where  AT  is  the  difference  between  the  ambient 
and  surface  temperatures  (200OF  - lOQOF) , and 

A = Surface  Area  = 29  ft^ 

The  above  .formula  is  valid  for  a cylinder  with  a 
2-foot  diameter. 

(3)  For  radiation: 

= Aea  |jt^  + 460°f)'^  - {t^  + 460°f)'^ 

= 8564  Btu/hr 

t,  = °F  of  Oil-Out  = 200°  (Assumed  Surface 
Temperature) 

t2  = °F  of  Ambient  Air  (100°F) 

= Stefan-Boltzman  Constant  (Black-Grey 
Area) 

= .173  X 10~®  Btu/hr/ft2/°R^ 
e = Emissivity /Absorptivity  Constant  = .9 
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Therefore,  the  oil  cooling  is  the  most  significant 
means  of  heat  removal;  and  on  a percentage  basis: 

Oil  Cooling  = 93.5% 

Air  Cooling  = 1.7% 

Radiation  = 4.8% 

A thermal  analysis  of  the  cooling  fin  conceptual  de- 
sign was  completed  using  NASTRAN  for  three  conditions: 
Temperature  specified  at  root,  heat  input  at  root, 
and  a combination  of  these  two  conditions.  The  theo- 
retical and  computer  solutions  were  obtained.  The 
model  consisted  of  a metal  fin  of  triangular  cross- 
section  attached  to  a plane  surface  to  help  carry  off 
heat  for  the  latter  (Figure  53) . Further  detail 
including  the  computer  output  is  provided  in  Appendix  D. 

A heat  transfer  analysis  was  made  of  the  entire  CH-47C 
forward  transmission  case.  Included  were  forced 
convective  oil  cooling,  thermal  conductivity  and  heat 
generated  by  the  input  pinion/sun  gear  bearings.  The 
NASTRAN  run  was  made  in  order  to  print  and  punch  out 
the  temperature  distribution.  This  was  followed  by  a 
thermal  stress/distortion  analysis.  The  following 
was  determined; 

(1)  NASTRAN  Level  15.5.1  (Reference  20)  will  not 
worJc  when  starting  from  a structural  check- 
point tape.  It  is  necessary  to  punch  out  the 
input  data  using  "ECH0  = PUNCH,  S0RT"  in  the 
case  control,  and  "ALTER  4,138",  "ENDALTER" , 
in  the  executive  control. 

(2)  The  cost  of  a heat  transfer  analysis  is  only  8% 
that  of  a full  stress  analysis  and  is  only  2% 
that  of  a dynamic  analysis. 

(3)  Computer  core  requirements  are  half  those  of  a 
statics  or  dynamics  solution. 

(4)  Correlation  (Figure  54)  yielded  somewhat  low 
values  where  oil  cooling  was  applied. 

(5)  Forced  convective  oil  cooling  in  NASTRAN  is 
limited.  For  a given  oil  flow,  downstream  temper- 
atures of  the  oil  are  the  same  as  the  upstream 
temperatures  (Reference  21) . 


20.  NASTRAN  User's  Manual  Level  15.5,  1976 

21.  Thornton,  E.,  and  Wieting,  A.,  COMPARISON  OF  NASTRAN  AND 
MITAS  THERMAL  ANALYSIS  OF  A CONVECTIVELY  COOLED  STRUCTURE, 
1975  Fourth  NASTRAN  User's  Experiences. 
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NOTE:  ALL  NUMBERS  SHOWN 
ARE  TEMPERATURES,!,  IN 
DEGREES  FAHRENHEIT  (AT 
SHUTDOWN). 

TEST  STAND  (HEAT  SINK)  = 
96,  98, 118 

TEST  CELL  AIR  =89,  95,98 


BEARING  2 
160,164 


BEARING  1 
162, 164 


157, 156 


162, 160  AT  90“  ■ 


158  AT  270“,  156 


(6)  Thermal  stresses  calculated  using  the  punched 

output  were  similar  to  those  obtained  previously 
and  shown  in  Table  9.  The  maximums  were  +2000 
psi . 

After  conducting  the  thermal  analysis,  an  assessment 
of  the  effects  of  the  dimensional  instability  of  the 
critical  housing  points  on  the  bevel  gear  mesh  was 
attempted.  The  exact  effect  of  housing  distortion  on 
the  load  capacity  of  the  bevel  gear  mesh  is  an  ex- 
tremely complex  subject  and  has  been  the  topic  of 
discussion  among  gear  experts  for  some  time.  Axial 
movement,  lateral  movement,  and  induced  misalignment, 
of  both  the  pinion  and  gear  members  must  be  determined, 
then  the  results  of  all  of  these  changes  for  each 
individual  member  must  be  evaluated  to  determine  the 
relative  effect  on  each  member.  Further  complications 
are  introduced  by  the  changes  in  backlash  and  clear- 
ances due  to  the  thermal  growth  of  the  gears  themselves. 
In  order  to  provide  some  indication  of  the  effect  of 
temperature  distortion  on  the  gear  mesh,  the  backlash 
and  mesh  pattern  (i.e.,  load  distribution)  were  assessed 

The  transmission  housing  flange  which  serves  as  the 
mounting  surface  for  the  input  pinion  cartridge  experi- 
enced an  outward  thermal  growth  of  0.029  to  0.032  inch 
when  the  housing  was  subjected  to  the  185°F  temperature 
distribution  from  Reference  12.  Considering  the  bearing 
stack-up,  mounting,  and  method  of  transfer  of  axial  load 
from  housing  to  gear  shaft,  the  result  is  an  axial  out- 
ward movement  of  the  bevel  pinion  of  0.019  to  0.022  inch 
Using  the  backlash  versus  axial  movement  plot  in  Figure 
14,  this  axial  movement  causes  an  increase  in  the  back- 
lash of  the  gear  mesh  of  about  0.007  to  0.008  inch, 
which  is  on  the  order  of  a 100%  increase. 


Actual  bevel  gear  mesh  patterns  corresponding  to 
varying  axial  positions  of  the  input  pinion  are  shown 
in  Figure  55.  It  is  evident  from  these  patterns  that 
the  axial  movement  has  a significant  effect  on  the 
mesh  pattern.  Figure  5 indicates  that  even  a properly 
pat'erned  ge?  ' mesh  with  a full  contact  pattern  experi- 
ences a peaked  stress  distribution  which  exceeds 
theoretical  predictions  based  on  an  assumed  uniform 
stress  distribution.  Any  deviation  from  this  correct 
pattern  will  aggravate  this  stress  peak  and  could  lead 
to  premature  failure. 


NO-LOAD  (BENCHI  PATTERN 


FULL-LOAD  PATTERN 


"ZERO"  POSITION 


C 


■ZERO"  POSITION 


V//LLLLUJ^ 

0.022  INCH  AXIAL  MOVEMENT 
OF  PINION  OUT  OF  MESH 


0.004  INCH  AXIAL  MOVEMENT 
OF  PINION  OUT  OF  MESH 


WHA 


CONTACT  PATTERN 


Figure  55.  Bevel  Gear  Mesh  Patterns  at  Different  Axial  Positions  of  Pinion 
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IV.  STATIC  AND  DYNAMIC  TRANSMISSION  STRESS  ANALYSIS  INCLUDING 
IN-FLIGHT  DYNAMIC  ROTOR  LOADS  AND  LOAD  PATH  DETERMINATION 


An  analytical  method  for  accurately  defining  the  stress 
distribution  and  load  paths  in  a helicopter  transmission 
housing  has  been  investigated.  Previously,  the  designer 
had  little  guidance  for  selection  of  the  design  with  best 
structural  efficiency.  With  continually  increasing  power 
requirements,  the  weight  penalties  imposed  by  nonoptimum 
structural  configuration  may  be  significant. 

Because  of  the  many  functions  performed  by  a transmission 
housing  and  its  complex  geometry,  analysis  is  difficult. 
Using  a NASTRAN  finite  element  housing  model,  however,  a 
static/dynamic  stress  analysis  may  be  conducted  and 
structural  deformations  predicted  (Figure  56) . Further- 
more, methods  for  structural  optimization  can  be  applied 
to  reduce  stress,  vibration  and  weight. 

The  work  documented  herein  includes  static  and  dynamic 
stress  analyses  and  considers  rotor  loads,  g-loads,  and 
steady-state  gear/bearing  loads  imposed  on  a transmis- 
sion housing.  The  NASTRAN  computer  program  can  also 
readily  handle  fabricated  and/or  composite  structures 
as  well  as  conventional  cast  metal  materials.  This  is 
accomplished  through  the  ability  of  NASTRAN  to  use  as 
input  a 6 X 6 material  property  matrix  and  an  orienting 
angle  for  each  element  to  define  the  direction  of  the 
property  orientation.  The  flow  diagram  for  the  stress 
analysis  is  shown  in  Figure  57. 

By  applying  representative  loads  to  the  housing  model, 
the  stress  distribution  throughout  the  housing  can  be 
calculated  for  varying  conditions.  The  static  and  dynamic 
stress  thus  calculated  can  be  superimposed  upon  the  thermal 
stress  distribution  to  provide  an  accurate  overall  picture 
of  both  the  steady-state  and  the  time-dependent  (fatigue 
producing)  stresses  occurring  in  the  housing  of  an  operat- 
ing transmission.  From  this  combined  stress  distribution, 
the  structural  load  paths  can  be  identified  and  the  struc- 
tural portions  of  the  housing  segregated  from  the  non- 
structural  portions.  By  utilization  of  structural  optimi- 
zation methods,  wall  thickness  changes  can  be  recommended 
and  weight  reduction  evaluated. 

Vibratory  3-per-rev  hub  loads  and  steady  Ig  loads  were 
calculated  using  a proprietary  Boeing  Vertol  rotor  loads 
analysis  computer  program  (L-02)  . A seunple  of  the  com- 
puter output  is  summarized  in  Figure  58.  The  calculated 
3-per-rev  dynamic  hub  loads  which  were  applied  to  the 
housing  model  via  the  rotor  shaft  support  bearings  are 
shown  in  Figure  59. 
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UPPER  COVER 


GRID  POINTS 
ELEMENTS 


RING  GEAR 

GRID  POINTS  = 216 
ELEMENTS  = 192 


CASE 

GRID  POINTS 
ELEMENTS 


• AUTOMATIC  GENERATION  OF  GRID 
POINTS  AND  ELEMENT  CONNECTIONS. 

■ HOUSING  MADE  OF  PLATE  ELEMENTS. 

■ METAL  OR  COMPOSITE  MATERIAL. 

• STATIC  STRESS  AND  DISTORTIONS. 

• LOAD  PATH  IDENTIFICATION. 

• STRUCTURAL  OPTIMIZATION. 

• REDUCE  STRESS  CONCENTRATIONS  BY 
CHANGING  WALL  THICKNESS  AND 
GEOMETRY. 

• SURVIVABILITY  AND  FAILSAFETY 
STUDIES. 

• VALIDATION  OF  ANALYTICAL  TOOLS 
FOR  APPLICATION  TO  ADVANCED 
TRANSMISSION  DESIGN. 


Figure  56.  Model  for  Dynamic/Static  Stress  Analyses 
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MESH  EXCITATION 
AND  COMPLIANCES 


Figure  57.  Flow  Diagram  of  NAST RAN  Stress  Analysis 
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3/Rev  Cosine  Ccxnponent  ^ 


ADVANCED  ROTOR  AEROELASTICITY  PROGRAM  L02 
FWD  RTR  72-123-01 
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Figure  58.  Sample  Rotor  Loads  Program  Output  ICH-47C  Forward  Rotor  at  243  RPM, 
50,000  LB  Gross  Weight,  150  KTS  - 3-Per-Rev  and  Steady  ( 1g)  Loads. 
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Another  Boeing  Vertol  finite  element  computer  program 
(D82)  determines  shaft  deformations  due  to  mesh  excit- 
ation and  translates  these  deformations  into  vibratory 
internal  bearing  loads  for  use  in  the  structural  analysis 
program.  These  vibratory  internal  loads  have  been  thor- 
oughly investigated  in  Reference  9 and  will  not  be  pursued 
further  here,  but  the  results  will  be  applied.  Natural 
frequencies  and  mode  shapes  for  the  housing  have  been 
calculated  and  stored  on  magnetic  tape  and  are  used  for 
both  of  the  above  dynamic  stress  analyses.  Applying  the 
vibratory  3-per-rev  hub  loads  and  the  dynamic  loads  due  to 
the  internal  components,  the  NASTRAN  vibratory  response 
analysis  (Rigid  Format  11)  was  used  to  calculate  dynamic 
response  and  stresses  of  the  housing. 

Steady-state  loads  of  two  types  were  considered:  Ig  steady 
level  flight  loads  and  transient  g-loads  for  three  flight 
maneuvers  including  the  ultimate.  All  of  these  loads, 
which  include  both  forces  and  moments,  are  applied  to  the 
housing  model  via  the  bearings.  The  three  maneuver  condi- 
tions are  calculated  similarly  to  the  inertia  relief 
capability  of  NASTRAN.  The  static  stress  distribution  for 
the  housing  using  each  of  these  load  conditions  is  calcu- 
ilated  using  NASTRAN  static  analysis  (Rigid  Format  1)  . 
Application  of  the  calculated  steady  Ig  hub  loads  is 
shown  in  Figure  60. 

Furthermore,  the  effect  of  bearing  nonlinearity  was 
investigated.  Since  the  bearings  deflect  nonlinear ly 
under  load  (Figure  61) , the  load  paths  generated  under  a 
static  loading  would  be  different  if  the  bearings  were 
assumed  to  be  linear  elastic.  This  could  be  analyzed 
using  the  nonlinear  feature  of  NASTRAN  called  piecewise 
linear  analysis  (Rigid  Format  6) . The  bearings  retaining 
the  CH-47C  forward  transmission  pinion  and  sun  gear  have 
nonlinear  stiffnesses.  The  bearings  (13  to  16,  9 and  17) 
are  shown  in  Figure  62.  Other  nonlinear  bearings  are 
shown  in  Figure  63.  These  are  bearings  1 and  2 for  the 
rotor  shaft,  the  upper  planetary  spherical  roller  bear- 
ings 4 and  the  lower  planetary  spherical  roller  bearing 
10.  The  stiffness  of  the  bearings  is  a function  of  ro- 
tor torque  and  may  be  calculated.  This  variation  in 
stiffness  as  a function  of  percentage  of  torque  for  the 
pinion-sun  gear  bearings  is  shown  in  Figures  64  and  65. 

Since  the  bearings  are  a small  portion  of  the  transmission 
housing  model,  a separate  computer  run  could  be  made 
assuming  the  bearings  to  be  linear.  The  load  paths  for 
the  linear  case  can  be  compared  to  the  nonlinear  load 
paths  to  ascertain  the  differences.  The  costs  of  the  two 
computer  runs  would  be  considered  in  the  evaluation.  It 
is  probable  that  the  nonlinearity  of  the  bearings  is  of 
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NOTE:  AS  THE  LOAD  INCREASES,  THE  RATE  OF  THE 

INCREASE  OF  DEFLECTION  DECREASES, THERE- 
FORE PRELOADING  (TDP  LINE)  TENDS  TO 
REDUCE  THE  BEARING  DEFLECTION  UNDER 
ADDITIDNAL  LOADING. 


Figure  61.  Deflection  Versus  Load  Characteristics 
for  Ball  Bearings 


INPUT  PINION  SUPPORT  BEAM 


Figure  62.  CH-47C  Forward  Transmission  Support  Bearings 
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Figure  64.  Bearing  Spring  Rate  (K)  for  CH-47C  Forward  Transmission  Input  Pinion 
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Figure  65.  Bearing  Spring  Rate  (K)  for  CH-47C  Forward  Transmission  Sun  Gear 
(Pinion  RPM:  7,460) 


no  consequence  in  static  and  maneuver  loadings.  However, 
the  nonlinearity  of  the  bearings,  which  is  a function  of 
rotor  torque,  does  affect  the  natural  frequencies  of  the 
internal  components.  These  natural  modes,  when  excited 
by  the  mesh  frequencies,  would  change  the  dynamic  loading 
on  the  case.  But  even  in  this  dynamic  situation,  it  would 
be  erroneous  to  overemphasize  the  nonlinear  stiffnesses 
since  the  bearings  act  only  in  compression,  not  in  ten- 
sion, which  is  significantly  more  nonlinear. 

Three  maneuver  conditions  were  selected  for  sample  static 
stress  analyses  of  the  CH-47C  forward  transmission  case. 
The  sign  convention  for  these  hub  loads  is  shown  in 
Figure  69  and  the  loads  are: 

• Symmetric  dive  and  pullout,  noseup  pitching 
(Figure  70) 

• Yawing  (Figure  71) 

• Recovery  from  rolling  pullout  (counterclockwise) 
(Figure  69) 

Three-per-rev  dynamic  loads  and  stresses  on  the  case  were 
determined.  Since  the  natural  frequencies  of  the  case 
start  at  about  600  Hz,  there  is  no  magnification  of  the 
12  Hz  exciting  frequency  (3/rev) . The  coupled  bending/ 
torsion  natural  frequencies  of  the  internal  components 
start  at  560  Hz;  so  once  again,  there  is  no  magnification 
of  the  loads.  Figure  55  shows  that  the  sine  and  cosine 
components  as  output  by  the  rotor  loads  computer  progrcun 
(Figure  54)  had  to  be  applied  to  the  hub  in  separate  runs. 
A bearing  computer  program  used  these  separate  load  condi- 
tions to  calculate  nonlinear  bearing  loads.  In  Figure  73, 
the  bearing  dynamic  loads  on  the  case  are  illustrated 
except  for  the  pinion  ball  bearing  in  the  X-direction 
(64  lb)  and  the  gear  duplex  ball  bearing  load  in  the  Z- 
direction  (24  lb).  However,  these  loads  were  applied  to 
the  NASTRAN  model.  The  sine  and  cosine  bearing  loads 
were  combined  into  single  polar  loads  for  the  NASTRAN 
dynamic  stress  analysis.  At  a rotor  shaft  speed  of  243 
rpm,  the  dynamic  stresses  due  tc  mesh  excitations  were 
calculated . 
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Figure  67.  Ultimate  Forward  Rotor  Head  Loads;  Symmetric  Dive  and  Pullout,  Noseup 
Pitching  (Case  4 - Flight  Condition  2B,  33,000  Lb  Gross  Weight) 


= 1,750,000  IN. -LB 
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Figure  68.  Ultimate  Forward  Rotor  Head  Loads;  Yawing  (Case  5 — Flight  Condition  68, 
33,000  LB  Gross  Weight! 


o 


115 


MESH  FREQUENCIES 
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Figure  70.  Bearing  Loads  Applied  to  Case  Due  to  Mesh  Excitations  ( lb)  Calculated  by 
Beam  Model  Dynamic  Analysis  of  Internal  Components 


■The  vibratory  loads  (exclusive  of  phasing)  are  shown 
on  Figure  70  for  the  above  four  mesh  excitations.  Other 
mesh  excitations  and  sidebands  exist,  but  microphone 
data  indicates  the  predominance  of  the  above.  Ihe  max- 
imum dynamic  stresses  indicated  by  NASTRAN's  prodigious 
output  are  approximately: 

(1)  600  psi  for  the  LPl  Excitation 

(2)  500  psi  for  the  LP2  Excitation 

(3)  200  psi  for  the  SBl  Excitation 

(4)  250  psi  for  the  LP3  Excitation 

This  is  based  on  an  assumed  equivalent  viscous  damping 
of  3%  (6%  structural) . The  values  of  these  loads  are 
contingent  on  the  natural  frequencies  of  the  internal 
components  which  also  include  the  nonlinear  bearing 
stiffnesses  and  the  gear  tooth  mesh  stiffness  (spiral 
bevel  pinion/gear) . 

Bearing  loads  for  Ig  (high-speed  level  flight)  and 
three  maneuver  conditions  were  calculated  from  hub 
loads.  The  hub  loads  for  these  conditions  are  given  in 
Figures  60,  67,  68  and  69.  All  bearing  loads  calculated 
using  these  hub  loads  as  well  as  the  bearing  number  con- 
vention are  given  in  Appendix  E. 

Maximum  stresses  for  the  static  analyses  of  the  case 
are  summarized  here: 

1.  Ig..  Steady  Flight,  + 3000  PSI 

2.  Ultimate  Maneuver,  + 15,000  to  25,000  PSI 

3.  Yawing  Maneuver,  + 20,000  PSI 

4.  Recovery  from  Rolling  Pullout,  + 5,000  PSI 

These  are  "eyeball"  numbers  and  are  isolated  areas  on 
the  case.  However,  assuming  an  allowable  of  +26,000 

psi  for  AZ91  magnesium,  there  exists  the  possibility 
of  high  stresses  on  the  case  for  certain  maneuvers. 
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Deflections  for  the  critical  shaft  support  bearing/ 
housing  interface  locations  were  evaluated.  Figure  71 
is  a schematic  of  the  bevel  pinion  and  sun/bevel  gear 
shaft  showing  the  predicted  deflections  and  slopes  due 
to  an  imposed  ultimate  load  condition  for  a magnesium 
and  steel  housing,  respectively.  The  steel  housing  was 
evaluated  in  order  to  establish  a basis  for  comparison. 
The  displacements  allowed  by  the  magnesium  housing  are 
four  to  five  times  that  of  the  steel  housing.  Table  12 
summarizes  the  critical  deflection  information. 

It  is  evident  from  Figure  72,  which  plots  deflection 
versus  torque  for  various  housing  materials,  that  the 
magnitude  of  the  housing  displacements  can  be  reduced 
substantially  by  the  use  of  stiffer  materials.  Actual 
deflection  test  data  is  shown  for  the  magnesiimi  housing. 
A steel  housing  is  apparently  very  desirable  from  the 
stiffness  aspect,  although  obviously  unacceptable  for 
aircraft  application  because  of  its  weight.  However, 
the  metal  matrix  material  provides  good  stiffness  char- 
acteristics at  only  a small  weight  penalty.  Further- 
more, the  slight  weight  increase  can  be  traded-off 
against  the  much  improved  properties  of  the  composite 
material  and  selective  stiffening  can  be  utilized,  with 
the  net  result  of  a siibstantially  stiffer  housing  with 
no  weight  penalty. 
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MAGNESIUM  CASE 


STEEL  CASE 


Figure  71.  Displacement  of  Internal  Components  Due  to  Ultimate  Load  Condition 
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TABLE  12.  DEFLECTION  AND  STRESS  SUMMARY 


LOAD  CONDITION 

SHAFT 

SLOPE 

MESH  DISPLACEMENT  - 
cm  (in. ) 

PINION 

GEAR 

PINION 

GEAR 

ULTIMATE 

i 

Magnesium 

.0017 

.0005 

.0429 

(.0169) 

.0091 

(.0036) 

Steel 

.0004 

.0001 

.0094 
( .0037) 

.0018 
( .0007) 

STEADY  FLIGHT  (Ig) 

Magnesiim 

.0006 

.0002 

1 

.0147 
( .0058) 

.0030 

(.0012) 

Steel 

.0001 

.0000 

.0033 

(.0013) 

.0005 
( .0002) 

LOAD  CONDITION 

TYPICAL  MAGNESIUM  HOUSING  STRESS  - 
kPa  (PSI) 

ULTIMATE 

+103425  to  +172375  (+15000  to  +25000) 

STEADY  FLIGHT  (Ig) 

+ 20685  (+3000) 

YAWING  MANEUVER 

+137900  (+20000) 

RECOVERY  FROM 

+34475  (+5000) 

ROLLING  PULLOUT 
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CH-47C  FORWARD  TRANSMISSION 
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Patton,  H.,  MODEL  YHC-IB  FORWARD  TRANSMISSION  DEFLECTION 
TEST,  Gleason  Works  Test  Report  1842,  1962. 


STRUCTURAL  OPTIMIZATION 


Structural  optimization  has  been  accomplished  utilizing  strain 
density  and  stress-ratio  resizing  techniques.  A strain  energy 
analysis  was  developed  by  Boeing  Vertol  under  ARO  sponsorship 
(Reference  6) . Expanding  upon  this  work,  a post-processor 
program  (S-83)  compatible  with  NASTRAN  was  developed  for 
analysis  of  the  strain  energy  density  distribution  throughout 
a structure  and  is  based  on  the  concept  that  for  a given  static 
load  a uniform  strain  density  throughout  yields  a maximum 
strength/minimum  weight  structure.  The  program  tabulates  each 
structural  element  from  highest  to  lowest  strain  density  and 
with  this  guidance  the  elements  of  high  strain  density  can  be 
altered  to  modify  the  structure.  Furthermore,  NASTRAN  Level 
16  includes  both  a strain  energy  and  a "fully  stressed  design" 
capability,  which  resizes  all  elements  in  a structure  for  a 
given  load  condition  using  a specified  allowable  stress  and  a 
stress-ratio  resizing  algorithm. 

A review  of  mathematical  optimization  methods,  such  as  described 
in  AGARD  LS-70,  was  conducted  and  the  applicability  of  these 
methods  to  this  program  was  assessed  (Appendix  A) . Although 
numerous  analyses  are  available,  few  are  presently  suitable  for 
application.  Hence,  the  strain  energy  and  fully  stressed 
design  methods  were  selected  as  the  best  available  methods  for 
use  herein. 

A strain  energy  density  distribution  analysis  for  structural 
optimization  was  conducted  using  S-83  for  the  transmission  case 
structure  to  determine  a maximum  strength/minimum  weight  design. 
The  ultimate  load  condition  was  used.  These  results  may  be 
used  for  either  adding  material  or  removing  material.  The 
area  for  adding  material  for  a more  uniform  strain  density  was 
found  to  be  in  the  vicinity  of  the  junction  of  the  input 
cartridge  and  case. 

A stiffness/weight  optimization  for  the  upper  cover  was  also 
conducted.  Nonlinear  and  azimuthally  varying  roller  and  ball 
bearing  loads  for  the  ultimate  flight  condition  were  used. 

The  upper  cover  is  a better  candidate  for  the  application  of 
optimization  techniques  because  of  its  relatively  high  weight 
(145  lb  as  compared  to  55  lb  for  the  case) , and  because  it 
transmits  the  hub  loads.  Figure  73  shows  the  model  as  well  as 
the  nonlinear  bearing  loads.  Selective  structural  element 
plots  generated  by  NASTRAN  for  30  structural  areas  of  high- 
est and  lowest  strain  density  are  indicated  in  Figure  74. 

NASTRAN  distorts  some  of  the  structural  elements  in  ortho- 
graphic plotting;  however,  these  plots  are  still  of  use  in 
visualizing  the  general  areas  of  the  higher  and  lower  strain 
densities.  Two  criteria  were  used  as  indicators  for  stiffness 


Figure  73.  CH-47C  Forward  Transmission  Bearing  Load  Distribution. 

Structural  Elements  With  Highest  Strain  Density  indicated. 
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Figure  74a.  NASTRAN  Orthographic  Plot  of  Areas  of 
Highest  Strain  Density  for  Upper  Cover 


Figure  74b.  NASTRAN  Orthographic  Plot  of  Areas  of 
Lowest  Strain  Density  for  Upper  Cover 
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improvement:  1)  stress  reduction  and  2)  bearing  outer  race 

misalignment.  There  are  different  types  of  misalignment,  such 
as  out  of  roundness  of  the  outer  race,  outer  race  warped 
azimuthally,  the  line  of  centers  not  geometrically  centered 
for  pairs  of  gears,  or  outer  race  tilted.  The  latter  criterion 
was  chosen  because  of  its  ease  of  calculation  from  the  NASTRAN 
output. 

The  optimality  principle  of  maximum  stiffness  for  a minimal 
weight  states  that  for  two  similar  structures  under  the  same 
loading  with  the  same  weight,  the  structure  with  the  more 
uniform  strain  density  is  stiffen.  It  follows  that  when  all 
the  strain  densities  are  equal  the  structure  is  the  stiffest 
possible.  For  the  typical  case  of  multiple  loadings,  the 
optimality  principle  must  be  modified  such  that  the  largest 
strain  density  in  each  element  for  all  load  conditions  is  the 
same  throughout  the  structure.  The  largest  strain  density  is 
not  necessarily  caused  by  the  same  loading  condition  in  all 
the  elements  (Reference  23) . In  addition,  the  analytical 
process  is  iterative.  There  may  also  be  side  constraints 
such  as  displacements,  member  sizes  and  stress  limitations. 
Finally,  the  final  optimized  design  may  not  be  readily  produc- 
ible, or  it  may  be  too  expensive  to  manufacture. 

This  implies  that  some  compromise  is  needed  for  all  these 
conditions.  For  simplicity  it  is  assumed  that  one  iteration 
with  the  addition  of  material  in  the  areas  of  highest  strain 
density  using  S-83  with  NASTRAN  is  the  best  practical  start. 

The  sample  S-83  output  for  the  upper  cover  at  the  ultimate 
flight  condition  is  given  as  Appendix  F (Table  F-1) , which 
includes  the  principal  stresses  as  well  as  the  element  numbers 
and  strain  densities.  This  output  is  used  to  illustrate  both 
the  addition  and  the  removal  of  material  for  stiffness  optimi- 
zation on  a weight  basis.  For  2014-T6  aluminum,  the  ultimate 
stress  in  tension  and  compression  is  65,000  psi;  the  yield  stress 
in  tension  is  55,000  psi;  and  the  yield  stress  in  compression 
is  58,000  psi.  From  Table  F-1  all  the  stresses  are  within  these 
limits.  The  principal  stresses  are  proportional  to  the 
strain  density,  or  more  precisely  the  "square"  of  the  principal 
stresses  in  the  matrix  sense  ( ) is  proportional  to  the 

strain  density.  Hence  the  strain  density  method  is  somewhat 
similar  to  the  fully  stressed  design  method  (FSD) . 

The  results  of  adding  material  are  summarized  in  Table  F-2. 

The  weight  penalty  was  15.7  lb.  The  thickness  of  the 
30  most  highly  strained  elements  was  approximately  doubled. 

Most  stresses  were  reduced  and  the  misalignment  was  reduced 
in  general.  For  excunple,  the  maximum  stress  went  from 
39,781  psi  to  33,500  psi  and  the  fore-to-aft  roller  bearing 


23.  Venkayya,  Knot,  and  Reddy,  ENERGY  DISTRIBUTION  IN  AN 
OPTIMUM  STRUCTURAL  DESIGN,  AFFDL-TR-68-156 . 
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misalignment  went  from  .0031  in. /in.  to  .0027  in. /in.  For 
stiffness  both  added  and  removed  (approximately  half  removed 
from  the  least  strained  elements) , the  weight  penalty  was 
2.64  lb.  The  maximxim  stress  was  34,367  psi,  and  the  fore-to- 
aft  misalignment  was  .0027  in. /in.  as  before.  This  is  shown  in 
Table  F-3.  For  weight  off  only  (Table  F-4) , the  cover  was 
13  lb  lighter,  the  maximum  stress  went  up  to  40,293  psi,  and 
the  fore-to-aft  misalignment  was  .0033  in. /in.  It  appears 
that  adding  stiffness  is  the  best  choice. 

A fully  stressed  design  analysis  of  the  CH-47C  forward  trans- 
mission upper  cover  was  conducted  using  NASTRAN  Level  16. 

Two  iterations  were  done  and  a value  of  .001  was  selected 
for  e , where : 


a is  a calculated  stress,  and  is  a defined  stress  limit 

(tension,  compression  and  shear) . Another  input  parameter, 

Y , was  selected  as  unity,  where: 


^NEW 


P 

OLD 

a 


l^a  + (1 


a ) Y 


is  the  new 

NEW 


property,  is  the  old  property  and 

UjjJj 


a = MAX 


il. 


for  all  structural  elements. 


Y is  an  iteration  factor  which  limits  the  property  change  in 
a single  iteration,  and  if  less  than  unity,  improves  the 
stability  of  the  iterative  process. 


The  maximum  change  in  any  property  is  limited  by  and 

ISuN-  «»ere 

< -|SW  < *SlAX* 
i 

and  P^  is  the  initial  value  of  the  property.  If  and  K^in 

are  unspecified  no  limits  are  imposed.  Kj>j;^X  and  Kjymj  were 
not  specified  for  the  test  run  here;  hence,  some  unrealistic 
thicknesses  for  the  upper  cover  resulted.  Nevertheless,  it 
is  felt  that  this  new  feature  in  NASTRAN  is  useful  for  weight 
reduction.  In  the  test  case  for  example,  the  results  were; 

Original  Weight  = 145  lb 

Iteration  no.  1 = 26.5  Lb 

Iteration  no.  2 = 34.8  Lb 

A sample  of  the  output  thicknesses  for  each  iteration  is  given 
in  Figure  75.  Since  no  restrictions  were  imposed  on  the 
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Figure  75.  Sample  Output  for  Fully  Stressed  Design  of  CH-47C  Forward  Transmission  Upper  Cover 
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Figure  75.  Continued. 


deflections  in  the  analysis,  some  deflections  became  unreal- 
istically large  (e.g.,  .6  in.). 

The  new  strain  energy  capability  of  NASTRAN  Level  16  has  been 
applied  to  the  CH-47C  forward  transmission  upper  cover.  The 
output  results  are  given  in  Table  F-5.  The  strain  energy  out- 
put is  calculated  for  Rigid  Format  1 (Static  Analysis)  only  by 
placing  "ESE  (Print,  Punch)  = ALL"  in  the  case  control  of  the 
NASTRAN  deck.  This  not  only  prints  out  the  strain  energy  for 
each  element,  but  also  punches  it  on  cards.  This  deck  may  then 
be  used  in  a post-processor  to  sort  the  strain  energies,  to 
list  the  percentage  of  total  strain  energy,  and  to  calculate 
the  accumulated  percentage.  The  limitation  that  strain  energies 
can  be  calculated  only  in  Rigid  Format  1 can  be  circumvented 
so  that  they  may  also  be  calculated  in  Rigid  Format  3 (normal 
mode  analysis)  by  utilizing  a post-processor  to  reformat 
punched  "ASET"  displacements  for  a given  mode  shape  from  Rigid 
Format  3 into  "SPC"  input  for  a Rigid  Format  1 run. 
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APPLICATIONS 


COMPOSITE  MATERIALS 


Current  cast  light  alloy  (magnesium)  transmission  housing 
technology  does  not  provide  an  optimum  support  structure  for 
power  train  dynamic  components  under  operating  loads.  These 
structures,  limited  by  current  materials  and  processing 
techniques  which  do  not  permit  structural  design  optimization, 
exhibit  excessive  deflections  and  displacements  under  load. 

The  metal  industry,  especially  the  light  metals  industry,  has 
been  seeking  alloys  with  higher  strength-to-weight  and  modulus- 
to-weight  values.  New  alloys  and  new  processing  methods  have 
achieved  small  improvements,  but  the  gains  are  no  longer 
proportionate  to  the  effort  required  to  attain  them.  Since 
all  of  the  widely  used  structural  metals  reach  a limit  of 
specific  strength  at  about  1 million  inches  and  a limit  of 
specific  modulus  at  about  100  million  inches,  research  is  being 
devoted  to  ways  of  getting  around  these  specific  strength  and 
specific  modulus  barriers.  High-modulus  fiber- reinforced 
composite  materials  offer  promise  in  providing  the  solution. 
These  materials  provide  the  needed  capability  for  selectively 
stiffening  the  housing  structure  for  reduced  deflection  and  for 
detuning  to  reduce  vibration/noise  levels. 

Studies  conducted  by  Boeing  have  indicated  the  overall  desir- 
ability and  payoff  to  be  gained  from  using  these  high  specific 
strength,  specific  modulus  materials  for  helicopter  transmis- 
sion housings.  For  example,  increasing  the  housing  wall 
stiffness  reduces  the  resulting  static  and  dynamic  displace- 
ments for  a specified  load  condition.  This  is  evident  from 

F = Kx 

where  F = applied  static  or  dynamic  load 
K = stiffness 
X = displacement 

Plots  of  displacement  versus  load  for  various  materials 
(Figure  72)  indicate  that  the  magnitude  of  housing  displacements 
can  be  reduced  substantially  by  using  the  stiffer  metal  matrix 
materials.  Steel  is  also  shown  in  the  figure  as  a point  of 
re ference . 

A number  of  composite  systems  have  evolved,  such  as  boron/ 
epoxy,  graphite/epoxy,  boron/aluminum,  graphite/aluminum,  FP/ 
aluminum  (FP  is  a trade  name  of  Dupont's  polycrystalline 
AI2O3  fiber)  and  FP/magnesium.  Each  of  these  has  its  own 
peculiar  advantages  and  limitations.  Several  metal  matrix 
candidate  systems  exhibit  specific  strength  in  the  range  of 
2 to  3 million  inches  and  specific  moduli  of  400  to  500  million 
inches.  Metal  matrix  composites  offer  unique  combinations  of 
improved  performance  and  reliability,  in  comparison  to  organic 
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composite  systems,  due  to  improved  shear  strength,  compressive 
strength,  resistance  to  environmental  degradation,  and  improved 
design  flexibility.  Also,  the  metal  matrix  composites  offer 
better  high-temperature  capability  for  elevated  temperature 
applications  such  as  helicopter  transmission  housings.  With 
improved  design  concepts,  more  sophisticated  use  of  materials 
can  be  made,  and  greater  weight  savings  realized  by  selectively 
strengthening  critical  areas  of  the  casting.  Reinforcing 
these  highly-loaded  areas  by  imbedding  high  strength  filaments 
would  produce  a more  efficient  design  and  a more  serviceable 
part.  Weight  limitations  and  the  potential  of  significantly 
reduced  maintenance  dictates  that  the  use  of  advanced  metal 
matrix  composites  for  many  critical  structural  applications 
be  examined  carefully. 

The  work  reported  herein  indicates  the  use  of  finite  element 
methods  utilizing  NASTRAN  in  conjunction  with  other  programs 
for  the  analysis  of  composite  materials  to  define  optimum  ma- 
terial configurations  and  orientations.  Micromechanics,  lamina 
theory,  and  the  total  constitutive  relation  for  a laminated 
plate  provide  the  basis  for  analysis  and  design  of  composite 
structures.  The  theory  of  composites  also  includes; 

• Thermal  stress  calculation 

• Equivalent  coefficients  of  thermal  expansion 

• The  determination  of  the  strains  and  stresses  in  each 
layer  of  the  laminate 

• Transverse  shear  stress  analysis,  interlaminar  shear 
stress 

• Laminate  interaction  diagram  depictions  based  on  the 
maximum  strain  theory  of  failure 

• Optimization  of  layups 

These  analyses  are  available  in  the  form  of  operational 
computer  programs.  Two  of  these  programs,  which  may  be  used 
as  pre-  and  post-processors  to  NASTRAN,  are  sxammarized  below. 

The  "Point  Stress  Laminate  Analysis"  (S71)  computer  program 
(Reference  24)  defines  material  properties  of  anisotropic 
composite  materials  by  computing  equivalent  orthotropic 
material  properties  suitable  for  use  in  a NASTRAN  analysis. 

This  may  also  be  used  after  the  NASTRAN  analysis  as  a post- 
processor to  obtain  interlaminar  and  laminar  stresses. 

Figure  76  illustrates  the  lamina  or  layer  coordinate  system 
(1-2)  that  is  transformed  to  the  laminate  (X-Y)  axis  system. 


24.  Reed,  D.  L.,  POINT  STRESS  LAMINATE  ANALYSIS,  Document 

FZM-5494,  AFML,  Advanced  Composite  Division,  WPAFB,  Ohio, 
April  1970. 
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The  resultant  stresses  and  moments  of  the  laminate  are  also 
shown.  These  represent  a system  that  is  statistically 
equivalent  to  the  stress  system  that  is  acting  on  the  laminate. 
The  program  will  accommodate  up  to  400  layers,  and  a point 
stress  analysis  can  be  performed  and  thermal  loads  and  point 
stresses  may  be  calculated. 

A second  computer  program  for  the  analysis  and  optimization 
of  laminated  composites  is  "COOP"  (Reference  25) . Optimization 
is  achieved  by  minimizing  an  objective  function,  which 
consists  of  terms  involving  laminate  weight,  cost,  stiffness 
and  strength.  A materials  property  table  is  included  in  the 
program  with  an  option  to  input  any  material  and  its  properties. 
Margins  of  safety  are  calculated  using  four  different  failure 
criteria.  Buckling  loads  for  inplane  end  load  and  shear  can 
also  be  calculated.  The  program  prints  out  the  inplane  and 
bending  matrices  for  the  selected  lamina.  The  program  takes 
the  elastic  moduli  Ell,  E22/  G12,  ®1,  ®2  and  the  Poisson's 

ratios,  ^12»  *^21»  of  the  individual  layers  within  the  laminate 
and  computes  the  elastic  moduli  Ex»  Ey,  Gxy»  ®x»  ®y,'^xy  and 
^^yx  for  the  final  laminate.  The  1-2  axes  refer  to  the  layer 
coordinate  system  and  the  x-y  axes  to  the  laminate  system. 

The  x-y  axes  correspond  to  the  x-y  loading  directions.  The 
program  uses  the  classic  laminated  plate  theory  (References  26 
and  27)  in  the  calculation  of  the  elastic  constants,  stresses 
and  strains. 

An  analysis  of  the  CH-47C  forward  transmission  housing  using 
a finite  element  NASTRAN  model  to  determine  the  optimum 
locations  and  orientations  of  composite  materials  for  selective 
stiffening  of  the  lower  housing  structure  was  conducted. 

Static  loads  representative  of  operating  conditions  were 
applied  to  the  model.  Based  on  a combination  of  a strain 
energy  and  fully  stressed  design  analyses  the  goal  is  to 
achieve  maximvim  stiffness.  The  NASTRAN  preprocesser.  Point 
Stress  Laminate  Analysis,  was  used  to  study  the  optimum  fiber 
configuration  and  to  define  the  anisotropic  material  proper- 
ties. The  determination  of  the  areas  for  stiffening  is 
contingent  on  the  load  condition.  Equivalent  orthotropic 
material  properties  for  a quasi-isotropic  graphite  composite 
material  (0°/+45°/90°  lamina)  were  obtained  using  computer 
program  S-71.  A sample  output  for  this  run  is  attached  as 
Appendix  G.  The  selective  addition  of  this  quasi-isotropic 
composite  material  to  the  transmission  case  stiffening  was 


25.  Dobyns,  A.,  COMPUTER  PROGRAM  FOR  ANALYSIS  AND  OPTIMIZATION 
COMPOSITES,  Users  Document,  October  1976. 

26'.  Tsai,  S.  W.  , Adcuns,  D.  F.,  and  Doner,  D.  R.,  ANALYSIS  OF 
COMPOSITE  STRUCTURES,  NASA  Report  CR-620,  November  1966. 

27.  Ashton,  J.  E.,  Halpin,  J.  C.,  Petit,  P.  H.,  PRIMER  ON 
COMPOSITE  MATERIALS:  ANALYSIS,  Technomic  Publishing 
Company,  1969. 
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based  on  the  strain  density  distribution.  The  finite  element 
model  was  augmented  to  incorporate  the  orthotropic  material 
properties  and  a static  analysis  was  conducted.  Using  the  ul- 
timate load  condition.  Table  13  summarizes  the  comparison  that 
was  made. 


TABLE  13.  STRESS  COMPARISON  - ORIGINAL  AND 
COMPOSITE  AUGMENTED  HOUSING 

WITH  QUASI-ORTHOTROPIC 
COMPOSITE  IN  AREAS  OF 
HIGHEST  STRAIN  DENSITY 

ELEMENT  ORIGINAL  STRESSES  (E  = 10.5  X 10^)  - 

NUMBER Mg  (E  = 6.5  x 10^  PSI) HIGH- DENSITY  GRAPHITE 


max 

min 

max 

min 

2272 

448 

-957 

460 

-781 

2242 

728 

-1219 

812 

-1084 

2273 

-330 

-871 

-303 

-630 

149 

-3.6 

-720 

-34 

-735 

1065 

-78 

-963 

-36 

-943 

2062 

13702 

-10206 

6489 

-3670 

2092 

12921 

-9229 

6552 

-2667 

2097 

2283 

-4533 

1984 

-2562 

112 

997 

-1115 

1266 

-896 

2100 

-3630 

-30426 

-989 

-14833 

2063 

25040 

-7062 

11530 

-1181 

The  original  stresses  have  generally  been  cut  to  about  half 
their  original  values.  Another  observation  is  that  the  addition 
of  material  must  be  done  very  selectively  in  patches,  rather 
than  over  broad  areas,  in  order  to  achieve  maximum  benefit 
from  the  strain  energy  method. 

The  use  of  composite  materials  for  transmission  housings  will 
•allow  stiffer  yet  lighter  structures.  Further,  vulnerability 
is  reduced  due  to  the  better  ballistic  tolerance  of  composite 
materials  and  survivability  is  improved.  It  is  possible  to 
extend  the  operation  of  a marginally  lubricated  bearing  if 
thermal  gradients  from  inner  to  outer  races  can  be  reduced, 
thus  preventing  loss  of  internal  clearance.  The  low  thermal 
conductivity  of  the  composite  housing  will  impede  heat  flow 
from  the  outer  race  and  tend  to  equalize  the  outer-inner  race 
temperatures.  If  material  properties  can  also  be  adjusted  to 
approximate  the  coefficient  for  expansion  of  the  bearing  race, 
an  ideal  condition  for  maintaining  bearing  clearance  exists. 
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VULNERABILITY/SURVIVABILITY  OF  A HELICOPTER  TRANSMISSION 


The  minimiim  vulnerability /survivability  standard  for  a contem- 
porary transmission  requires  continued  safe  operation  for  at 
least  30  minutes  after  damage  from  any  single  hit  by  a 7.62mm 
projectile  at  a range  of  100  meters.  More  specifically,  they 
are  designed  to  operate  at  the  power  required  for  flight  at 
the  speed  for  maximum  range  at  sea  level  standard  conditions 
and  primary  mission  gross  weight  for  not  less  than  30  minutes 
after  depletion  of  the  main  lube  system  lubricant. 

The  work,  conducted  herein  will  prove  valuable  in  designing 
transmissions  to  meet  both  the  current  goals  and  future,  more 
stringent  goals.  The  vulnerability/survivability  benefits  to 
be  derived  are  concentrated  in  two  areas  - the  housing  and  the 
internal  gear/bearing  system.  The  analytical  methods  provided 
by  NASTRAN  provide  more  accurate  load  path  definition  and 
hence  allow  the  designer  to  build-in  improved  load  path  redun- 
dancy. Also,  the  housing  can  be  designed  for  the  efficient 
use  of  composite  materials  with  properties  that  eliminate  the 
brittle  fracture  characteristics  displayed  by  magnesium  when 
si;i)jected  to  hydraulic  shock  loads.  Through  NASTRAN,  a stiffer 
case  can  be  designed  which  will  improve  the  load  capacity  of 
the  gears  and  bearings  by  decreasing  misalignment  and,  there- 
fore, theoretically  allow  the  use  of  smaller  components.  This 
would  reduce  vulnerability  by  reducing  the  vulnerable  area  of 
the  transmission.  This  size  reduction  is  not  likely  to 
materialize  in  practice  for  sometime  since  the  effect  of  mis- 
alignment on  load  capacity  is  not  yet  defined  with  sufficient 
confidence  to  permit  trade-offs  between  size  and  misalignment. 
The  immediate  practical  benefits  will  thus  be  confined  to 
somewhat  higher  capacity  at  the  sizes  determined  by  current 
design  methods,  which  will  contribute  to  better  survivability. 

Further  significant  advances  in  survivability  result  from  an 
improved  understanding  using  the  NASTRAN  analyses  of  the 
thermal  conditions  existing  during  normal  and  emergency  loss- 
of-l\abricant  conditions.  In  addition  the  gears  must  be 
designed  with  sufficient  clearance  to  allow  thermal  expansion 
at  the  gear  tips  and  roots.  It  is  possible  to  extend  the 
operation  of  a marginally  lubricated  bearing  if  thermal  gra- 
dients from  inner  to  outer  races  can  be  reduced,  thus  preventing 
loss  of  internal  clearance.  The  low  thermal  conductivity  of  a 
composite  housing,  for  example,  impedes  heat  flow  from  the 
outer  race  and  tends  to  equalize  the  outer-inner  race  tempera- 
tures. By  using  a NASTRAN  model  to  approximate  the  coefficient 
of  expansion  of  the  bearing  race,  an  ideal  condition  for 
maintaining  bearing  clearance  exists. 
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Survivability /vulnerability  in  transmission  systems  is 
concerned  with  the  loss  of  lubrication,  radar  cross  section, 
and  fail-safety.  The  latter  two  areas  are  discussed  further 
in  subsequent  sections.  It  has  been  shown  that  the  heat 
rejection  rate  of  fins  can  be  calculated  by  the  NASTRAN 
thermal  analyzer  for  an  elimination  or  a reduction  of  external 
oil  lines  for  cooling.  Elimination  of  the  oil  lines  and 
placement  of  the  cooler  in  proximity  to  the  case  would  reduce 
the  vulnerable  area  of  the  CH-47  from  11.4  to  6.5  square  feet. 
Using  a system  that  could  by-pass  the  cooler  in  an  emergency 
would  further  reduce  the  vulnerable  area  to  4.1  square  feet. 

If  the  sump  is  armored  the  vulnerable  area  becomes  2.3  square 
feet,  and  a nonlubrication  capability  would  reduce  this  to  0.6 
square  feet.  Since  experience  has  shown  that  no  ballistic 
damage  has  caused  drive  system  severances,  the  potential  exists 
for  a highly  survivable  helicopter  with  the  elimination  of  the 
vulnerability  of  the  lubrication  system.  The  NASTRAN  tht'xmal 
analyzer  is  a valuable  analytical  tool  for  such  a design  problem, 
problem. 

It  is  clear  from  Table  14,  which  relates  the  CH-47  vulnerable 
area  to  combat  kills,  that  the  drive/lube  system  is  the  most 
vulnerable  area.  Analysis  indicates  that  the  larger  the  area, 
the  greater  the  potential  for  a "kill" . Complete  loss  of 
oil  from  a transmission  due  to  a .30-  or  .50-caliber  hit  in 
the  cooling  system  will  occur  in  approximately  20  to  30 
seconds.  Because  of  the  high  heat  rejection  rate  and  the 
large  loads  carried  by  the  transmission,  complete  loss  of  oil 
is  considered  to  be  an  A-kill,  while  total  cooling  loss  with 
no  oil  loss  is  considered  to  be  a B-kill.  Of  particular 
concern  is  the  ability  of  a helicopter  to  continue  its  flight  for 
a limited  period  after  loss  of  transmission  lubricant.  Boeing 
Vertol  has  acc\imulated  considerable  experience  with  non- 
lubricated  operation  of  components.  Much  of  this  has  been 
under  actual  conditions  of  oil  starvation,  while  other  experi- 
ence was  obtained  with  tests  under  simulated  oil  system 
failures.  These  tests  indicate  that  the  bearings  are  the 
first  items  to  show  distress  and  that  tapered  roller  bearings 
and  ball  thrust  bearings  in  particular  are  susceptible  to 
damage  at  high  speeds  and  high  loads.  Spiral  bevel  gears 
appear  to  be  the  next  most  vulnerable  item,  although  their 
damage  may  be  partially  ascribed  to  adjacent  bearing  failure. 

The  establishment  of  an  emergency  nonlubricated  capability 
of  30  minutes  or  more  will  reduce  the  K-,  A-,  and  B-kill  heli- 
copter transmission  categories  to  almost  zero  and  thereby  save 
personnel  and  aircraft. 
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TOTAL  28.84  100.0  100 
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RADAR  CROSS-SECTION 


The  computer-generated  plotting  capability  afforded  by  NASTRAN 
is  useful  for  assessing  the  cross-sectional  area  projected  by 
a structure  when  observed  from  any  viewpoint.  Typical  radar 
cross-sectional  areas,  assuming  a transparent  fuselage,  are 
plotted  in  Figure  77  for  the  CH-47C  forward  transmission  at 
various  angles.  Any  other  orientation  may  also  be  evaluated 
simply  by  specifying  the  angular  orientation  about  three  axes. 
Utilizing  the  plotting  scale  factor,  the  cross-sectional  areas 
presented  can  be  calculated  from  the  plots. 

The  reinforcement  and  modification  to  the  wall  thickness 
considered  herein  have  no  effect  on  the  radar  cross-section  of 
the  housing  since  the  basic  size  and  geometry  are  unchanged. 
However,  the  modeling  capability  would  be  valuable  for 
determining  the  best  housing  geometry  for  minimum  overall 
cross-section  during  the  design  or  major  redesign  of  a trans- 
mission housing. 
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Figure  77.  Typical  Plots  of  CH-47C  Forward  Transmission  for  Evaluating  Radar 
Cross-Sectional  Area 
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Figure  77.  Continued 
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Figure  77.  Continued 


FAIL-SAFE/SAFE-LIFE  DESIGN 


Many  helicopter  fatigue  critical  components  are  designed  for 
safe  life,  whereby  they  are  assigned  a service  life  in 
operational  hours  and  removed  from  service  at  this  elapsed 
time  to  preclude  catastrophic  failure.  This  approach  is 
effective  in  preventing  failure  due  to  statistical  variations 
in  the  flight  spectrum,  loads  and  strength.  However,  most 
helicopter  component  failures  are  caused  by  the  unknown,  the 
unpredictable,  or  the  unexpected  (Figure  78) , such  as  material 
defects,  manufacturing  and  maintenance  errors,  pilot  errors 
or  battle  damage.  Current  philosophy  is  to  design  helicopters 
to  be  tolerant  of  such  defects  so  that  the  risks  from  unantic- 
ipated causes  may  be  minimized.  Defect  tolerance,  a design 
concept  whereby  an  incipient  or  partial  failure  will  not 
result  in  catastrophic  failure,  can  be  achieved  by  either 
providing  an  alternate  load  path  with  a short  but  sufficient 
life  (fail-safe)  or  permitting  the  planned  and  timely  detection 
of  fatigue  damage  or  operational  deterioration  (safe  crack 
growth) . 
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Figure  78.  The  Unexpected  Causes  Most  Helicopter  Component  Failures 


A fail-safe  structure  can  be  either  active  or  passive.  The  pas- 
sive system  uses  one  load  path  always  carrying  a load  and  one 
path  normally  not  loaded  but  becoming  loaded  if  the  first 
path  fails.  The  active  system  uses  two  paths  always  sharing 
the  load,  but  each  is  capable  of  taking  over  the  entire  load. 

The  results  of  the  NAS TRAN  model  analysis  can  be  used  to  study 
the  fail-safety  of  the  housing  by  evaluating  load  path  redundancy, 
that  is  the  ability  of  the  housing  to  continue  to  transmit  all 
loads  even  when  one  load  path  is  severed.  Load  paths  can  be 
artificially  severed  within  the  model,  and  the  stresses  in  the 
remaining  load  paths  recalculated  under  this  condition  using 
the  modified  model.  Methods  of  exact  analysis  of  modified 
structures,  based  on  matrix  manipulations  to  modify  the  initial 
analysis,  have  been  developed  in  an  effort  to  preclude  a 
complete  reanalysis  and  the  associated  cost  (Reference  28) . 

28.  Melosh,  R.  J.,  Johnson,  J.  R. , and  Luik,  R.,  STRUCTURAL 
SURVIVABILITY  ANALYSIS,  Philco  Ford  and  AFFDL,  Paper 
Based  on  Contract  AF33 (615) -5039 . 
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This  approach  is  an  extensive  subject  in  itself,  however,  and 
has  been  bypassed  herein.  Checking  of  stresses  and  deflections 
determines  if  adequate  gear  contact  and  bearing  alignment  is 
maintained.  Comparison  of  stresses  with  allowables  indicate 
the  degree  of  fail-safety. 

The  alternative  approach,  safe  crack  growth,  is  based  on 
fracture  mechanics.  It  may  utilize  a periodic  inspection  or 
an  integral  detection  system  with  an  obvious  failure  indicator. 
In  both  approaches,  parts  are  sized  so  that  the  rate  of 
fatigue  crack  growth  from  a defect  would  be  slow  enough  to 
assure  inspection  and  detection  before  any  component  failure. 
Defect  tolerance  is  achieved  through  a sub-threshold  crack 
growth  concept  for  an  assxraied  defect  size.  This  defect  may  be 
an  inclusion  inherent  in  the  material,  or  it  may  be  introduced 
during  manufacture  or  in-service.  Using  this  concept,  the 
component  is  sized  to  a stress  level  sufficiently  low  to  assure 
that  the  defect  will  not  propagate  during  the  life  of  the 
component.  No  inspections  other  than  those  specified  for 
normal  aircraft  maintenance  are  required. 

The  fracture  mechanics  analysis  uses  the  mean  of  scatter  crack 
growth  data,  a crack  simulation  model,  and  the  maximum  load  in 
the  mission  profile.  Crack  growth  rate  characteristics  of 
metals  are  measured  by  testing  standard  compact  specimens. 
Typical  test  results  are  shown  in  Figure  79,  where  data  is 
presented  as  crack  growth  rate  versus  the  stress  intensity 
factor  range  (AK)  where: 

AK  = Aa  M TTa  f (a) 

and 

Aa  = range  of  applied  fatigue  stress, 

(o  max  - o min) 

a = crack  size  parameter 

f (a)  = function  of  crack  size  and  specimen  or 
component  geometry 

A significant  characteristic  of  t-.ie  data  is  that  there  is  a 
point,  designated  the  AK  threshold,  below  which  fatigue  crack 
growth  is  extremely  slow  or  negligible.  Generally,  the  thres- 
hold stress  intensity  factor  range  will  differ  with  material. 
For  a given  material,  the  threshold  stress  intensity  is 
primarily  influenced  by  stress  ratio  and  temperature. 

Data  indicates  that  environment  and  loading  frequency  do  not 
significantly  influence  the  threshold  stress  intensity  range 
versus  stress  ratio  curve.  This  "mean-of-scatter " curve  is 
representative  of  medium  strength  steels  and  titaniums.  The 
ultimate  tensile  strengths  of  the  materials  tested  to  define 
the  curve  are  in  the  approximate  range  of  85  ksi  to  180  ksi. 
Substitution  of  the  threshold  stress  intensity  factor  range 
in  the  equation  for  AK  allows  determination  of  the  maximum 
crack  size  which  will  not  propagate  at  a given  stress  level. 
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This  maximuin  crack  size  will  vary  depending  on  the  crack  model 
(i.e.,  on  f (a))  selected  for  analysis.  Parametric  charts 
may  be  generated  for  selected  crack  models  which  will  define 
the  size  of  nonpropagating  cracks  at  various  stress  range 
conditions.  It  is  useful  to  convert  AO”  to  steady  and 
alternating  stresses,  and  this  type  of  presentation  for  the 
surface  flaw  defect  is  shown  in  Figure  80. 


FIIACTUIIi  MECHANICS  TYPE  TESTING  IDENTIFIES 


THE  THRESHOLD  STRESS  INTENSITY  LEVEL  IS  DEPENDENT  ON  STRESS  RATIO 


A MATERIAL  PROPERTY  CALLED  $ 

THRESHOLD  STRESS  INTENSITY  LEVEL  (AKy^) 

AT  WHICH  FATIGUE  CRACK  GROWTH  IS 

EXTREMELY  SLOW  OR  NEGLIGIILE  | 


CRACK  PROPAGATION 
RATE  aWaM- INCHES 
PER  CYCLE 


STRESS  INTENSITY  FACTOR  RANGE.  aK  - PSI  5/iC X 1,000 


STRESS  RATIO,  R 


Figure  79.  Crack  Growth  Data  for  the  Application  of  Fracture  Mechanics  to  Defect  Tolerant  Design 


USING  THE  EQUATIONS  OF  FRACTURE  MECHANICS,  THE  THRESHOLD 
STRESS  LEVELS  ASSOCIATED  WITH  VARIOUS  CRACK  GEOMETRY 
MODELS  CAN  BE  DISPLAYED  IN  FORMATS  AS  SHOWN 


0.0001  0.001  0.01  0.1  1.0 
THRESHOLD  DEFECT  DEPTH  b (INCHES) 


Figure  80.  Analytical  Model  for  Application  of  Fracture  Mechanics  to  Defect  Tolerant  Design 
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Realizing  that  additional  safety  normally  compromises  perform- 
ance, studies  were  made  to  assess  this  penalty.  One  example 
of  this  weight  penalty  due  to  the  safe  crack  growth  approach 
is  presented  for  a pitch  housing.  The  30  hour  post- indication 
crack  period  requires  an  11%  heavier  section  than  would  a safe 
life  version  of  the  same  part  as  shown  in  Figure  81. 


COMPONENT  WEIGHT 
LIFE  (3,600  HR) 
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Figure  81.  Weight  Study  for  Safe  Crack  Growth  - HLH  Pitch  Arm 


With  the  above  background,  it  is  apparent  that  application  of 
the  NASTRAN  finite  element  model  to  the  stress/fatigue  analysis 
of  transmission  housings  should  lead  to  better  understanding 
of  crack  propagation  and  the  safe  crack  growth  characteristics 
of  a structure.  Previous  experience  with  stress  determination 
at  a crack  head  using  finite  elements  has  indicated  a need  for 
an  extrapolation  post  processer  computer  program,  since  the 
stress  results  for  quadrilateral  or  triangular  structural 
elements  are  assumed  to  exist  at  their  C.G.  Using  this  extrap- 
olation scheme  has  led  to  good  correlation  with  analytical 
results  for  circular,  elliptical  and  V-notch  cracks.  The 
ability  to  include  temperature  in  the  analysis  is  also  useful. 

A simulated  .1  inch  longitudinal  crack  was  introduced  into  the 
NASTRAN  model  at  a highly  stressed  area  under  an  ultimate 
loading.  The  crack  could  have  resulted  from  manufacturing 
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imperfection,  fatigue,  or  missile  impact.  In  general,  local 
stresses  increased  reflecting  load  redistribution  by  the  crack, 
and  the  crack  surfaces  separated  radially  by  .032  inch.  A sum- 
mary comparison  of  stresses  before  and  after  introduction  of 
the  crack  is  given  in  Figure  82. 
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Figure  82.  Crack  Study 
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Previous  fail-safe/safe-life  considerations  were  primarily  ap- 
plied to  the  fuselage,  hub,  and  rotor  blades  (References  29  and 
30) . The  criteria  should  be  more  specifically  defined  with 
respect  to  transmissions  for  improved  future  designs.  As 
applied  to  helicopter  transmissions,  safe-life  criteria 
(fatigue  failure,  crack  propagation)  should  be  primarily 
applied  to  the  cover  and  retention  mountings.  Experience  has 
determined  that  the  retention  and  mounting  hardware  have  a 
relatively  high  incidence  of  fatigue  failure  for  the  CH-47C 
forward  transmission  (Reference  31) . However,  with  the  advent 
of  new  composite  materials,  safe-life  criteria  should  be 
extended  to  the  transmission  case.  Since  the  cover  transmits 
the  steady  and  vibratory  hub  loads  to  the  aircraft  airframe 
through  the  retention  lugs,  a NASTRAN  analysis  of  the  cover/ 
retention  lugs  yields  useful  safe-life  information.  In 
addition,  since  the  drive  system  is  one  of  the  heaviest  compo- 
nents of  a helicopter,  a post-processor  program  such  as  S-83 
to  calculate  the  strain  density  distribution  for  weight 
optimization  is  useful  for  weight  reduction.  The  strain 
density/stress  analysis  of  NASTRAN  would  also  indicate  the 
load  paths  for  fail-safe  considerations  and  would  improve 
crack-growth  characteristics  by  tending  to  equalize  the  stress 
distribution. 


29.  Feldt,  G.  V.,  and  Russell,  S.  W.,  FAIL-SAFE/SAFE-LIFE  IN- 
TERFACE CRITERIA,  Technology  Incorporated,  USAAMRDL  TR 
75-12,  Eustis  Directorate,  U.S.  Army  Air  Mobility  Research 
and  Development  Laboratory,  Fort  Eustis,  Virginia,  April 
1975,  AD  A009519. 

30.  Needham,  J.  F.,  FAIL-SAFE/SAFE-LIFE  INTERFACE  CRITERIA, 
Hughes  Helicopters,  USAAMRDL  TR  74-101,  Eustis  Directorate, 
U.S.  Army  Air  Mobility  Research  and  Development  Laboratory, 
Fort  Eustis,  Virginia,  January  1975,  AD006131. 

31.  USAAMRDL  TR  75-56B. 
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CRASHWORTHINESS 


Two  indices  which  have  been  utilized  to  assess  the  overall 
crashworthiness  of  an  aircraft  are  the  levels  of  acceleration 
experienced  by  the  occupants  and  the  preservation  of  the  integ- 
rity of  the  occupied  cabin  areas  during  a specified  crash  con- 
dition. The  former  index  reflects  such  characteristics  as 
fuselage  design  and  occupant  restraint  systems  and  is  not  with- 
in the  scope  of  the  work  reported  herein.  However,  the  latter 
index  is  directly  affected  by  the  design  of  a helicopter  drive 
system.  The  attachment  of  large  masses  such  as  engines,  trans- 
missions, and  rotors  to  the  upper  fuselage  aggravates  the  col- 
lapse of  the  structure  and  often  results  in  loss  of  occupiable 
volume  and  crushing  injuries  or  entrapment  of  occupants. 

Helicopter  transmissions,  due  to  their  location  typically  being 
above  or  adjacent  to  the  crew  or  passenger  cabin  areas,  pose 
a direct  hazard  to  occupants  during  a crash  situation  and  must 
be  designed  with  inherent  retention  strength  in  excess  of  the 
ultimate  crash  load  requirements.  Several  design  criteria 
relative  to  restraining  the  transmission  from  entering  occupied 
areas  must  be  satisfied  to  insure  crashworthiness:  retention  of- 
the  transmission  in  the  airfrcime,  sufficient  strength  of  the 
airframe  support  structure  to  preclude  total  collapse,  and 
integrity  of  the  transmission  structure  including  retention 
of  the  internal  components  in  the  housing.  The  first  two 
criteria  necessarily  address  topics  such  as  strengtii  of  attach- 
ment points,  strength  of  mounting  bolts  and  dissipation  of 
kinetic  energy  by  the  helicopter  structure  by  elastic  and 
plastic  energy  absorption.  The  strength  of  the  fuselage 
structure  can  be  accurately  assessed  using  a NASTRAN  model 
with  loads  representative  of  those  occurring  at  the  transmis- 
sion mounting  legs  placed  at  the  airframe  attachment  points. 

The  third  criterion,  structural  integrity  of  the  housing  itself, 
is  the  focus  of  the  work  discussed  herein.  The  NASTRAN  model  of 
the  transmission  housing  was  used.  Hub  loads  and  inertia  loads 
were  applied  to  the  transmission,  these  loads  were  converted 
to  bearing  reactions  on  the  housing  model,  and  a NASTRAN  static 
stress  analysis  was  performed.  The  strength  of  the  transmis- 
sion supporting  legs  and  the  structural  integrity  of  the 
housing  when  loads  representative  of  crash  conditions  are 
applied  were  evaluated. 

The  contractor's  current  experience  with  crashworthiness 
testing  has  indicated  that  loads  between  40  and  lOOg  should 
be  considered.  Assuming  a specified  initial  flight  condition 
for  a given  helicopter,  the  velocity  at  impact  can  be  calcu- 
lated. Then,  after  establishing  the  duration  of  the  crash  load 
impulse  based  on  recent  extensive  crashworthiness  testing, 
the  acceleration  and  hence  the  force  at  impact  may  be  determined. 
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utilizing  these  loads  in  the  NASTRAN  inertia  relief  analysis, 
the  housing  displacements  and  stresses  due  to  crash  loads  may 
be  calculated.  A comparison  of  stresses  with  allowables  provides 
an  evaluation  of  the  crashworthiness  of  the  housing  and  its 
support  structure.  The  strain  energy  distribution  after  impact 
is  a source  of  information  for  structural  improvement  for  crash- 
worthiness  . 

A crashworthiness  study  was  conducted  for  a CH-47C  forward 
pylon  using  the  inertia  relief  feature  of  NASTRAN  (Rigid  Format 
2) . The  weight  of  the  fuselage  section  was  5882  lb  and 
extended  rearward  to  Station  160.  The  forward  transmission 
was  attached  at  4 points  (cover  legs)  and  weighed  1081  lb. 

Above  this  was  a 683-lb  hub  assembly  with  three  blades  attached 
(1000  lbs) . The  angle  of  impact  was  10°,  and  the  impact  time 
was  AT  = .25  sec.  The  initial  altitude  was  55  feet,  giving 
a deceleration  of  7.4  g's.  These  numbers  reflect  a recent 
crash  test  at  Boeing  Vertol  conducted  under  Contract  DAAJ02- 
76-C-0015.  The  inertia  relief  capability  of  NASTRAN  is  linear 
elastic,  whereas  a true  crashworthy  analysis  is  elastoplastic . 
However,  NASTRAN  gives  a first-cut  solution.  The  model  was 
completely  free  and  "SUPORT"  cards  were  used  to  avoid  stiffness 
matrix  singularity.  Six  nonredundant  degrees  of  freedom 
were  constrained.  "C0NM2"  cards  were  used  for  the  mass 
distribution.  The  original  model  is  shown  in  Figure  83;  the 
model  after  impact  is  shown  in  Figure  84,  indicating  the  nose 
plowing  into  the  ground  but  with  no  resultant  transmission 
housing  impacting  into  the  crew  area.  Figure  85  shows  this 
earth  gouging  pictorially.  Figure  86  illustrates  the  need  for 
a crashworthiness  study  of  the  transmission  housing  due  to 
its  location  above  the  occupied  area.  The  analytical  results 
obtained  in  this  first- cut  solution  agreed  with  the  test 
results  under  Contract  DAAJ02-76-C-0015  (Crashworthiness 
Study)  in  that  the  housing  remained  intact  with  only  the  nose 
collapsing.  NASTRAN  (Rigid  Format  2)  was  difficult  to  inter- 
pret because  the  inertia  loads  of  the  concentrated  masses  of 
the  forward  pylon  were  not  printed  out.  In  addition,  "SPC" 
cards  could  have  been  used  instead  of  "SUPORT"  cards  since  the 
inertia-impact  loads  are  self-equilibrating . A second  inter- 
pretation of  crashworthiness  would  be  the  integrity  of  the 
transmissions  internal  components  after  impact. 
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Figure  85.  Earth  Gouging  (Plowing)  of  Fuselage  Under  Longitudinal  Impact 


Figure  86.  Schematic  Diagram  of  Idealized  Aircraft  With  Transmission  Concentrated  in 
Upper  Fuselage 
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EXPERIMENTAL  PROGRAM 


INTRODUCTION 


The  experimental  work  performed  under  this  contract  involved 
two  separate  test  programs.  The  uniform  temperature  testing 
of  a baseline  transmission  case  was  described  previously. 

The  testing  described  in  this  section  was  concerned  with  the 
measurement  of  housing  stresses  due  to  operating  loads. 

OBJECTIVE 


The  stresses  occurring  in  both  a baseline  and  a modified  ver- 
sion of  a CH-47C  forward  rotor  transmission  when  subjected  to 
steady-state  operating  loads  were  measured  to  experimentally 
determine  the  effect  of  selectively  stiffening  a transmission 
housing.  The  results  were  evaluated  to  determine  the  effec- 
tiveness of  the  analytically  predicted  structural  changes  for 
reducing  peak  stress  levels.  The  criterion  for  this  stiffening 
process  was  the  uniformity  of  the  strain  density  distribution. 

Test  Stand 


The  transmission  was  statically  and  dynamically  tested  in  the 
Boeing  Vertol  closed-loop  test  stand  (Figure  87) . This  stand 
employs  four  components  to  close  the  torque  loop.  First,  a 
set  of  helical  gears  increases  the  output  or  rotor  shaft  speed 
to  the  input  or  synchronization  shaft  speed.  A torque  device 
connects  this  gear  shaft  to  a bevel  gearbox.  The  bevel  gear- 
box closes  the  loop  to  the  input  shaft  of  the  transmission 
and  also  connects  to  a variable  speed  clutch  and  an  electric 
motor  which  drives  the  system.  This  closed-loop  test  stand 
provides  the  capability  of  running  a transmission  over  its 
full  design  torque  and  speed  range  under  controlled  conditions, 
including  rotor  lift,  drag,  and  pitching  moments. 

Data  Acquisition 

Strains  were  recorded  utilizing  twelve  strategically  located 
strain  gages  which  were  affixed  to  the  exterior  surface  of  the 
transmission  housing.  The  strain  gage  locations  are  presented 
in  Figure  88.  The  instrumentation  console  is  shown  in 
Figure  89. 

Test  Configuration 

The  transmission  used  in  this  program  was  a CH-47C  forward 
rotor  transmission,  except  as  structurally  modified  per  the 
strain  density  analysis  described  below.  The  baseline  trans- 
mission was  installed  in  the  closed-loop  test  stand  and 
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Figure  89.  Test  Stand  Instrumentation  Console 


load  representative  of  steady-state  flight  loads  (lift,  drag, 
pitching  and  torque)  were  applied.  Strain  gage  data  was 
recorded  for  both  static  (nonoperating)  and  dynamic 
conditions  with  the  transmission  operating  at  a synch  shaft 
speed  of  7460  rpm.  This  procedure  was  then  repeated  except 
that  the  transmission  case  was  selectively  stiffened  by  the 
addition  of  contour  doubler  plates  to  allow  for  a more  uniform 
strain  density. 

Areas  for  improving  the  existing  housing  were  determined 
analytically  by  the  finite  element  and  strain  energy  methods 
developed.  The  optimization  of  a housing  structure  should 
consider  both  the  addition  and  the  removal  of  material.  Never- 
theless, no  material  was  removed  from  the  housing  for  the  pur- 
pose of  the  testing  presented  herein,  since  at  the  time  of  the 
design  process  there  was  no  logical  criteria  for  selecting  the 
areas  for  weight  reduction  in  an  existing  transmission.  A fully 
stressed  design  (FSD)  capability  now  exists  in  NASTRAN  Level  16. 
An  FSD  analysis,  conducted  for  the  housing  subsequent  to  the 
testing,  indicated  several  areas  for  reducing  wall  thickness 
and  hence  weight  but  no  additional  hardware  was  built.  The 
casting  process  dictates  the  minimum  practical  thickness 
allowed.  Also,  a structure  is  generally  designed  using  limit 
loads  and  requiring  that  there  is  no  permanent  deformation. 

These  two  criteria  set  a limit  on  the  minimum  weight.  The 
latter  stress  analysis  would  determine  a strain  density 
distribution,  which,  when  made  uniform  would  yield  a stiffer 
cover/case . 

The  criterion  for  a maximum  stiffness  structure  for  a given 
weight  is  the  uniformity  of  the  strain  density  distribution 
(Reference  32) . In  practice,  this  strain  density  distribution, 
which  is  a tabular  computer  listing  of  the  structural  element's 
strain  density  from  highest  to  lowest,  indicates  the  most 
efficient  areas  to  add  material  in  order  to  obtain  a more  uni- 
form distribution  for  greater  strength.  The  areas  of  highest 
and  lowest  strain  density  were  determined  using  the  S-83  com- 
puter program.  The  areas  in  highest  strain  were  determined 
for  all  the  flight  conditions  calculated  (Ig  high-speed  level 
flight,  ultimate,  yawing,  recovery  from  rolling  pullout) , and 
the  areas  of  the  housing  which  were  modified  are  shown  in 
Figure  90.  Furthermore,  the  analysis  confirmed  that  both  the 
undesirable  structural  and  vibration/noise  characteristics  were 
closely  associated  with  similar  areas  of  high  strain  energy. 
Hence,  the  structural  changes  recommended  herein  are  compatible 
with  those  recommended  for  vibration/noise  reduction  studies 
in  Reference  9.  The  modified  transmission  housing  hardware  is 
shown  in  Figure  91. 


32.  Howells,  R.  W. , Sciarra,  J.  J.,  and  Ng,  G.  S. , THERMAL 
AND  STRUCTURAL  ANALYSIS  OF  HELICOPTER  TRANSMISSION 
HOUSINGS  USING  NASTRAN,  NASA  TMX-3428,  October  1976. 
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Figure  90.  Test  Housing  Stiffened  by  Attachment  of  Doubler  Plates 


‘ 


Figure  91.  CH-47  Forward  Transmission  Case  With  Modifications  (Cross  Hatched  Areas) 
to  Wall  Thickness  for  Selective  Stiffening 
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Test  Results  and  Correlation 


The  static  loading  condition  simulated  rotor  shaft  loads  at 
150  knots  high-speed  level  flight  (Ig)  with  25,000  lb  lift, 

260  lb  drag,  84,000  in. -lb  pitch  (positive  noseup)  and 
-720,000  in- lb  of  rotor  shaft  torque.  The  measured  stress 
values  for  the  baseline  housing  at  this  static  load  condition 
are  summarized  in  Table  15  and  compared  with  the  values 
predicted  by  the  NASTRAN  model.  Except  for  element  number 
156  and  the  two  questionable  strain  gages  noted,  the  predicted 
and  measured  values  correlate  quite  well.  The  sense  of  the 
stress,  tension  or  compression,  is  indicated  correctly  and 
the  magnitudes  are  within  a reasonable  experimental  scatter. 

The  values  predicted  by  NASTRAN  for  the  ultimate  load  condition 
are  also  shown  in  this  table. 


A comparison  of  the  magnitudes  of  the  measured  static  stresses 
for  the  baseline  housing  and  the  stiffened  housing  modified 
according  to  the  strain  energy  method  is  shown  in  Table  16.  Of 
the  ten  gages  that  functioned  during  the  test,  six  showed 
substantial  reduction,  three  showed  essentially  no  change,  and 
only  one  increased.  The  net  result  was  a significant  overall 
decrease  in  the  stress  levels.  Since  the  strain  energy  method 
utilizes  the  matrix  equivalent  of  the  square  of  the  strains 
(e.g.,  analogous  to  1/2  Kx^  for  a simple  elastic  system),  it  is 
appropriate  to  examine  the  stress  magnitudes. 


The  baseline  and  stiffened  housings  were  also  tested  dynamically 
at  an  input  pinion  speed  of  7460  rpm  and  with  the  same  Ig  load 
condition  applied.  The  dynamic  strains/stresses  for  all  ten 
strain  gages  were  reduced  by  the  addition  of  the  stiffener 
plates,  which  were  located  by  strain  density  principles.  Readings 
were  made  by  taking  peak-to-peak  readings  from  an  oscilloscope. 

The  stress  magnitudes  for  the  baseline  and  modified  housings, 
which  are  compared  in  Table  17,  indicate  an  overall  reduction 
of  stress  levels. 

In  addition  to  the  normal  experimental  error  due  to  equipment 
tolerances  and  data  recording,  the  placement  of  the  strain 
gages  was  quite  critical.  The  finite  element  model  predicts 
average  stresses  at  the  centroid  of  each  element.  Since  these 
stresses  can  vary  substantially  across  the  element,  the  strain 
gages  were  placed  as  close  to  the  center  of  each  element  as 
possible.  Furthermore,  the  magnitudes  of  the  strains  measured 
were  quite  small  and  hence  relatively  sensitive  to  small 
inaccuracies. 
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TABLE  15 


STATIC  STRESS  SUMMARY  - MEASURED  AND 
PREDICTED  VALUES  FOR  BASELINE  HOUSING 


LOCATION 
(ELEMENT  NO. ) 

BASELINE 

HOUSING  STRESS  (PSI) 

NASTRAN 
PREDICTED  (Ig) 

MEASURED  (Ig) 

NASTRAN 

PREDICTED 

(ULTIMATE) 

141 

-157 

-214 

-449 

149 

-228 

» 

79 

Gage  Erratic 

-667 

153 

-455 

-150 

-1297 

156 

-124 

780 

-317 

165 

807 

598 

2319 

1065 

127 

Gage  Lost 

371 

2062 

2095 

624 

6123 

2092 

1892 

-78 

Gage 

Questionable 

5442 

2097 

1485 

1404 

4265 

2242 

146 

150 

544 

2272 

41 

195 

675 

2273 

568 

598 

1641 
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TABLE  16 


STATIC  STRESS  StJdMARY  (MEASURED  VALUES) 
FOR  BASELINE  AND  STIFFENED  HOUSING 


LOCATION 

Ig  MEASURED  STRESS 

MAGNITUDE  (PSI) 

(ELEMENT  NO. ) 

BASELINE  HOUSING 

STIFFENED  HOUSING 

CHANGE 

141 

214 

215 

1 

149 

79 

Gage  Lost 

153 

150 

72 

-78 

156 

780 

156 

-624 

165 

598 

514 

-84 

1065 

Gage  Lost 

Gage  Lost 

2062 

624 

150 

-474 

2092 

78 

78 

0 

2097 

1404 

1125 

-279 

2242 

150 

215 

65 

2272 

195 

195 

0 

2273 

598 

514 

-84 

Average  Change  of 

Stress  Magnitude  = 

-156 

TABLE  17.  DYNAMIC 

STRESS  SUMMARY 

LOCATION 

MEASURED  STRESS  MAGNITUDE  (PS I) 

(ELEMENT  NO.) 

BASELINE  HOUSING 

STIFFENED  HOUSING 

CHANGE 

141 

260 

208 

-52 

153 

293 

260 

-33 

156 

325 

260 

-65 

165 

195 

130 

-65 

2062 

260 

234 

-26 

2092 

260 

234 

-26 

2097 

325 

260 

-65 

2242 

260 

208 

-52 

2272 

325 

182 

-143 

2273 

260 

156 

-104 

Average  Change  of 

Stress  Magnitude  = 

-63 
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CONCLUSIONS 


1.  Based  on  the  Thermal/S tatic/Dynamic  Analyses  conducted, 
it  is  apparerx  that  NASTRAN  can  be  used  as  an  effective 
tool  for  transmission  analysis  and  design.  In  fact,  there 
presently  exists  no  other  comprehensive  analytical  tool. 

The  heat  transfer/thermal  stress  capability  of  NASTRAN 

is  applicable  to  lubrication/ cooling  analysis.  The  stress/ 
dynamic  capability  is  applicable  to  analyzing  existing 
configurations  and  in  optimizing  new  configurations. 

2.  It  is  possible  to  construct  an  accurate  finite  element 
model  of  a transmission  utilizing  engineering  drawings  as 
evidenced  by  the  plotting  and  weight  generator  capabilities 
of  NASTRAN.  Essentially  the  same  model  can  be  used  to 
evaluate  static,  dynamic,  and  thermal  load  conditions. 

3.  The  ability  of  NASTRAN  to  accurately  predict  thermal 
distortions  of  a transmission  housing  has  been  verified 
by  correlation  with  test  data. 

4.  When  analyzing  the  housing  structure,  the  effect  of  the 
internal  components  (i.e.  gears,  bearings,  shafts)  must 
be  considered.  It  may  be  necessary  to  model  these  compo- 
nents either  in  detail  or  in  a simplified  manner.  Further- 
more, in  some  instances  such  as  the  thermal  growth  analysis 
it  may  be  possible  to  ignore  the  internal  components. 

5.  By  evaluating  the  displacements  of  the  housing  model  (grid  v, 

points)  at  the  bearing/housing  interfaces,  the  shaft  slopes 

and  displacements  at  the  gear  mesh  can  be  determined. 

Although  the  magnitude  of  these  displacements  can  exceed 
.010,  further  evaluation  is  needed  to  establish  the  precise 
effect  on  life  and  performance. 

6.  There  are  numerous  theoretical  techniques  for  structural 
optimization,  but  practical  methods  suitable  for  applica- 
tion to  the  design  process  are  minimal.  The  strain 
energy  method  and  the  Fully  Stressed  Design  (FSD)  method 
appear  to  be  the  only  methods  available  which  are  readily 
useable  by  a designer. 

7.  The  NASTRAN  model  predicts  stresses  which  correlate  with 
the  experimentally  measured  stress  values. 
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8.  The  structural  modifications  determined  from  the  strain 
energy  method  were  effective  in  reducing  the  overall 
stress  levels.  This  was  confirmed  from  measured  data  and 
analysis. 

9.  NASTRAN,  used  in  conjunction  with  other  pre-  and  post- 
processor programs,  has  shown  that  composite  materials 
can  be  effectively  employed,  at  least  in  theory,  to 
stiffen  transmission  housings.  Also,  the  thermal 
characteristics  of  a composite  housing,  which  have 
been  the  topic  of  concern  for  sometime,  can  be 
defined. 

10.  Numerous  aspects  of  the  vulnerability /survivability 

problem  can  be  investigated  using  finite  element  modeling. 
This  has  been  demonstrated  for  such  topics  as  radar  cross- 
section,  loss-of-lubricant  emergency  operation,  fail- 
safety,  crack  growth,  and  crash  worthiness. 
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RECOMMENDATIONS 


Based  on  the  results  of  this  effort,  it  is  recommended  that; 

1.  Further  study  be  conducted  to  establish  the  precise  effect 
of  displacement  and  misalignment  of  gears  on  the  life  and 
tierformance  of  a helicopter  transmission. 

2.  A further  and  quite  extensive  program  be  conducted  to  re- 
view existing  optimization  methods,  to  select  the  best  or 
develop  a new  method,  and  to  formulate  a practical 
computer-aided  design  procedure. 

3 . Research  be  concentrated  to  perfect  the  manufacturing  tech- 
nology necessary  to  utilize  composite  materials,  par- 
ticularly metal  matrix  materials,  for  the  selective 
stiffening  of  helicopter  transmission  housings. 

4.  Analyses  for  the  prediction  of  heat  generation  at  gear 
meshes  and  rolling  element  bearings  be  refined  and  extended. 

5.  Designers  be  encouraged  to  apply  the  finite  element 
methods  indicated  herein  to  future  helicopter  transmissions. 
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APPENDIX  A 


REVIEW  OF  OPTIMIZATION  TECHNIQUES 


The  goal  of  the  helicopter  designer  has  always  been  a structure 
that  meets  all  operational  requirements  and  also  is  of  minimum 
weight.  Loads  are  first  determined  for  the  critical  conditions 
(e.g.,  maneuver,  gust,  landing,  etc.)  and  the  structure  is  sized 
for  strength  based  upon  these  loads.  Modification  of  this 
original  "first  cut"  structure  will  generally  be  required,  and 
thus  an  iterative  procedure  is  involved.  A key  step  in  this 
procedure  is  the  optimization  of  the  major  structural  components 
for  strength,  based  upon  the  loads  applied  during  a given 
cycle . 

The  classical  approach  to  structural  optimization  has  been 
engineering  judgement.  Since  the  development  of  finite  element 
techniques  has  made  possible  the  analysis  of  complex  structures, 
the  enormous  number  of  structural  elements  and  side  constraints 
make  the  number  of  possibilities  for  structural  changes  large. 
The  design  trend  now  evolving  is  for  automated  sizing  methods 
which  will  optimize  the  complex  structures  for  some  specified 
constraint  criteria  (e.g.,  weight,  stiffness,  displacements,  and 
natural  frequency) . Typically  the  general  layout,  structural 
component  construction,  and  materials  are  assumed  to  be 
already  selected;  the  optimum  member  sizes  are  to  be  determined. 
Other  factors  such  as  weight  penalty,  location,  and  ease  of 
manufacture  must  also  be  considered. 

There  are  many  schemes  which  have  been  put  forth  for  both 
static  and  dynamic  optimization.  These  encompass  both  direct 
and  indirect  (optimality  criteria)  methods.  AGARD  LS-70 
(Reference  33)  reviews  several  of  these  methods,  including 

1.  Linear  programming 

2.  Nonlinear  programming  (direct  search) 

3 . Geometric  programming 

4.  Sieve-search 

5.  Inscribed  hyperspheres 

6.  Fully  stressed  design 

7 . Energy  methods 


33.  AGARD  Lecture  Series  No.  70  on  Structural  Optimization, 
AGARD-LS-70,  October  1974. 
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In  addition,  several  other  publications  covering  structural 
optimization  (see  the  List  of  References)  were  reviewed  with 
the  objective  of  assessing  the  applicability  of  these  methods 
and  the  feasibility  of  integrating  them  into  the  analytical 
technique  for  this  program.  Generally,  the  theoretical  de- 
velopment of  these  optimization  techniques  is  much  further 
advanced  than  their  practical  application  to  the  design 
of  large,  complex  structures.  Sample  applications  of  these 
methods  are  typically  limited  to  simple  structures  such  as 
trusses  since  the  methods  become  unwieldy  for  structures 
representative  of  those  found  in  practice.  As  a consequence 
of  the  abundance  of  analytical  techniques  and  lack  of  a solid 
practical  design  method,  the  need  for  a separate  and  quite 
extensive  program  is  indicated  which  would  have  the  objectives 
of  reviewing  existing  optimization  methods,  selecting  the  best 
or  developing  a new  method,  and  formulating  a practical 
computer-aided  design  procedure  based  on  the  method  selected. 

The  purpose  herein  is  not  to  develop  new  methods  of  structural 
optimization  but  rather  to  select  a method  from  those  currently 
available  and  apply  it.  No  conclusive  evidence  is  available 
to  show  that  one  method  of  optimization  is  clearly  superior. 
However,  it  is  currently  believed  that  while  mathematical 
programming  methods  are  suitable  for  detailed  component  opti- 
mization they  are  not  practical  for  the  overall  optimization 
of  large,  complex  structures,  where  optimality  methods  can  be 
applied  more  efficiently.  Two  of  the  many  existing  optimality 
methods,  strain  energy  density  (SED)  and  fully  stressed  design 
(FSD) , appear  to  be  at  least  as  good  as  the  others  and  have 
been  applied  in  practice  to  the  design  of  some  relatively 
large  structures.  These  two  methods  have  been  selected  for 
consideration  in  this  program. 

Fully  Stressed  Design  (FSD) 

NOTE;  Selected  portions  of  this  discussion  have  been  taken 
from  NASTRAN  documentation. 

The  Fully  Stressed  Design  (FSD)  technique,  using  the  stress- 
ratio  algorithm  to  drive  the  design  toward  a fully  stressed 
state,  is  probably  the  most  popular  resizing  method  for 
strength  optimization  (Reference  34) . NASTRAN  Level  16 
includes  a method  of  design  optimization  for  linear  static 
analysis  (Rigid  Format  1)  based  on  this  relatively  simple 
fully  stressed  design  strategy.  The  traditional  FSD  iterative 
procedure  is  based  upon  the  intuitive  belief  that  a given 
structural  configuration,  subject  to  stress  constraints  only, 
is  of  minimum  weight  when  the  stresses  in  all  the  elements 
are  at  prescribed  limits  under  at  least  one  design  loading 


34.  Dwyer,  W. , Rosenbaum,  J. , Shulman,  M.  and  Pardo,  H., 

FULLY-STRESSED  DESIGN  OF  AIRFRAME  REDUNDANT  STRUCTURES, 
AFFDL-TR  68-150,  October  1968. 
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condition  (Reference  35) . According  to  this  concept,  the 
cross-sectional  properties  of  each  structural  element  are 
resized  independently  at  each  design  iteration  to  produce  a 
limit  stress  (zero  margin  of  safety)  somewhere  within  the 
element  — assuming  that  in  each  iterative  cycle  the  internal 
load  distribution  is  unaffected  by  the  changes  in  its  cross- 
sectional  properties.  This  assumption  is  strictly  true  only 
for  statically  determinate  structures.  This  was  demonstrated 
by  Schmit  (Reference  36)  who  found  some  optimum  designs  which 
were  also  fully  stressed.  Most  practical  structures,  including 
a transmission  housing /ring/cover/internal  components  assembly, 
are  redundant. 

For  indeterminate  structures  of  low  redundancy,  the  above 
assumption  is  not  badly  in  error,  and  a few  repetitions  of  the 
algorithm  will  produce  a stress  distribution  throughout  the 
structure  which  has  very  nearly  a zero  margin  of  safety  in 
every  element  (i.e.,  a "fully  stressed  design").  The  procedure 
will  converge  more  slowly  (if  at  all)  in  structures  of  high 
redundancy,  and  modifications  of  the  basic  strategy  may  be 
required  to  achieve  convergence.  Furthermore,  there  is  no 
assurance  that  the  fully  stressed  design  of  a highly  redundant 
structure  will  be  an  optimum  design  in  any  meaningful  sense. 

It  is  relatively  easy  to  construct  examples  in  which  the 
procedure  converges  to  a "pessimum"  design.  Consider,  for 
excunple,  the  simple  case  of  two  parallel  rods  that  are 
rigidly  connected  together  at  their  ends  and  which  differ  only 
in  their  allowable  stresses.  Since  in  this  case  the  stresses 
in  the  two  rods  are  equal  regardless  of  their  areas,  the 
algorithm  will  increase  the  area  of  the  weaker  rod  at  the 
expense  of  the  stronger,  and  in  the  limit  only  the  weaker  rod 
will  remain. 

Practical  modifications  such  as  minimum  element  size  require- 
ments are  introduced  into  the  basic  iterative  resizing 
procedure.  Additional  complexities  associated  v/ith  members 
carrying  combined  bending  and  membrane  loads  and  biaxial  stress 
states  require  special  attention.  When  there  are  no  minimum 
element  size  requirements  to  be  satisfied,  the  nonuniqueness 
of  fully  stressed  designs  is  most  evident  since  there  are  at 
least  as  many  FSD's  as  there  are  combinations  of  statically 
determinate  member  formations  capable  of  supporting  the  spec- 
ified loads.  This  raises  some  doubt  as  to  the  convergence  of 


35.  Dwyer,  J.  W. , Emerton,  R.  K.  and  Ojalvo,  I.  V.,  AN  AUTO- 
MATED PROCEDURE  FOR  THE  OPTIMIZATION  OF  PRACTICAL  AERO- 
SPACE STRUCTURES  (Volume  I - Theoretical  Development  and 
User's  Information),  AFFDL-TR  70-118,  March  1971,  Air  Force 
Flight  Dynamics  Lab,  Wright-Patterson  AFB,  Ohio. 

36.  Schmit,  L.  A.,  STRUCTURAL  DESIGN  BY  SYSTEMATIC  SYNTHESIS, 
Proceedings  of  the  Second  Conference  on  Electronic 
Computation,  ASCE,  Volume  89,  August  1963. 

171 


related  synthesis  methods.  Although  attempts  have  been  made 
to  develop  general  convergence  criteria,  no  definitive  state- 
ments concerning  the  failure  of  convergence  to  a single  design 
by  FSD  based  iterative  methods  can  be  made. 

From  the  above  discussion  it  is  seen  that  a fully  stressed 
design  algorithm  cannot  be  used  uncritically.  Nevertheless, 

FSD  is  very  attractive  because  of  its  basic  simplicity  and  it 
will  produce  excellent  designs  in  many  practical  cases.  FSD 
based  procedures  have  delivered  practical  and  efficient  aircraft 
structural  designs.  Furthermore,  in  many  instances  where 
theoretical  optimums  were  computed,  it  was  found  that  they 
were  also  fully  stressed.  The  relative  simplicity  of  the 
application  of  FSD  techniques  also  accounts  for  their  further 
development  and  use. 

In  a sample  case  for  a three-spar,  five-bay,  swept  box  beam 
shown  in  Reference  34,  convergence  is  rapid  for  two  types  of 
FSD  design:  an  average  stress  method  and  a nodal  stress 
method.  Results  for  both  types  of  analysis  were  similar. 

The  weights  of  the  two  structures  change  very  little  after  the 
initial  resizing,  as  indicated  in  the  following  table. 


NODAL  STRESS  (lb)  AVERAGE  STRESS  (lb) 

Start 

179.5 

179.5 

1 Cycle 

119.9 

118.9 

3 Cycles 

119.5 

118.5 

20  Cycles 

119.2 

119 

Three  load  conditions 
was  required  that  the 

were  used.  For 
Henky-Von  Mises 

an  allowable  stress,  it 
effective  stress 

ae= 

-Jx 

3T  ) 

be  less  than  60,000  psi.  The  structure  and  loadings  are 
shown  in  Figures  A-1  and  A-2. 


For  each  of  the  two  fully  stressed  designs  the  skin  thicknesses 
were  initially  set  at  a uniform  0.16  in.  Twenty  redesign 
cycles  were  run  in  each  case,  and  the  resulting ‘cover  gages 
are  plotted  in  Figure  A-3. 

A fully  stressed  design  may  not  always  be  the  lightest  possible 
structure  satisfying  the  strength  and  minimum  gage  constraints. 
Work  is  being  done  in  this  area  using  true  optimization 
procedures  (i.e.,  those  which  focus  directly  upon  the  weight 
of  the  structure  and  attempt  to  minimize  it  subject  to  the 
aforementioned  constraints) . Currently,  the  more  promising 
approaches  include  the  penalty  fvinction  emd  the  gradient 
search  techniques.  To  date,  however,  those  fully  stressed 
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Figure  A-1.  Example  Structure  Geometry  Figure  A-2.  Example  Structure  Applied  Loads 
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Figure  A— 3.  Element  Thickness  (Inches)  as  a Function  of  the  Number  of  Cycles 
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designs  which  have  been  shown  ta  be  nonop tinnim  have  fallen 
into  a very  restricted  class,  and  even  for  this  group  they 
have  been  fairly  close  to  optimum.  In  view  of  the  much  greater 
computer  time  required  for  the  penalty  function  and  the 
gradient  search  techniques,  the  FSD  approach  is  considered  to 
be  an  acceptable  practical  compromise  for  the  present. 

The  physical  quantities  involved  in  the  NASTRAN  FSD  design 
algorithm  are  properties,  stresses  and  stress  limits.  The 
properties  may  include  thicknesses,  cross-sectional  areas  or 
moments  of  inertia.  Most  NASTRAN  elements  have  several, 
independent  properties.  They  also  have  several  types  of 
stresses  and  several  places  where  stresses  can  be  evaluated. 

The  stress  limits  include  those  for  tension,  compression,  and 
shear.  As  an  example,  the  design  iteration  algorithm  of  the 
simple  case  of  an  element  with  one  property  is  presented: 


Let 


(1) 


where  (T  = stress,  OJ  = stress  limit,  and  where  the  search 
for  a maximum  value  is  extended  over  all  user-identified 
stress  components  and  locations  and  also  over  all  designated 
loading  cases.  The  new  property  for  the  element  is  evaluated 
from  the  old  property  by  the  formula 


where  V is  a parameter  selected  by  the  user.  For  Y = \ (the 
default  value).  Equation  2 becomes 


A = oC.  A , , 
new  ^ old 


(3) 


If  the  product  Q-h  were  invariant.  Equation  3 would  give 


O^new  tTold  ^old 

new  ^ 

so  that  the  value  of  CT  ^ew  would  just  be  equal  to  the  limit 
stress  in  this  special  case. 

For  Y*  = 0,  it  is  seen  that  A^ew  = ^old»  for  values  inter- 

mediate between  zero  and  one  the  property  is  changed  by  less 
than  a factor  ofpC.  Thus  ■jp  is  a parameter  which  moderates 
the  property  changes  at  each  iteration  and  may  be  used  tr 
improve  the  convergence  of  the  algorithm. 

For  NASTRAN  elements  with  more  than  one  cross-sectional 
property,  all  of  the  properties  are  changed  according  to  a 
fixed  rule.  For  example,  in  the  case  of  the  BAR,  the  moments 
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of  inertia  are  changed  in  direct  proportion  to  the  change  in  the 
area.  This  is  equivalent  to  the  assumption  that  each  BAR  has 
a thin-walled  cross-section  whose  thickness  is  being  changed 
uniformly.  The  procedures  for  elements  with  more  than  orte 
cross-sectional  property  cannot  be  used  for  the  detailed 
design  of  individual  elements.  The  incorporation  of  more 
elaborate  procedures  in  NASTRAN  has  been  foregone  due  to  the 
inherent  limitations  of  the  fully  stressed  design  algorithm. 
Indeed,  it  is  not  clear  that  any  fully  automated  general 
purpose  design  procedure  can  successfully  cope  with  the 
simultaneous  requirements  of  overall  and  detailed  design. 

Strain  Energy  Density 

Strain  energy  techniques  have  been  applied  mainly  to  the 
dynamic  optimization  of  structures  by  shifting  eigenvalues 
(natural  frequencies)  away  from  exciting  frequencies.  The 
objective  for  an  aircraft  in  particular  is  a maximum  eigenvalue 
shift  for  a minimum  weight  change.  The  mode  shape  is  used  to 
find  the  strain  energy  content  of  the  components  of  the 
structure.  It  can  be  shown  analytically  that  a complex 
structure  can  most  efficiently  be  designed  dynamically  by 
ensuring  that  the  modal  density  differential  is  uniform 
throughout  the  structure.  The  density  differential  of  a 
structural  element  is  the  difference  between  the  strain  density 
and  the  kinetic  energy  per  unit  volume  (kinetic  density) . 

In  most  cases,  the  strain  density  may  be  used  as  an  approxi- 
mation of  the  density  differential  since  the  kinetic  density 
is  relatively  small.  This  objective  may  also  be  stated  alter- 
natively as:  (1)  Find  the  least  weight  structure  with  the 
largest,  lowest  natural  frequency,  or  (2)  Find  the  least 
weight  structure  for  a specified  natural  frequency. 

A method  for  optimizing  a structure  for  a given  dynamic  loading 
has  been  described  by  Sciarra  (Reference  37) . This  method 
has  been  applied  to  a medium  size  helicopter  (Boeing  Vertol 
Model  347)  in  high-speed  level  flight.  The  results  are  shown 
in  Figure  A-4.  The  excitations  of  the  fuselage  are  the  n per 
rev  hub  loads  which  are  approximately  10%  of  the  gross  weight. 

It  is  seen  that  the  heights  of  the  sticks  that  are  representa- 
tive of  the  vibratory  response  levels  are  reduced  after  proper 
optimization.  Figure  A-5  shows  that  the  second  natural  frequen- 
cy (73.4)  has  been  moved  away  from  the  exciting  frequency  (79.5). 
These  results  are  for  the  forward  pylon  only.  This  type  of 
experience  is  directly  applicable  to  the  design  of  an  optimum 
transmission  housing. 


37.  Sciarra,  J.  J. , USE  OF  THE  FINITE  ELEMENT  DAMPED  FORCE 

RESPONSE  STRAIN  ENERGY  DISTRIBUTION  FOR  VIBRATION  REDUCTION, 
presented  at  ARO-D  Military  Theme  Review,  The  Helicopter 
and  V/STOL  Aircraft  Research  Conference,  U.S.  Army  Re- 
search Office,  NASA-Ames  Research  Center,  Moffett  Field, 
California,  AD-751809,  September  1972. 
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Figure  A —5.  Forward  Pylon  Spectrum  — Original  and  DFR  Modified 
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The  dependency  of  the  static  strain  energy  distribution  on  the 
loading  condition  makes  trade-offs  necessary.  However,  if  the 
first  eigenvalue  (natural  frequency)  is  maximized  for  a minimum 
weight,  this  would  make  the  optimization  independent  of  the 
loading.  This  can  be  seen  in  a violin  string  as  its  tension  is 
increased  yet  its  weight  is  constant.  Also,  different  eigenvalues 
(first,  second,  etc.)  give  different  optimal  designs.  As  an 
example,  optimal  designs  for  different  eigenvalues  of  the  portal 
frame  are  shown  in  Figures  A-6  and  A-7.  The  finite  element 
models,  the  first  four  mode  shapes,  and  the  optimal  designs  for 
the  first  four  modes  are  Shown. 

The  Strain  density  distribution  concept  is  also  utilized 
statically  to  identify  structural  load  paths  and  to  evaluate 
the  stress  efficiency  of  the  structural  design  (stiffness/ 
weight) . Venkayya  (Reference  38)  uses  the  finite  element 
method  in  an  iterative  cycle  to  find  a minimum  weight,  maximum 
strength  design  for  a given  static  load  condition  of  a structure. 
The  resizing  of  the  structural  elements  of  the  complex  model 
is  contingent  on  the  strain  density  distribution.  The  optima- 
lity principle  is  that  the  strain  density  should  be  uniform. 

A resizing  mechanism  establishes  this  uniformity.  Since 
competing  objectives  are  involved  in  the  optimization,  trade- 
offs are  necessary.  For  example,  different  static  load 
conditions  give  different  optimal  designs.  The  optimization 
analysis  of  a structure  is  also  nonlinear,  and  an  iterative 
loop  should  be  employed. 

The  beauty  of  the  strain  energy  method  is  that  strain  energy 
is  easily  obtained  using  finite  elements.  If  the  displacements 
and  rotations  of  a structural  element's  nodes  are  obtained  and 
called  {X},  a column  matrix,  then  1/2  {x}"^  [K]  {X}  is  the 
strain  energy  of  the  particular  structural  element  where  {K} 
is  the  stiffness  matrix. 


38.  Venkayya,  V.  B. , Khot,  N.S.,  and  Reddy,  V.  S. , OPTIMIZATION 

OF  STRUCTURES  BASED  ON  THE  STUDY  OF  STRAIN  ENERGY  DISTRIBUTION, 
AFFDL-TR  68-] 50,  1968,  Pages  111-153. 
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CONCLUSIONS  REGARDING  STRUCTURAL  OPTIMIZATION 

1.  There  are  numerous  theoretical  techniques  for  structural 
optimization,  but  practical  methods  for  application  to  the 
design  process  are  minimal. 

2.  Structural  optimization  for  this  contract  was  performed 
using  strain  energy.  FSD,  as  it  appears  in  NASTRAN  Level 
16  was  also  evaluated  herein.  The  FSD  method  appears  to 
be  a feasible  approach  to  structural  optimization. 

3.  A further  and  quite  extensive  program  is  necessary  to 
review  existing  optimization  methods,  to  select  the  best 
or  develop  a new  method,  and  to  formulate  a practical  com- 
puter-aided design  procedure. 
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APPENDIX  B 


SAMPLE  CALCULATIONS  FOR  HEAT  GENERATION 


Sample  calculations  for  CH~47  input  pinion  spiral  bevel 
angular  contact  ball  bearing  no.  15  (Reference  15) . 


= 400  i Z D cos 


2 f^  (1  - y) 


2 f.  - 1 
1 


1/2 


1 
Z 
f , 


Basic  Static  Capacity,  lb 
Number  of  Rows, 

Number  of  Rolling  Elements  Per  Row 
r^  , Inner  Race  Curvature  (From  S04) 

d“ 

D cos  a 
d 


M, 


m 


f . Fft  d 
1 S m 


= Z 


Bearing  Friction  Torque  (From  Palmgren) 
in. -lb 
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Fg  = Static  Equivalent  Load,  Lb  (From  S04) 
z,  y from  Table  B-1 

F^  “ “ Applied  Axial  Load,  Lbs  (From  S04) 

TABLE  B-1.  VALUES  OF  z AND  y 


Ball  Bearing  Type 

s 

y 

Deep  groove 

a = 0° 

0.0009 

0.55 

Angular  contact 

a = 30" 

0.001 

0.33 

Angular  contact 

o » 40" 

0.0013 

0.33 

Thrust 

o = 90" 

0.0012 

0.33 

Self-aligning 

a = 10" 

0.0003 

0.4 

For  Bearing  no.  15  at  100%  torque,  7460  rpm 

r 

^ / A ,T\  /T,\  /n  .%,..0  J 2 ( 


Cs  = 

(400) 

(1)  (14) 

= 

30287 

lb 

where  i = 

400 

Z = 

1 

D = 

1.125 

a = 

250 

f = 

52 

^m  ~ 

6.1 

f.  = 

.001 

r 133 

6111  * 

X 

30287 

M = 

(.00059)  (6050) 

fi-  (1-125?  (.! 
.52)1^  6.1 

2(.52)  -1 


= .00059 


= Mechanical  Friction  Torque 

M^  = Viscous  Friction  Torque  (in. -lb) 

= (1.42)  (10"^)  f (V^  N)^/^ 

o o m 

V^  = Centistokes  (e.g.  3 cs  at  200°F,  MIL-L-7808) 

N = RPM  (e.g.  , 7460) 

f = Factor  (Table  B-2) 
o 

M^  = (1.42)  (10"^)  (2)  (3  * 7460)^^^  d^^^  = 5.1  in. -lb 
'V)TAL  = M,  + M^  = 21.8  + 5.1  = 26.9  in. -lb 
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HP  = (26.9)  (7460)  = 3.18  HP 

63025 

Heat  Generated  by  Bearing  no.  3 = (42.42)  hp  = 134.9  Btu/min 


TABLE  B-2.  VALUES  OF 


Bearing  Type 

Mist 

Lubrication 

OUBath 

Lubrication* 

Grease 

Lubrication 

Venical  Mounting 
Flooded  Oil 
Lubrication 

Jet  Lubrication 

Deep-groove  ball 
bearings  (single  row), 
Self-aligning  ball 
bearings  (double  row). 
Ball  thrust  bearings 

0.7- -.It 

1.5- • -2 

3-  -4 

Filling-slot  ball 
bearings  (single  row). 

1 

2 

4 

bearings  (single  row) 

Angular-contact  ball 
bearings  (double  row) 

2 

4 

8 

Single-row,  tapered 
roller  bearings 

Spherical  roller  thrust 
bearings 

1.5--2 

3... 4 

6--8 

Cylindrical  roller 
bearings  (single  row) 

1 •••1.5 

2--3 

4--6 

Spherical  roller 

bearings  (double  row) 
Tapered  roller  bearings 
(double  row) 

2-  • -3 

4 - -6 

8-12- 

* Oil  level  reaches  center  of  lowest  rolling  element  in  a horizontal  mounting, 
t Lower  values  pertain  to  light  series;  higher  values  pertain  to  heavier  dimension 
scries. 
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Sample  calculation  for  CH-47  input  pinion  roller  bearing  heat 
generation  (Reference  15)  . 


Assume  for  roller  bearing  no.  1 the  Palmgren  relation 
(Bearing  Friction  Torque)  = f^  (bearing  factor)* 

Fg  (bearing  shaft  force)  * djj,  (pitch  diameter) 

or  = (.0003)  (15485  lb)  (6.5  in.)  = 30.4  in. -lb  The  fj^  chosen 

is  for  cylindrical  roller  bearings. 


Now  assume  for  lubricated  roller  bearing  viscous  friction 
torque  (Mv)  the  Palmgren  relation; 

= 1.42  X 10“5  fQ  (VqN) d^  where 

m 

fjj  is  a factor  depending  on  the  type  of  bearing  and  the  method 

of  lubrication  (e.g.,  2.5  for  single-row  cylindrical  roller 
bearings) , 

Vq  = centistokes  (e.g.,  3 cs  (§  200OF,  MIL-T-7808) 

N = rpm  (e.g.,  7067) 

or  My  = (1.42)  (10"5)  (2.5)  (3  x 7067) (6.5)^ 

My  = 7.47  in. -lb 

MipoTAli  ~ ^ ~ 37.87  in  • — lb 

(37.87)  (7067)=  a o K>  = 2TrNMTOTAL 

^ ~ 63025  T330MyTT2) 

Heat  Generated  = (42.42)*  = 178  Btu/min 

(1  HP=  42.42  Btu/min) 
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APPENDIX  C 

SAMPLE  OUTPUT  FOR  INDUCED  MISALIGNMENT  PROGRAM 
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2060  5.25000  60.00000  14.93400  2240  5.25000  240.00000  14.93400 
2090  5.25000  90.00000  14.93400  2270  5.25000  270.00000  14.93400 
2120  5.25000  120.00000  14.93400  2300  5.25000  300.00000  14.93400 
2150  5.25000  150.00000  14.93400  2330  5.25000  330.00000  14.93400 
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SAMPLE  OUTPUT 


V=  5.0  XL2=  12.500  T-  0.500 

0.80107E  00  1 0.43353E  00  o' 

0.76902E  00  1 0.41365E  00  0 

0.76902E  00  0 0.11534E  01  0 

0.76902E  00  1 0.91479E  00  0 

0.80107E  00  1 0.8e003E  00  0 

0.76902E  00  0 0.59285E  00  0 

BTU/HR=  65.83  TOP-  0.40840E-01 

V=  10.0  XL2=  12.500  T-  0.500 

0.10582E  01  1 0.60672E  00  0 

0.10159E  01  1 0.57636E  00  0 

0.10159E  01  0 0.12751E  01  0 

0.10159E  01  1 0.58592E  00  0 

0.10582E  01  1 0.54592E  00  0 

0.10159E  01  0 0.41158E  00  0 

BTU/HR=  114.85  TOP-  0.40844E-01 

V=  15.0  XL2=  12.500  T-  0.500 

0.12455E  01  1 0.75161E  00  0 

0.11957E  01  1 0.71127E  00  0 

0.11957E  01  0 0.13907E  01  0 

0.11957E  01  1 0.43751E  00  0 

0.12455E  01  1 0.40489E  00  0 

0.11957E  01  0 0.32037E  00  0 

BTU/HR=  159.08  TOP-  0.40849E-01 

V=  20.0  XL2-  12.500  T-  0.500 

0.13983E  01  1 0.88452E  00  0 

0.13424E  01  1 0.83415E  00  0 

0.13424E  01  0 0.15038E  01  0 

0.13424E  01  1 0.34951E  00  0 

0.13983E  01  1 0.32164E  00  0 

0.13424E  01  0 0.26296E  00  0 

BTU/HR-  200.44  TOP-  0.40853E-01 

V-  25.0  XL2-  12.500  T-  0.500 

0.15296E  01  1 0.10114E  01  0 

0.146848  01  1 0.95068E  00  0 

0.14684E  01  0 0.16162E  01  0 

0.14684E  01  1 0.29033E  00  0 

0.15296E  01  1 0.26588E  00  0 

0.14684E  01  0 0.22277E  00  0 

BTU/HR-  239.79  TOP-  0.40857E-01 

V-  30.0  XL2-  12.500  T-  0.500 

0.16460E  01  1 0.11350E  01  0 

0.1580 IE  01  1 0.10637E  01  0 

0.15801E  01  0 0.17286E  01  0 

0.15801E  01  1 0.24747E  00  0 

0.16460E  01  1 0.22S64E  00  0 

0.15801E  01  0 0.19280E  00  0 

BTU/HR-  277.62  TOP-  0.40862E-01 

V-  35.0  XL2-  12.500  T-  0.500 

0.17513E  01  1 0.12571E  01  0 

0.16812E  01  1 0.11747E  01  0 

0.16612E  01  0 0.18417E  01  0 

0.16812E  01  1 0.21487E  00  0 

0.17513E  01  1 0.19514E  00  0 

0.16812E  01  0 0.16948E  00  0 

BTU/HR-  314.21  TOP-  0.40866E-01 

V-  40.00  XL2-  12.500  T-  0.500 

0.18479E  01  1 0.13786E  01  0 

0.17740E  01  1 0.1284aE  01  0 


BESSEL  FUNCTION  INPUT/OUTPUT 

* OF  FINS 

BOT-  0.12490E  01  0.31102E  04 

BOT-  0.94584E  00  0.17827E  04 

BOT-  0.80387E  00  0.12871E  04 

BOT-  0.71629E  00  0.10205E  04 

BOT-  0.65501E  00  0.85386E  03 

BOT-  0.60888E  00  0.73752E  03 

BOT-  0.57244E  00  0.eS162E  03 


HEIGHT  OF  FINS  (FT) 
0.12959E  03 


0.74280B  02 


0.53630E  02 


0.42563E  02 


0.35577E  02 


0.30730E  02 


0.27151E  02 
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THEORETICAL  SOLUTION 

A metal  fin  of  triangular  cross  section  is  attached  to  a plane  surface  to  help  carry  offbeat  from  the  latter.  Assum- 
ing dimensions  and  coordinates  as  shown  in  the  accompanying  figure,  find  the  steady -state  temperature  distribu- 
tion along  the  fin  if  the  root  (i.e.,  wall)  temperature  is  and  if  the  fin  radiates  freely  into  air  of  constant 
temperature 

We  shall  base  our  analysis  upon  a unit  length  of  the  fin,  and  shall  assume  that  the  fm  is  so  thin  that  temperature 
variations  parallel  to  the  base  can  be  neglected.  We  also  assume  that  0 is  so  small  that  cos  6 may  be  replaced 
by  1. 

Now  consider  the  heat  balance  in  the  element  of  the  fin  between  x and  x + A*.  This  element  gains  heat  by  internal 
flow  through  its  right  face  and  loses  heat  by  internal  flow  through  its  left  face  and  by  radiation  through  its  upper 
and  lower 


surfaces.  Through  the  right  face  the  gain  of  heat  per  unit  time  is 
Area  X thermal  conductivity  X temperature  gradient 


(i  x^)  fk/  h. 

_ 

bkxdu 

\ a 1 \i 

ix) 

X + Ax 

a dx 

X V- Ax 


Through  the  left  face  the  element  loses  heat  at  the  rate 
\bkxdu\ 


b 

i 


a dix  J X 

. Arough  the  surfs  . exposed  to  the  air  the  element  loses  heat  at  the  rate 
Area  X outer  conductivity  X (surface  temperature  - air  temperature) 


(h)(u-  UgJ  = 2h(u  - uJAx 


= /2  X 1 ~^\ 

\ cos  6 1 

litions  the  rate  ol 

hfcydnl  _ I bkxdu 

a dx  x+Ax  a dx 


Under  steady-state  conditions  the  rate  of  gain  of  heat  must  equal  the  rate  of  loss,  and  thus  we  have 

r* 


L 


+ 2h  (u-  uJAx 
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Writing  this  as 


-x^ 

dx  x+Ax  dx  X _ 2ah  „ i = q 
Ax  bk  “ 

and  letting  Ax  O,  we  obtain  the  differential  equation 
dx  bk  “ 


If  we  set 


U = u -u„  and 

“ bk 


this  becomes 

dx 

The  general  solution  of  this,  according  to  the  theory  of  Sec.  8.3,  is 
U = u -Ug  = CiJo(2od  + CiYq  (Tai 
or,  using  the  modified  Bessel  functions, 

u - =C3/o(2a  + c^Koilas/x) 

Since  Ko(2oc  x)  is  infinite  when  x = 0/:^  must  be  zero,  leaving 
u - Ua  = C3/o(2a  y/xj 
When  X =a,  u = and  thus 

Therefore 

^ /„(2a  a) 
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NASTRAN  EXECUTIVE 


CONTROL 


DECK  ECHO 


ID 

HEAT,  FLOW 

APP 

HEAT 

SOL 

1.0 

TIME 

5 

CEND 

COOLING  FIN 

JANUARY 

21,  1976  NASTRAN  2/1/73 

CASE 

c 

0 N T R 0 L 

DECK  ECHO 

CARD 

COUNT 

1 

TITLE#COOLING  FIN 

2 

OUTPUT 

3 

THERMALtALL 

4 

ELFORCE#ALL 

5 

SPCFtALL 

6 

SUBCASE  1 

7 

LABEL #rEMPERATUBE  SPECIFIED 

8 

SPC#100 

9 

SUBCASE  2 

10 

LABEL#HEAT  INPUT  AT  ROOT 

11 

LOAD# 10 

12 

SPC#7 

13 

SUBCASE  3 

14 

LABEL*TEMFERATURE  SPECIFIED 

AND 

HEAT  INPUT  AT 

ROOT 

15 

LOAD*  10 

16 

SPC#100 

17 

BEGIN  BULK 

PAGE  1 


PAGE  2 
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CARD 


COUNT 

1 

2 

3 

4 

5 

. . 6 

. . 7 

. . 8 

. . 9 . . 1 

1- 

CHBDY 

101 

200 

AREA4 

1 

5 

7 

2 

&CH1 

2- 

&CH1 

101 

105 

107 

102 

3- 

CHBDY 

102 

200 

AREA4 

2 

7 

9 

3 

&CH2 

4- 

&CH2 

102 

107 

109 

103 

5- 

CHBDY 

103 

200 

AREA4 

1 

4 

6 

2 

&CH3 

6- 

&CH3 

101 

104 

106 

102 

7- 

CHBDY 

104 

200 

AREA4 

2 

6 

8 

3 

&CH4 

8- 

&CH4 

102 

106 

108 

103 

9- 

CHBDY 

105 

200 

AREA4 

5 

11 

13 

7 

&CH8 

10- 

&CH8 

105 

111 

113 

107 

11- 

CHBDY 

106 

200 

AREA4 

7 

13 

15 

9 

&CH7 

12- 

&CH7 

107 

113 

115 

109 

13- 

CHBDY 

107 

200 

AREA4 

4 

10 

12 

6 

&CH5 

14- 

&CH5 

104 

110 

112 

106 

15- 

CHBDY 

108 

200 

AREA4 

6 

12 

14 

8 

&CH6 

16- 

&CH6 

106 

112 

114 

108 

17- 

CHEXA2 

201 

1 

10 

12 

13 

11 

4 

6 &E201 

18- 

&E201 

7 

5 

19- 

CHEXA2 

202 

1 

12 

14 

15 

13 

6 

8 &E202 

20- 

&E202 

9 

7 

21- 

CWEDGE 

301 

1 

1 

5 

4 

2 

7 

6 

22- 

CWEDGE 

302 

1 

3 

8 

9 

2 

6 

7 

23- 

GRID 

1 

0.0 

2. 

.2 

24- 

GRID 

2 

1. 

2. 

.2 

25- 

GRID 

3 

2. 

2. 

.2 

26- 

GRID 

4 

0.0 

1. 

.1 

27- 

GRID 

5 

0.0 

1. 

.3 

28- 

GRID 

6 

1. 

1. 

.1 

29- 

GRID 

7 

1. 

1. 

.3 

30- 

GRID 

8 

2. 

1. 

.1 

31- 

GRID 

9 

2. 

1. 

.3 

32- 

GRID 

10 

0.0 

0.0 

.0 

33- 

GRID 

11 

0.0 

0.0 

.4 

34- 

GRID 

12 

1. 

0.0 

.0 

35- 

GRID 

13 

1. 

0.0 

.4 

36- 

GRID 

14 

2. 

0.0 

.0 

37- 

GRID 

15 

2. 

0.0 

.4 

38- 

GRID 

101 

0.0 

2. 

.2 

39- 

GRID 

102 

1. 

2. 

.2 

40- 

GRID 

103 

2. 

2. 

.2 

41- 

GRID 

104 

0.0 

1. 

.1 

42- 

GRID 

105 

0.0 

1. 

.3 

43- 

GRID 

106 

1. 

1. 

.1 

44- 

GRID 

107 

1. 

1. 

.3 

45- 

GRID 

108 

2. 

1. 

.1 

46- 

GRID 

109 

2. 

1. 

.3 

47- 

GRID 

110 

0.0 

0.0 

.0 

48- 

GRID 

111 

0.0 

0.0 

.4 

49- 

GRID 

112 

1. 

0.0 

0.0 

50- 

GRID 

113 

1. 

0.0 

.4 
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S 0 

R T E D 

B U 

L K DA 

T A 

ECHO 

CARD 

COUNT 

i 

2 

3 

4 

5 . . 6 

7 ..  8 

51- 

GRID 

114 

2. 

0.0 

0.0 

52- 

GRID 

115 

2. 

0.0 

.4 

53- 

MAT4 

1 

7.6 

54- 

NAT4 

2 

2. 

55- 

PHBDY 

200 

2 

56- 

QHBDY 

10 

AREA4 

1.063 

10 

12 

13 

57- 

QHBDY 

10 

AREA4 

1.063 

12 

14 

15 

58- 

SPC 

7 

101 

70. 

59- 

SPC 

7 

102 

70. 

60- 

SPC 

7 

103 

70. 

61- 

SPC 

7 

104 

70. 

62- 

SPC 

7 

105 

70. 

63- 

SPC 

7 

106 

70. 

64- 

SPC 

7 

107 

70. 

65- 

SPC 

7 

108 

70. 

66- 

SPC 

7 

109 

70. 

67- 

SPC 

7 

110 

70. 

68- 

SPC 

7 

111 

70. 

69- 

SPC 

7 

112 

70. 

70- 

SPC 

7 

113 

70. 

71- 

SPC 

7 

114 

70. 

72- 

SPC 

7 

115 

70. 

73- 

SPC 

70 

10 

1 

400. 

11 

1 

400. 

74- 

SPC 

70 

12 

1 

400. 

13 

1 

400. 

75- 

SPC 

70 

14 

1 

400. 

15 

1 

400. 

76- 

SPCADD 

ENDDATA 

100 

7 

70 

**NO  ERRORS  FOUND  - EXECUTE  NASTRAN  PROGRAM** 


*USER  INFORMATION  MESSAGE  3023, 


***USER  INFORMATION  MESSAGE  3027,  SYMMETRIC  REAL  DECOMPOSITION  TIME  ESTIMATE  IS 
***USER  INFORMATION  MESSAGE  3035 


0 SECONDS. 
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APPENDIX  E 


CH-47C  FORWARD  TRANSMISSION 
BEARING  NUMBERING  SYSTEM  CONVENTION  AND  LOADS 

This  appendix  gives  the  bearing  numbering  system  convention 
and  loads  for  the  CH-47C  forward  transmission.  Figure  E-1 
shows  the  forward  rotor  hub  loads  converted  to  bearing  loads 


Sll 


till 


Figure  E-1. 


Forward  Transmission  CH-47C  Hub  Loads 
Converted  to  Bearing  Loads. 
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Hub  loads  are  converted  to  bearing  loads  and  applied  to  the 
transmission. 

1.  Cover,  Ring  and  Transmission  Case,  Cases  3,  4,  5,  and  6. 

2.  Transmission  Case  Only,  3 Per  Rev  (Cases  1 and  2). 


LOADS  ON  COVER  S-70  SIGN  CONVENTION 


COVER 


© 

ROLLER 

F 

Y 

My 

Case 

3 

1918.5 

0. 

-11.5 

45.584 

1-g 

4 

130420. 

4746.7 

0. 

-35.3 

1158.3 

Symmetric 
Dive  and 
Pullout, 
Noseup 
Pitching 

5 

0. 

-11511. 

0. 

102. 

0. 

Yawing 

6 

-45960. 

-50717. 

0. 

446.7 

-408 . 7 

Recovery 

From  Rolling 
Pullout  (cc) 

-^z 

-F 

Y 

p 

-My 

© 

BALL 

Case 

3 

-7222.7 

-1621.5 

25038. 

-5851.5 

25766. 

4 

-106100. 

-3018.8 

97282. 

-11006. 

381890. 

5 

0. 

22536. 

64855. 

87003. 

0. 

6 

38390. 

41421. 

54066. 

146850. 

-136000. 

203 


SHAFT  SYSTEM  LOADS  ON  CASE 


SUN  GEAR 


ROLLER 

1 

100%  TORQUE ‘ 

-5802.4 

+7082.9 

0. 

-27.757 

-30.377 

Z 

Case  1 3C 

-59.8 

73. 

0. 

-.3 

-.3  i 

4 

2 3S 

-9.9  i 

1 

12. 

0. 

-.05 

-.05  ' 

X 

3 

-4471.3 

5458. 

0. 

21.4 

-23.4 

Y 

V 

/ 

4 

-13055. 

15936. 

0. 

62.4 

-68.3 

1/ 

5 

6 

-10880.  ; 

-8703.6 

i 13280. 

! 10624. 

0.1 

1 1 

1 52.0 

41.6 

-57.  , 

-45.6  i 

SHAFT 

SYSTEM 

-F  i 

Z 

1 

' F 

Y j 

! 

^Z 

-My 

After 

changing 

S-04  signs. 

J 

1 

^ BALL 

(Upper) 

i 

j 

100%  TORQUE  ^ 

+52.220 

+2349.7 

2364.0 

2457.3 

-124.05 

Case  1 3C  , 

.54  1 

24.2 

24.3 

25.3 

-1.3 

2 3S  1 

i .09  j 

1 4.0 

4.02 

4.2 

-.2 

3 1 

! 40.2  i 

1 1 

^ 1810.7 

1821.7 

1893.6 

* -95.6 

4 

i 117.5 

* 5286.8 

5319. 

5528.9 

-279.1 

5 

97.9 

1 

1 4405.7 

4432. 

4607.4 

-232.6 

6 

j 78.3 

3524.5 

3546. 

3686. 

-186.1 

^ BALL 

(Lower) 

1 

» 

100%  TORQUE 

-12.942 

68.971 

53.020 

-67.641 

-15.741 

1 

Case  1 3C 

-.13 

.71 

.65 

-.7 

-.16 

1 

1 

2 3S 

-.02 

.11 

.11 

-.11 

-.03 

1 

3 

-10. 

53.1 

48.6 

-52.1 

-12.13 

4 

-29.1 

155.2 

141.8 

-152.2 

-35.42 

5 

-24.2 

129.3 

118.2 

-126.8 

-29.51 

6 

-19.4 

103.4 

94.5 

-101.5 

-23.6 
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SHAFT  SYSTEM  LOADS  ON  CASE 


PINION 


ROLLER 

F 

X 

F 

Y 

^Z 

My 

100?4  TORQUE 

0. 

-12826. 

1187. 

336.07 

51.847 

Case  1 

0. 

-132. 

12.2 

3.5 

.53 

2 

0. 

-21.8 

2.02 

.6 

.09 

3 

0. 

-9884. 

914.7 

259. 

40. 

4 

0. 

-28860. 

2670.8 

756. 

116.6 

0. 

-24050. 

2225.6 

630. 

97.2 

6 

0. 

-19240. 

1780.5 

504.1 

77.8 

-F 

X 

-F 

Y 

“^Z 

+Mz 

-My 

BALL 

100%  TORQUE 

-61.523 

0. 

0. 

36.355 

10.315 

Case  1 

-.63 

0. 

0. 

.37 

.11 

2 

-.10 

0. 

0. 

.06 

.018 

3 

-47.4 

0. 

0. 

28. 

8.0 

4 

-138.4 

0. 

0. 

81.8 

23.2 

5 

-115.4 

0. 

0. 

68.2 

19.3 

6 

-92.3 

0. 

0. 

54.5 

15.5 

0)  BALL 

100%  TORQUE 

6111.5 

0. 

0. 

1260.0 

359.48 

Case  1 

63. 

0. 

0. 

13. 

3.7 

2 

10.4 

0. 

0. 

2.1 

.61 

3 

4709. 

0. 

0. 

970.9 

277. 

4 

13750. 

0. 

0. 

2835. 

808.8 

5 

11460. 

0. 

0. 

2362. 

674. 

6 

9167. 

0. 

0. 

1890. 

539. 

Z 

A 


X 


SHAFT 

SYSTEM 
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SHAFT  SYSTEM  LOADS  ON  CASE 


PINION 


ROLLER 

F 

X 

F 

Y 

100%  TORQUE 

0. 

3325. 

-2558. 

122.50 

70.313 

Case  1 

0. 

34.2 

-26.3 

1.3 

.724 

2 

0. 

5.6 

-4.3 

.2 

.12 

3 

0. 

2562. 

-1971. 

94.4 

54.2 

4 

0. 

7481. 

-5756. 

275.6 

158.2 

5 

0. 

6234. 

-4796. 

229.7 

131.8 

6 

0. 

4987. 

-3837. 

183.8 

105.4 
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APPENDIX  F 
UPPER  COVER  S-83 

TABLE  F-1.  ORIGINAL  UPPER  COVER  S-83  OUTPUT  INCLUDING 
PRINCIPLE  STRESSES  Original  145.16  Lb 


ELFMFNT 

in 


PRINCIMt  STRESSES 
RAJOR  RtNOR 


S TRAIN  S I T Y 

STRAIN  DENSITY 

NuMaEiS. 


6010 

0.1929363E 

05 

-.67260129 

04 

0.2374239 

02 

6017 

0.1991 4348 

05 

-.33596769 

04 

0.2141359 

02 

6016 

0.1974498E 

05 

-.32298559 

04 

0.2096379 

02 

5020 

0.1775009E 

05 

-.20692709 

04 

0.1630029 

02 

51 

0.11 03 95 7P 

05 

-.10759719 

05 

0.1485159 

02 

53 

0.1371005F 

05 

-.56163599 

04 

0.1274459 

02 

5018 

0.1645783F 

05 

0.65831029 

04 

0.1173739 

02 

5016 

0.1188723E 

05 

-.55250089 

04 

0.1013tl9 

02 

6018 

0.1482290F 

05 

0.54923839 

04 

0.9476269 

01 

5022 

0.2591676F 

04 

-.12974799 

05 

0.9337079 

01 

52 

-.91861339 

03 

-.1425274^ 

05 

0.9323889 

01 

5021 

0.89098639 

04 

-.77393839 

04 

0.8684849 

01 

'21 

0.10384359 

05 

-.59221569 

04 

0.8635369 

01 

6010 

0.25551339 

04 

-.12440559 

05 

0.8626819 

01 

5005 

0.10343709 

05 

-.53407939 

04 

0.8097309 

01 

50 

0.90098249 

04 

-.60594349 

04 

0.7238829 

01 

4017 

-.15412509 

03 

-.11537959 

05 

0.6287459 

01 

5017 

0.11856579 

05 

0. 27846 05^ 

04 

0.6080829 

01 

26 

0.11446719 

C5 

0.56317979 

03 

0.6062639 

01 

6005 

-.40573679 

04 

-.11600009 

05 

0.5790789 

01 

305 

-.55391529 

04 

-.11192569 

05 

0.5581309 

01 

4016 

0.78625399 

03 

-.10457079 

05 

0.5481509 

01 

320 

0.10874799 

05 

0.58041809 

03 

0.5459679 

01 

322 

0. 366591 8c 

04 

-.84740129 

04 

0.4983979 

01 

41 

0.10620329 

05 

0.16997469 

04 

0.4971339 

01 

6011 

-.69681649 

03 

-.1007954E 

05 

0.4652C59 

01 

5011 

0.71535559 

04 

-.47013249 

04 

0.4490259 

01 

56 

0.97661489 

04 

0.10947119 

04 

0.4280679 

01 

55 

0.89837199 

04 

0.63463289 

03 

0.3692699 

01 

25 

0.7C793019 

04 

-.3410eT99 

04 

0.3659159 

01 

58 

0.89667629 

04 

O.I 1161599 

04 

0.3590169 

01 

310 

0.94731729 

04 

-.32587509 

03 

0.3506039 

01 

42 

0.58928639 

04 

-.46093449 

04 

0.3473739 

01 

38 

0.65920129 

04 

-.34574959 

04 

0.3316849 

01 

501  C 

0.67043599 

04 

-.32185009 

04 

0.3275879 

01 

4018 

-.15403179 

04 

-.8464039P 

04 

0.3136399 

01 

317 

0.76505479 

04 

0.1 1213199 

04 

0.2591749 

01 

44 

0.72267669 

04 

0.34757889 

04 

0.2314679 

01 

29 

0.62616059 

04 

-.15162819 

04 

0.2259089 

01 

318 

0.67151459 

C4 

0.46308209 

04 

0.2242969 

01 

401  1 

0.60306569 

04 

-.18332459 

04 

0.2220929 

01 

316 

0.67030049 

04 

0.41997819 

04 

0.2141629 

01 

6015 

0.68359889 

04 

0.2123509' 

04 

0.2007979 

01 

805 

0.17241539 

C4 

-.57393169 

04 

0.2004629 

01 

80  3 

0.65979889 

04 

0.37364159 

04 

0.2004109 

01 

311 

-.34858249 

04 

-.66465909 

04 

0.1992749 

01 

4005 

0.67486689 

04 

0.21738709 

04 

6.1957199 

01 

315 

0.58626919 

04 

0.48623289 

04 

0.1914149 

01 

49 

0.56924069 

04 

-.10811819 

04 

0.1781869 

01 

47 

0.6$2319tt 

^4 

0.63515091 

03 

i 

1 

ml 

6l 

5014 

0.62639499 

04 

0.19025009 

04 

0.1686119 

01 

43 

0.4938375E 

04 

-.2046125E 

04 

0.1661409 

01 

54 

o.4i4Ti9i? 

04 

-.2t48627t 

04 

■ o.lYUdi? 

01 

27 

0.83524129 

03 

-.52966689 

04 

0.1500829 

01 
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TABLE  F-2 


. UPPER  COVER  S-83  RESULTS  WITH  STIFFNESS  ADDED 
TO  EQUALIZE  THE  STRAIN  DENSITY  160.83  Lb 


FLFMFNT 
T_D 


PRINCIPAL  STRESSFS 

NAjriR minor 


STRAI'N  OFNSITy 
STRAIN  density 


" CXeiA-L 


6010 

0.  12361R7E 

05 

-.48465746 

04 

0.1017876 

02 

6018 

0.15223808 

05 

0.3706184F 

04 

0.1001126 

02 

S'' IS 

1460096E 

05 

0.57915516 

04 

0.9232329 

01 

5016 

0. 1393799e 

05 

-.54415706 

04 

0.8878588 

01 

5020 

0. 11437778 

05 

-.39S7061F 

04 

0.8343876 

J1 

6016 

0.  11863508 

05 

-.23926726 

04 

'’.7819458 

''1 

6017 

0. 1247090E 

05 

-.41679306 

03 

0.7568846 

01 

AOlO 

0.3317330F 

04 

-.10551196 

05 

0.6867078 

01 

5022 

0.30358958 

04 

-.10530756 

05 

0.6671178 

01 

6''17 

-.2714''63F 

03 

-.1148384F 

05 

0.6190686 

01 

S21 

0.8583773E 

04 

-.47619498 

04 

0.58049/8 

01 

51 

0.75336138 

04 

-.53507118 

04 

0.5265688 

01 

S'’ '’5 

'’.  72666818 

'’4 

-.48>9516F 

04 

''.4703609 

01 

4016 

0.4418789E 

03 

-.94957626 

04 

0.4427966 

01 

305 

-.49963638 

04 

-.95630236 

04 

0.4121549 

01 

5021 

0.46963448 

''4 

-.64274616 

''4 

''.3917819 

•ll 

5011 

0.67683528 

04 

-.42396296 

04 

0.3891416 

01 

26 

0.920059 88 

04 

0.80694146 

03 

0.3841048 

01 

5017 

0.9166375E 

04 

0.80181256 

'>3 

''.3812946 

“"I 

53. 

0.86896258 

04 

-.53212306 

03 

0.3746796 

01 

320 

0.89983288 

04 

0.54985166 

03 

0.3722856 

01 

32? 

0.31995298 

04 

-.70021026 

04 

0.3488989 

’I 

31'’ 

''.8'“  5656  36 

"4 

-.70333406 

'53 

0.32821 59 

01 

5 010 

0.73395515 

04 

-.20239576 

04 

0.32023/8 

01 

38 

3.6782473F 

04 

-.28829026 

04 

0.3168276 

01 

56 

'’.82295'' '’8 

'>4 

'’.43860726 

'’3 

''.3126726 

01 

52 

0.19837895 

04 

-.72416566 

04 

0.3112176 

01 

44 

0.83757898 

04 

0.39392246 

04 

0.3097626 

01 

41 

0.81558558 

''4 

''.19569956 

'*4 

''.2874878 

Ol 

6011 

-.81838116 

03 

-.79719066 

04 

0.2863986 

01 

6005 

-.3227536F 

04 

-.77545866 

04 

0.2614666 

01 

58 

0.74708838 

04 

0.86601326 

03 

'’.25'' 0976 

'’1 

4'’ 18 

-.  15915118 

04 

-.74670826 

04 

0.2422016 

01 

20 

0.66125868 

04 

-.11796996 

04 

0.2380649 

01 

50 

0.70512388 

04 

0.61092616 

03 

0.2257186 

01 

25 

53527546 

-’4 

-.2916354f 

04 

0.2233986 

01 

317 

0.68912346 

04 

0.76693956 

03 

0.2132106 

01 

4011 

0.56475518 

04 

-.14220886 

04 

0.1854136 

01 

4''05 

'’.6i7''i'’36 

1 

'?4_ 

0.1845346 

01 

311 

-.23977086 

04 

-.65367896 

04 

0.1342046 

01 

42 

0.37588106 

04 

-.38183156 

04 

0.1794216 

01 

6020 

0.54244846 

04 

-.162664'>F 

0.1789806 

01  . 

316 

0. 615220 7F 

04 

0.34357456 

94 

0.1735386 

01 

49 

6015 

0.50149696 

0.5200816F 

04 

04 

-.21295316 

-.14139396 

04 

04  

0.1731406 

0.1602086 

01 

01 

318 

0.5549996F 

04 

0.34023936 

04 

0.1456036 

01 

54 

0.26097356 

04 

-.40963486 

04 

0.1441546 

01 

50 

0.30406076 

04 

-.36227556 

04 

0,1393066 

01 

46 

0.5349S16E 

04 

0.30299346 

03 

9.1319016 

91 

6008 

0.1603876F 

04 

-.43957196 

04 

0.1252436 

01 

805 

0.16666146 

_2t. 

-.43014306 

04 

0.1226696 

91 

315 

'>.48731416 

04 

34556776 

«4 

0.1198296 

01 

5014 

0.50238836 

04 

0.13818426 

04 

0.1086196 

01 

55 

0.34980236 

04 

-.22639016 

04 

TABLE  F-3.  UPPER  COVER  S-83  RESULTS  WITH  STIFFNESS  BOTH 
ADDED  AND  REMOVED  147.80  Lb 

STRAIN  DENSITY 


ELEMENT 

ID 

PRINCIPAI. 

MAJOR 

STRESSES 

MINOR 

STRAIN  DENSITY 

CTRJAZ 

6010 

0.1236197E 

05 

-.4846574E 

04 

0.101787E 

02 

6018 

0.1522380E 

05 

0.3706184E 

04 

0.100112E 

02 

5018 

0.1460096E 

05 

0.5791551E 

04 

0.923232E 

01 

5016 

0.1093799E 

05 

-.5441570E 

04 

0.887858E 

01 

5020 

0.1143777E 

05 

-.3987061E 

04 

0.834387E 

01 

6016 

0.1186350E 

05 

-.2392672E 

04 

0.781945E 

01 

6017 

0.1247090E 

05 

-.4167930E 

03 

0.756884E 

01 

4010 

0.3317330E 

04 

-.1055119E 

05 

0.686707E 

01 

5022 

0.3035895E 

04 

-.1053075E 

05 

0.667117E 

01 

4017 

-.2714063E 

03 

-.1148384E 

05 

0.619068E 

01 

321 

0.8583773E 

04 

-.4761949E 

04 

0.580497E 

01 

51 

0.7533613E 

04 

-.53507115 

04 

0.526568E 

01 

5005 

0.7266691E 

04 

-.4879516E 

04 

0.470360E 

01 

4016 

0.4418789E 

03 

-.9495762E 

04 

0.442796E 

01 

305 

-.4996363F 

04 

-.9563023F 

04 

0.412154E 

01 

5021 

0.4699344E 

04 

-.6427461E 

04 

0.391781E 

01 

6011 

0.6768352E 

04 

-.4239629E 

04 

0.389141E 

01 

26 

0.9200598E 

04 

0.8069414E 

03 

0.384104E 

01 

5017 

0.9166375E 

04 

0.8018125E 

03 

0.381294E 

01 

53 

0.8689625E 

04 

-.5321230E 

03 

0.374679E 

01 

320 

0.8998328E 

04 

0.5498516E 

03 

0.372285E 

01 

322 

0.3199529E 

04 

70021 02E 

04 

0.348898E 

01 

310 

0.8056563E 

04 

'.7003340E 

03 

0.328215E 

01 

5010 

0.7339551E 

04 

-.2023957E 

04 

0.320237E 

01 

38 

0.6782473E 

04 

'.2882902E 

04 

0.316827E 

01 

56 

0.8229500E 

04 

0.4 38607 2E 

03 

0.312672E 

01 

52 

0.1983789E 

04 

-.7241656E 

04 

0.311217E 

01 

44 

0.8375789E 

04 

0.3939224E 

04 

0.309762E 

01 

41 

0.8155855E 

04 

0.1956995E 

04 

0.287487E 

01 

6011 

-.8183811E 

03 

-.7971906E 

04 

0.286398E 

01 

6005 

-.3227536E 

04 

-.7754586E 

04 

0.261466E 

01 

58 

0.7470883E 

04 

0.8660132E 

03 

0.250097E 

01 

4018 

-.1591511E 

04 

-.7467082E 

04 

0.242203E 

01 

29 

0.66125e6E 

04 

-.1179699E 

04 

0.238064E 

01 

59 

0.7051238E 

04 

0.6109241E 

03 

0.225718E 

01 

25 

0.5352754E 

04 

-.2916354E 

04 

0.223398E 

01 

317 

0.6891234E 

04 

0.7669395E 

03 

0.213210E 

01 

4011 

0.5647551E 

04 

-.1422088E 

04 

0.185413E 

01 

4005 

0.6370203E 

04 

0.3452069E 

04 

0.184534E 

01 

311 

-.2397708E 

04 

-.6536789E 

04 

0.184204E 

01 

42 

0.3758810E 

04 

-.3818315E 

04 

0.179421E 

01 

6020 

0.5424484E 

04 

-.1626640E 

04 

0.178980E 

01 

316 

0.6152207E 

04 

0.3435745E 

04 

0.173538E 

01 

49 

0.5014969E 

04 

-.2129531E 

04 

0.173140E 

01 

6015 

0.5200816E 

04 

-.1413939E 

04 

0 . 160208E 

01 

318 

0.5549996E 

04 

0.3402393E 

04 

0.145603E 

01 

54 

0.2609735E 

04 

-.4096348E 

04 

0.144154E 

01 

50 

0.3040607E 

04 

-.3622755E 

04 

0.139306E 

01 

46 

0.5349816E 

04 

0.3029934E 

03 

0.131901E 

01 

6008 

D.1603876E 

04 

-.4395719E 

04 

0.125243E 

01 

805 

0.1666614E 

04 

-.4301430E 

04 

0.122669E 

01 

315 

0.4873141E 

04 

0.3455677E 

04 

0.119829E 

01 

5014 

0.5023883E 

04 

0.1381842E 

04 

0.108619E 

01 

55 

0.3498023E 

04 

-.226390 IE 

04 

0.106242E 

01 
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TABLE  F-4.  UPPER  CX»VER  S-83  RESULTS  WITH  ONLY  STIFFNESS 
REMOVED  132.1  Lb 

STRAIN  DENSITY 

ELEMENT  PRINCIPAL  STRESSES  STRAIN  DENSITY 


ID 

MAJOR 

MINOR 

CTRJAZ 

6010 

0.1956978E 

05 

-.7544402E 

04 

0.253414E 

02 

6017 

0.2018899E 

05 

-.2928082E 

04 

0.215769E 

02 

6016 

0.2003551E 

05 

-.3226973E 

04 

0.215354E 

02 

5020 

0.1853384E 

05 

-.2371160E 

04 

0.179329E 

02 

53 

0. 14747 34E 

05 

-.5956566E 

04 

0.146603E 

02 

51 

0.1116325E 

05 

-.9763926E 

04 

0.137179E 

02 

5018 

0.1719305E 

05 

0.6939863E 

04 

0.128185E 

02 

5022 

0.3829840E 

04 

-.1440469E 

05 

0.122211E 

02 

5016 

0.1252716E 

05 

-.5260559E 

04 

0.107519E 

02 

6018 

0.1539315E 

05 

0.5878719E 

04 

0.102357E 

02 

52 

0.2666102E 

04 

-.1335970E 

05 

0.989767E 

01 

5021 

0. 888663 7E 

04 

-.8663668E 

04 

0.962622E 

01 

4010 

0.3135535E 

04 

-.1291580E 

05 

0.961717E 

01 

321 

0.1084443E 

05 

-.6065684E 

04 

0.930981E 

01 

4017 

-.2627344E 

03 

-.1287309E 

05 

0.779369E 

01 

5005 

0.1037427E 

05 

-.4913027E 

04 

0.779138E 

01 

50 

0.8713559E 

04 

-.5619422E 

04 

0.657653E 

01 

305 

-.5835727E 

04 

-.1184445E 

05 

0.p4505E 

01 

6005 

-.4422484E 

04 

-.1195999E 

05 

0.616865E 

01 

26 

0.1150239E 

05 

0.5356055E 

03 

0.613054E 

01 

320 

0.1148042E 

05 

0.5374883E 

03 

0.610630E 

01 

322 

0.4463922E 

04 

-.8797250E 

04 

0.580296E 

01 

4016 

0.5957070E 

03 

-.1070056E 

05 

0.565999E 

01 

6011 

-.8747656E 

03 

-.1051063E 

05 

0. 50234 3E 

01 

5017 

0.1060345E 

05 

0.1499875E 

04 

0.498775E 

01 

44 

0.1067296E 

05 

0.4420555E 

04 

0.495074E 

01 

5011 

0.7433371E 

04 

-.4606508E 

04 

0.466076E 

01 

5010 

0.8824934E 

04 

-.2439926E 

04 

0.463287E 

01 

41 

0.1012871E 

05 

0.1877117E 

04 

0.4487 2 IE 

01 

56 

0.9982266E 

04 

0.1293656E 

04 

0.444039E 

01 

310 

0.8977410E 

04 

-.9949414E 

03 

0.415077E 

01 

58 

0.9446078E 

04 

0.1058836E 

04 

0.400468E 

01 

38 

0.7407105E 

04 

-.3429916E 

04 

0.392896E 

01 

59 

0.9202398E 

04 

0.5126563E 

03 

0.390469E 

01 

25 

0.6988113E 

04 

338740 2E 

04 

0.357633E 

01 

4018 

-.1930869E 

04 

-.8794379E 

04 

0.335506E 

01 

311 

-.3484292E 

04 

-.8140113E 

04 

0.288929E 

01 

4005 

0.7777238E 

04 

0.4787438E 

04 

0.286354E 

01 

29 

0.7116492E 

04 

-.1420547E 

04 

0 . 280860E 

01 

4011 

0.6875289E 

04 

-.1656828E 

04 

0.272067E 

01 

317 

0.7770340E 

04 

0.1213876E 

04 

0.266460E 

01 

42 

0.5207668E 

04 

-.3975205E 

04 

0.266002E 

01 

316 

0.6974750E 

04 

0.4563543E 

04 

0.236093E 

01 

318 

0.6806883E 

04 

0.4351590E 

04 

0.222652E 

01 

6015 

0.7143563E 

04 

0.1989109E 

04 

0.219553E 

01 

315 

0.6537805E 

04 

0.4678027E 

04 

0.216723E 

01 

805 

0.1809063E 

04 

-.5818906E 

04 

0.208150E 

01 

803 

0.6553496E 

04 

0.3729751E 

04 

0.198012E 

01 

54 

0.4690391E 

04 

-.3151950E 

04 

0.196069E 

01 

49 

0.5907676E 

04 

-.1172797E 

04 

0.193364E 

01 

6020 

0.5779137E 

04 

-.1295853E 

04 

0.189325E 

01 

5014 

0.6533008E 

04 

0.1865928E 

04 

0.183538E 

01 

46 

0.5952234E 

04 

0.5796028E 

03 

0.160042E 

01 

314 

0.5103063E 

04 

-.1462528E 

04 

0.156404E 

01 

table  F-5. 


isa.BBs-.K“  ™ “ 

strain  energies 

irTEMENTSTR" 

o ..TT  TYPES  # 7,82Q7930E  03 

total  for  all  types  » ^ 

percent  of  total 
0.2887 


element-type  t 


element-id 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 


strain -ENERGY 
2.257979E  01 
8.476082E  00 
6.368781E  00 
6.808832E  00 
4. 95967 IE  00 
2.260995E  00 
3.314654E  00 
3.556217E  00 
5.658769E  00 
4.591954E  00 
3.398863E  00 
1.219447E  01 
9.520355E  00 
9.603554E  00 
0.827014E  01 
1.783812E  01 
2.961180E  01 
3.610527E  01 
3.091107E  01 
1.604617E  01 
1.110542E  01 
1.276445E  01 
1.785712E  01 
1.529993E  01 
4.895911E  01 
4.659048E  00 
2.098515E  01 
3.505397E  01 
1.20 587 4E  01 
4.163676E  01 
4.552327E  01 
1,097861E  02 
1.843174E  02 
2.751477E  02 
1.467004E  02 
5.067717E  01 
1. 99997 4E  01 
5.673958E  01 
3.864603E  01 
3.998790E  01 
7.116743E  01 
2.801283E  01 
5.886792E  01 
5.292776E  01 
2.711314E  01 
3.510774E  01 
9.691559E  01 
6.472322E  01 
9.587067E  01 
1.594803E  02 


0.1084 

0.0814 

0.0871 

0.0634 

0.0289 

0.0424 

0.0455 

0.0724 

0.0587 

0.0435 

0.1559 

0.1217 

0.1228 

0.2336 

0.2281 

0.3786 

0.4617 

0.3952 

0.2052 

0.1420 

0.1632 

0.2283 

0.1956 

0.6260 

0.0596 

0.2683 

0.4482 

0.1542 

0.5324 

0.5821 

1.4038 

2.3568 

3.5182 

1.8758 

0.6480 

0.2557 

0.7255 

0.4941 

0.5113 

0.9100 

0.3582 

0.7527 

0.6768 

0.3467 

0.4489 

1.2392 

0.8276 

1.2258 

2.0393 
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APPENDIX  G 


SAMPLE  S-71  RUN 


SAHPLE  PROBLEM  OUASI-ISOTBOPIC  SCIARRA  4/9/76 

**•  INPUT  DATA  *** 

KEYl  ■ 

0 

KEY2  « 

1 

KEY3  - 

0 

KEY4  « 

1 

KEYS  « 

0 

THE  NUMBER 

OF  LAYERS  IN  THB>  LAMINATE  IS  4 

THE  NUMBER  OF  NATERIAL  TYPES  Is  1 

THE  (WBER  Of  MABIHfl  COMMIIbNs  IS  I 

SAMPLE  S-71  RUN  TO  PROVIDE  INPUT 
TO  NASTRAN,  HIGH  MODULUS  GRAPHITE, 
QUASI-ISOTROPIC.  4/13/76 
COST  $1.32 

INCLUDING  TOERMAL  COEFFICIENT  OF 
EXPANSION  GOOD  TO  350»F. 


NATYPE  El  E2  Ul  G ALPHAl  ALPHA2  ALPKA6 

1 0.2500000E  08  O.I700000E  07  0.3000000E  00  0.6500000E  06  -3.3000000E-06  0.1949999E-04  0.0 


•**  LAYER  DATA  •** 


LAYER  NO. 

MATYPE 

ORIENTATION 

THICKNESS 

I 

1 

0.0 

0 . 08000 

2 

1 

45.00000 

0.09000 

3 

1 

-45.00000 

O.ObCCO 

4 

1 

90.00000 

0.0800 

***  ALLOWABLE  STRAIN  DATA  *** 


MATYPE 

LIMIT  STRAIN 

1 > DIRECTION 

LIMIT  STRAIN 

2 * DIRECTION 

LIMIT  STRAIN 

SHEAR 

LIMIT  STRAIN 

1 - DIRECTION 

LIMIT  STRAIN 

2 » DIRECTION 

LIMIT  STRAIN 

SHEAR 

COMPRESSION 

COMPRESS ION 

NEGATIVE 

POSITIVE 

POSITIVE 

POSITIVE 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

070 

***  OUTPUT  DATA  *** 


asHP^nmgspEmis 

A MATRIX 

B MATRIX 

D MATRIX 

0.336e6E  07 

0.10937E  07 

0.25789B-02  -0.22506E  06 

0.0 

-0.37509E  05 

0.34696E  05 

0.33842E  04 

0.38512E-04 

0.10937E  07 

0.3368SE  07 

0.S8628E  00  OTO 

0.22S0£b  06 

-0.1^564fc  oS 

0.33842B  04 

0.34696E  US 

U.H7b5ie-03 

0.58638E  00 

0.113758  07  -0.37509E  05 

-0.37S09E  OS 

0.0 

0.3a512E-04 

0.87551E-02 

0.37579E  04 

<A/T)  MATRIX (A/T)  INVERSE  MATRIX 

0.1^271  08  0.34179E~07  T,a0592E-02  0 . 106iaE-06  -0. 34474B-07  0.175278-13 
Q.341798  07  0.105271  OB  0.18321801  -0.34474E-Q7  Q.106iaa-O6  -0.5464Ba-l3 

O.e0S92»-O2  0.ie321B  01  0,3SM81  07  0.175278-13 


— wma  Lwnmg  bjisto  camnim 

U - 0.941771  07  n - S7  ux  - 66 — QXt  - d.3»4te  67 
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:rmal  expansion  data 


THERMAL  EXPANSION  COEFFICIENT 

X FOR 

COMPOSITE  - 0.1278569E-05 

THERMAL  EXPANSION  COEFFICIENT 

THERMAL  EXPANSION  COEFFICIENT 

Y FOR 

XY  FOR 

COMPOSITE  - 0.1278568E-05 

COMPOSITE  » -0.1340665E-11 

COEFFICIENT  OF  THERMAL 

FORCE 

FORCE 

MX  - 

NY  - 

0.5705626E  01 

0.5705626E  01 

COEFFICIENT  OF  THERMAL 

FORCE 

NXY  « 

-0.7721391E-06 

COEFFICIENT  OF  THERMAL  MOMENT  MX  - 

COEFFICIENT  OF  THERMAL  MOMENT  MY  » 

“"o729?Io49E“00 

-0m2951049E  00 

COEFFICIENT  OF  THERMAL 

MOMENT 

MXY  - 

0.9836811E-01 

***  INPUT  DATA  FOR  COMBINED  N ■ 

M ANALYSIS  *** __ 

LOAD  CASE  NIBSER  1 

NX  - 100. 

MX  « 

0. 

NY  « 0. 

NXY  » Q. 

MY  » 

MXY  - 

0. 

0. 

TEMPERATURE  » 0. 

***  BINDING  OUTPUT  DATA  * 

'** 

A-PRIME  MATRIX 


B'PRIME  MATRIX 


O.025O34OE-O6  -0.3479507E  06  0.2444781E-06  0.5421612E-05  0.1992494E-05  0.4761990E-05 


-0.3479508E-06  0.8250325E-06  -0. 2444779E-06  -0. 1992499E-05  -0.5421597E-05  0.4761975E-05 


0.2444781E-06  -0»2444779E-06  0.1042089E-0S 


71374E-05  0,2471366E-05  0.11569S6E«11 


0.5421612E-05  -0.1992499E-05  0. 2471374E-05  0.6576595E-04  0.9181552E-05  0. 3422757E-04 


0.I992494E-05  -0.5421597E-05 0.2471366E-05 0. 91815S2E-05 0.65765?8E-04  -0. 3422753E-0 

^4?61990E-0S  0,476I975E«05  0.1156956E-11  0. 34227S7E«04  -Q.34227S3E-04  0. 3611685E-03 


C^PRIME  MATRIX 


D>PRIME  MATRIX 


MID-PLAME  STRAINS  AND  CURVAITJRES  *** 


LOAD  CASE  NUMBER  ' 


• 0.82S0339E-04 


• -0.3479507E-04 


0.24447B1E>04 


0.S421611E>03 


0.1992494E-03 


0.4761987E-03 


***  COMBINED  BENDING  AMD 


STRESSES,  STRAINS,  AND  MARGINS  OT  SAFETY  FOR  BACH  LAYER  *** 


Z • -0.160000  THETA  • 0 

-Q.1409E  03  -0.1I62E  03  -0.3363E  02  -0.4242E-0S  -0.6667E-04  -0,S174E-04  -O.IOOOE  01  -O.IOOOE  -1  -O.IOOOE  01 


Z « -0.080000  THETA  > 4S. 

-0.3171E03  -0.4731E  01  -0.S841E  02  -0.1263E-04  0.1022E-05  -0.89B7B-04  -O.IOOOE  01  -O.IOOOE  01  -O.IOOOE  01 


Z ■ 0.0  THETA  • 45. 

0.913SE  03  0.38418  02  -0.7624B  02  0.3608E-04  0.1163B-04  -0.1173E-03  -O.IOOOE  01  -O.IOOOE  01  -O.IOOOE  01 


Z - 0.0  THETA  ■ -45. 

0.31118  03  0.6768E  02  0.76248  02  0.1163B-04 


- 0.080000  THETA  • -45. 

0.15648  03  0.94088  02  0.22248-04  0.84788-04  0.1447B-03  -0.10008  01  -0.10008  01  -0.10008  01 


Z > 0.160000  THBTA  • 90. 

5 0.1352802  0.2880E  03  -0.6542802  -0.29158-05  0. 16928-03  -0.1006B-03  -O.IOOOE  01  -O.IOOOE  01  -0.1000801 


0.1173E-03  -O.IOOOE  01  -O.IOOOE  01  -O.IOOOE  01 
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1547-78 


